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Preface 


This volume of the Transactions of the Society of Rheology 
contains papers which were presented at the Annual Meeting of the 
Society which was held at the Franklin Institute in Philadelphia, 
Pennsylvania, on November 5-7, 1958. For those papers which 
are not presented in full, abstracts are given, and in some cases an 
indication is included of where the complete paper has been, or is 
to be, published. 

The Bingham Medal, which is presented annually in recognition 
of the accomplishments of an outstanding scientist in the field of 
rheology, was presented at the 1958 Meeting to Ronald 8. Rivlin. 
The Bingham Meda! Presentation Address, by E. H. Lee, is given 
at the beginning of this volume. Professor Rivlin’s Medalist Ad- 
dress, “Mathematics and Rheology,’ has been published in the 
May 1959 issue of Physics Today. Former recipients of the Bingham 
Medal, since the inception of the award, are as follows: 


Melvin Mooney 1948 John D. Ferry 1953 
Henry Eyring 1949 Turner Alfrey, Jr. 1954 
W. F. Fair, Jr. 1950 Herbert Leaderman 1955 
Percy W. Bridgman 1951 Arthur V. Tobolsky 1956 
A. Nadai 1952 Clarence M. Zener 1957 


The appearance of the present volume has been delayed to a 
greater than usual extent by an unusually protracted reviewing 
period for the manuscripts. It is planned to decrease this time lag 
considerably in the future volumes, and Volume IV, containing 
papers for the 1959 Thirtieth Anniversary Meeting, under the editor- 
ship of E. H. Lee of Brown University, will appear soon after the 
present volume is issued. 

Officers of the Society, who are elected for a two-year term, are 
for the years 1957-59 as follows: 
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President: Vice-President: 


J. H. Dillon J. H. Elliott 
Textile Research Institute Research Center 
Princeton, N. J. Hercules Powder Co. 
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Presentation of the Bingham Medal to 
Ronald S. Rivlin 


E. H. LEE, Brown University, Providence, Rhode Island 


I would first like to congratulate 
the Committee of Award on a most 
appropriate choice for the Bingham 
Medalist for 1958. I know that you 
will all agree that Professor Ronald 
5. Rivlin’s outstanding contributions 
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Il. Experimental Procedure 


Butt joints of two types were used; they were designed to avoid, as 
far as possible, the complications caused by stress concentration at 
the 3-phase boundary line. In the type shown at the left of Figure 
1, the adherends were steel or glass cylinders, 5 cm. long and 1.9 
em. in diameter. Calculated amounts of polyethylene powder were 
placed on the basis face of one cylinder, three spacers (small pieces of 
metal foil) were put on this face at equal distances from each other, 
and the other cylinder was placed on them, with one of its bases 
facing down. The clearance between the two cylinders is greatly 
exaggerated in the figure. The amount of polyethylene was calculated 
so as to form, on melting, a small drop in the middle of the basis 
faces. As molten polyethylene wetted the solids, the drop acquired 
the shape indicated in the figure. After cooling below the freezing 
range, the adhesive still retained a similar shape. Because its cross 
section was smaller at the ‘‘waistline” than at the adherend-adhesive 
interfaces, stress concentrations obtaining in these interfaces were 
believed to be irrelevant; and indeed the rupture usually occurred 
so far from the interfaces that, after rupture, polyethylene was 
visible with the naked eye on both cylinders. 


, baie | 3 2 
! I 


Fig. 1. Vertical cross section of the adherends used, the clearance between 
them being exaggerated. (1) Adherends. (2) Spacers. (3) The adhesive 
film. 
































atthe 3-phase boundary were believed to be unimportant because of 
thé shape of the adhesive layer inditated in the figure; here again 
the thickness of this layer is greatly exaggerated (for instance, 
a hundredfold). 

The working surfaces of the metal cylinders were ground on 
Carborundum paper and polished with an aqueous suspension of 
alumina until flat within one or a few interference fringes, washed 
with distilled water, and degreased in toluene vapor. The degreasing 
was repeated before every application of the adhesive, but grinding 
and polishing a given cylinder were done again only if the surface 
ceased to be (nearly) optically flat. The glass cylinders were ground 
by the manufacturer. 

The joints were heated in an oven or with a heating tape above 
the melting range of the polyethylene used. They were cooled 
slowly (together with the oven), or at a fair rate (in air) or rapidly 
(in cold water), but apparently the number of experiments was 
insufficient to detect any difference between the results of these 3 
treatments. 

The joints were broken cither in a Suter tensile tester in which the 
load increased approximately 7 X 10° dynes in a second, or in an 
Instron-type machine in which the cross-head advanced either 0.001 
or 0.002 cm. in a second. The breaking stresses at the higher speed 
in this machine were about 5% greater than at the lower speed. The 
data represented in Figure 2 were obtained with the Suter tester, 
and those of Figures 3 and 4 with the Instron-type instrument at 
0.001 cm./sec. 

The specimens serving for the determination of the tensile strength 
of the polyethylenes were either cylindrical rods (for instance, 7 cm. 
long and 0.7 cm. in diameter) with thick ends, or ‘“‘double spades” 
whose straight part had the dimensions of, e.g., 5.3 K 1.3 K 0.3 em. 


Ill. Experimental Results 


No values can be listed for the strength of adhesive joints made 
with commercial polyethylene wax either as received or mixed with 
chlorinated paraffin wax or modified rosin, or vulcanized with sulfur; 
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Fig. 2. Breaking stress of polyethylene joints. Ordinate: Breaking stress. 
Abscissa: Thickness of the adhesive film. The figures at the experimental 
points indicate the number of repeat experiments. The dashed line indicates 
the tensile strength of the material in bulk. Purified Epolene N. Heated to 
130-140°C. Circles: Steel-to-steel joints. Crosses: Glass-to-glass joints. 
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Fig. 3. See legend to Fig. 2. Purified Dylan 3700. Heated to 150-185°C. Cir- 
cles: Steel-to-steel joints. Crosses: Steel-to-glass joints. 
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Fig. 4. See legend to Fig. 2. Purified Marlex 50. Heated to 150—-160°C. ‘Circles: 
Steel-to-glass joints. Crosses: Glass-to-glass joints. 


so many of these joints broke in handling that no average strength 
could be computed. The main data obtained for ‘‘adhesionable’”’ 
polyethylenes are represented in Figures 2, 3, and 4. Their abscissa 
indicates the thickness of the adhesive film, and the breaking stress 
is plotted along the ordinate. The numbers at the experimental 
points denote the number of repeat experiments on which the value 
given is based. The dashed line near the right-hand edge of each 
graph shows the tensile strength of the material, determined on rela- 
tively large samples. 

It is seen that the breaking stress of the joints in all instances is 
greater than the tensile strength of the material in bulk, and that at 
great thicknesses of the ‘glue line’ the former value approaches the 
latter. Thus, as soon as the weak boundary Jayers are removed, 
the strength of an adhesive joint is determined above all by the 
strength of the adhesive in it. The increase of strength in thin layers 
presumably is due chiefly to the smaller probability of bad flaws in 
thinner specimens‘ in the instance of polyethylene wax, Figure 2, 
because this wax is almost a Hookean solid. In the other two 
instances that effect presumably is intensified by the greater difficulty 
of plastic flow in narrower clearances. The two other possible effects, 
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namely, (a) the dependence of texture and of residual stresses on the 
rate of cooling and, thus, on the thickness of the heat-insulating 
adhesive film, and (b) the dependence of the local stresses caused by the 
external force (during the rupture test) on the above thickness, are 
believed to be less important in our experiments. 

In Figures 2 and 3 the circles represent results for steel-to-steel 
joints, while crosses mean glass-to-glass joints in Figure 2 and steel- 
to-glass joints in Figure 3. In Figure 4 circles are for steel-to-glass 
and crosses for glass-to-glass. It may be concluded from inspec- 
tion of the figures that the strength of a joint is independent of the 
nature of the adherend; thus, it is not affected by the molecular 
adhesion between adherend and adhesive; such adhesion obviously 
would depend cn the nature of both of the materials in contact. 


The authors thank Lord Manufacturing Company of Erie, Pennsylvania, and 
Owens-Corning Fiberglas Corporation of Ashton, Rhode Island, for research 
grants which made the above investigation possible. Eastman Chemical Products, 
Inc., Koppers Company, and Shell Chemical Corporation kindly supplied the 
polyethylenes used in this work. 
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Synopsis 


Polyethylenes, from which substances capable of producing weak boundary 
layers were removed, adhered to metals and glasses, and the breaking stress of the 
joints was closely related to the tensile strength of the polyethylenes in bulk, 
independent of the nature of the adherends. Thus, the “rheological” theory of 
adhesive joints was sustained. 


ADDENDUM 


Gladly complying with the suggestion of a referee, we expand in the 
following some of the statements made in the body of the paper. 

1. Experimentally, no significant difference between breaking 
stresses in the arrangement of Figure 1 left and Figure 1 right was 
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observed. This can be seen, for instance, in Figure 4 in which circles 
represent results on systems of ‘“‘Figure 1 right,’’ and crosses, on 
systems of ‘‘Figure 1 left’? type. The agreement between the two 
sets of values indicates that the stress concentrations were similar in 
the two instances. The polyethylene ‘‘pancake’ almost always 
broke relatively far from the two adherends, so that adhesive was 
visible on both rupture surfaces with the unaided eye. It may be 
asked why rupture did not start at the edge of the pin in Figure 1 
right, as there should be a considerable stress concentration at the 
3-phase boundary along this edge. The explanation, in all prob- 
ability, is given by the observation that molten polyethylene rose 
around the pin and, on solidification, formed a thin sleeve covering 
the latter. Thus, the joint was a butt joint below the basis of the 
pin and a tubular lap joint above it. If r is the radius of the pin, 
then as a first approximation the strength of the former would be 
ar’e,, and of the latter, 2xrlc,,! being the height of the pin, o, the 
tensile strength of the adhesive, and oa, its shearing strength. As 
2/1 was twice as great as r, the butt joint would be the weaker of the 
two as long as o, > o,/2. In reality the situation is much more 
complex because the polyethylene sleeve around a steel pin must be in 
a state of considerable radial and circumferential strain; poly- 
ethylenes contract about 12-15 volume-% during freezing while the 
contraction of steel in the same temperature interval is less than 1%. 
The theory of a sleeve shrunk on a rigid cylinder has been worked 
out,*t but its equations cannot be applied in the present instance 
because the mobility of polyethylene in the intermediate stages 
between definitely liquid and definitely solid is not well known. 

2. We attribute the decrease in the strength of polyethylene 
wax joints on increasing the “glue line’ thickness to the greater 
probability of a bad flaw in a thicker specimen. The number of our 
experiments was not sufficient to test this contention by statistical 
analysis, but at least in one instance the prediction of the theory 
seems to be confirmed. The 11 measurements on polyethylene 
wax films between glass cylinders with a film thickness of 13 uy, 
see Figure 2, covered the range of 77.6 to 105.9 bars. If the adhesive 
film had a truly cylindrical shape, then 77.6 bars should be the mean 
breaking stress of joints 13 XK 11 = 143 yu thick (see ref. 4 in the main 

* H. Poritsky, Physics, 5, 406 (1934). 
+t G. D. Redston and J. E. Stanworth, J. Soc. Glass Technol., 29, 48 (1945). 
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body of the paper). The main breaking stress of the 12 joints of 
127 » was 79.3 bars; that is, in accord with the expectation. The 
results on cylinder-pin systems (circles in Fig. 2) are less favorable 
but may be accounted for by the greater divergence of the shape of 
the adhesive film from the theoretical shape. In these systems the 
ranges of the breaking stresses were, e.g., 93.8-129.2 bars for the 33 
tests at a film thickness of 8 uw, and 41.64—-102.0 bars for the 11 
tests at the thickness of 155 u. 

The main argument for the probability theory of the curve ot 
Figure 2 is, however, that the alternative explanations seem to be 
unacceptable. Plastic flow in polyethylene wax may be ruled out; 
the stress-strain curves of this material are almost straight lines, and 
the distance between two marks on a wax rod is practically identical 
before and after extension to rupture. The effect of the film thick- 
ness on rate of cooling and, consequently, on the tensile strength of 
the adhesive cannot be considerable as was demonstrated by the 
experiments mentioned in Section II of the paper. 

On the contrary, the behavior of high-melting polyethylenes (ligs. 
3 and 4) would be expected to be influenced by their plasticity be- 
cause the bars increase in length 2 to 5 times before breaking, and 
only a part of this elongation was reversible. Unfortunately, this 
reversible part has not been measured on the samples referred to in 
the paper, but it was about 75% for thin polyethylene films (about 
0.01-0.03 em. thick); in these films, 100 units of length were ex- 
tended to 120 in the moment of rupture and became 105 after the rup- 
ture. Presumably, this elasticity during flow is responsible for the 
fact that the curves of Figures 3 and 4 are so different from those pre- 
dicted by the theory of squeezing a perfectly plastic drop between two 
rigid plates. This theory, usually attributed to Prandtl (1923) and 
reproduced in textbooks,* predicts the breaking stress (as long as the 
radius of the joint is constant and much greater than the thickness of 
the adhesive film) to be inversely proportional to this thickness. 
The actual effect is definitely smaller. For instance, in Figure 3, 
the 7 experiments at 23 « gave rupture stresses of 197.2 to 222.0 bars, 
and the 7 experiments at 114 y, 130.3 to 151.7 bars. Prandtl’s 
equations would give 39 to 44 bars for the strength of the thicker 
joints. 

* F. K. Th. van Iterson, Plasticity in Engineering. Blackie, London, Glasgow, 
1947, p. 88. 
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The Theory of Locking Materials* 


ARIS PHILLIPS, Department of Civil Engineering, Yale University, 
New Hawen, Connecticut 


1. Introduction 


According to the nomenclature introduced by Prager,' materials 
may be classified into soft and hard. For a soft material the stress 
increment required to produce a specified strain increment diminishes 
with increasing deformation. For a hard material this stress incre- 
ment grows as the deformation proceeds. Examples of the two types 
of materials are aluminum, which is a soft material, and cellular rub- 
ber, which is a hard material. 


A 8 
% 














Figure 1. 


An extreme case of a soft material is the rigid-plastic one, Figure 1 
(line A). For such a material the stress increases up to the yield 
stress o) without any strain, and then the stress stays at o» while the 
strain can increase to any value. An extreme case of a hard material 


* The results presented in this paper were obtained in the course of research 
sponsored by the Office of Naval Research. 
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is the ideal locking one, Figure 1 (line B). For such a material the 
strain increases without any stress up to the locking strain «, and 
then the strain locks at « and cannot be increased while the stress 
can be increased to any value. 

In this paper we shall consider possible stress-strain laws of ideal 
locking materials, and with the help of an example discuss a peculiar- 
ity of the concept of locking. The lack of experimental results for 
combined stresses concerning the deformational response of materials 
which in a simple tension or simple compression test exhibit the 
stress-strain curve of a hard material, obliges us to be guided in this 
investigation by theoretical considerations only. 

Prager,' in his discussion of the ideal locking materials, points out 
that theoretical speculations in plasticity about the structure of the 
yield condition or stress-strain law repeatedly guided investigators to 
unexpected experimental results and that the same may be expected 
to happen in the field of locking materials. The present writer con- 
curs with the ideas expressed above, and he also believes that an 
investigation of the peculiarities of the locking concept will help in 
improving our understanding of the structure of the stress-strain laws 
of real materials. 


Il. Stress-Strain Relations 


Very general stress-strain relations for an ideal locking material 
can be found by introducing a locking condition. This is a relation 
between the six strain components which must be satisfied at a point 
of a body when locking exists at this point. A locking condition can 
be represented by a locking surface in the six-dimensional strain space. 
Then locking occurs at a given point of the body when and only when 
the endpoint of the vector @€1, €22, €s3, ~/2e2, ~/ Zee, ~/ es), repre- 
senting the state of strain, lies on the locking surface. 

In addition to the locking condition there must be available a lock- 
ing rule. Consider the six-dimensional stress space in which the stress 
at a given point of a body is represented by the vector #(ou, o22, 033, 
V 201, V20n, ~/20n). By superimposing the stress space on the 
strain space we can express the locking rule as a relation between the 
directions of the two vectors zand &. In Prager’s formulation of the 
theory of locking materials it is assumed that the vector @ is normal 
to the plane which is tangent to the locking surface at the endpoint 
of z, Figure 2. Such a locking rule can be expressed analytically by 
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Gi = A(0g/0e4;) (1) 


where 


g(éu, €22, €33, €12, €23, én) = k (2) 


is the locking condition. In eq. (1) the factor \ is a proportionality 
factor. Equations (1) and (2) together with the relation o,; = 0 for 
g < k are a complete system of stress-strain relations for the ideal 
locking material. 











%y 


Figure 2. 


In this paper we shall assume that there are two independent rheo- 
logical equations of state, one for the change in volume and the other 
for the change in shape. This assumption is introduced because it is 
generally accepted that materials behave differently when changing 
shape than when changing volume.? Of course the assumption of the 
independence of the two rheological equations of state is strictly valid 
only for small strains.* This assumption of two independent rheo- 
logical equations of state has been used by Prager in a recent paper.‘ 
He introduced a material which follows Hooke’s law for changes in 
shape while it has a volumetric equation of state of the elastic-locking 
type, Figure 3. In the present paper we shall assume that both 
equations of state are of the ideal locking type, Figure 4, and we wish 
to explore the behavior resulting from such a choice. 

The volumetric equation of state for our material is 


tn =0 for  |en| < |knl (3) 
om = indefinite for ¢, = kp 
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where ¢,, and ¢, are the mean stress and mean strain, respectively, 
and k,, is the value which «,, must reach in order for the volumetric 
locking to occur. Obviously | ¢, | is always smaller or equal to | k,, |. 


om 











Figure 3. 


on 




















(0) 
Figure 4. 


To write the distortional equation of state we must introduce the 
distortional locking condition 


glen, €22, €33, €12, €23, a) =k (4) 


where ¢,, are the six components of the strain deviation, g is a function 
of these components, and k is the value which this function g must 
take for distortional locking to occur. When g < k no deviatoric 
stresses are possible; for g = k deviatoric stresses are possible; and 
finally g can never take values greater than k. From the locking con- 
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dition it is apparent that at locking the six components of strain de- 
viation do not need to remain constant; it is only required that any 
change in the values of these components will be such that the func- 
tion g will retain the value k. 

To introduce the locking rule we proceed as follows: 


From the stress-strain diagram B in Figure 1 it follows that at the 
locking stage, 


ode = 0 


Generalizing for combined stresses we can write, using the summation 
convention, 


8 de, = 0 (5) 


where deé,,; is the increment of the component e,,; of the strain deviation 
and s, is the corresponding component of the stress deviation. 
Equation (5) says that the distortional strain energy per unit volume 
cannot change at the locking stage. . 

Because at locking the function g has a constant value, it follows 
that 


dg = (dg/de,,)de,, = 0 (6) 


Hence, from eqs. (5) and (6) it follows that the vector with the com- 
ponents s, is parallel to the vector with the components g/0e,). 
Therefore we can write 


8iy = A(Og/ Oe) (7) 


where } is a positive factor of proportionality which as a rule depends 
on the coordinates x,y,z. We conclude that eqs. (4) and (7) represent 
the equation of state for change in shape. 

The argument could be advanced that in analogy with the theory 
of rigid plastic bodies we should use the incremental stress-strain 
relations 


ds = A(Og/de45) (8) 


instead of the relations (7). That such stress-strain relations as 
given by eq. (8) are not acceptable follows frem the following discus- 
sion. Let us consider eq. (7) first. In Figure 5 we give a schematic 
two-dimensional representation of the locking condition g = k. When 
distortional locking occurs, the endpoint A of the vector é@ with the 
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Figure 5. 


components ¢,; is on the surface g. At the same time the vector AB 
with the components 0g/de,; is normal to the surface at A, and the 
vector § with the components s, is parallel to AB. A change ds 
in the stress will produce a new stress vector § + ds, and the corre- 
sponding change in the strain will be found by obtaining a new normal 
CD to the surface g = k. This normal CD is parallel to the vector 
3 + ds. Then the change de in the strain due to the change ds in 
stress will be given by the vector AC. Any arbitrary infinitesimal 
change ds in the stress is generally permitted because, in general, we 
should always be able to find a normal CD to the surface parallel to 
the stress vector. 

Consider now the incremental stress-strain relations (8). In this 
case ds, and not §, is parallel to the normal AB. Therefore, the stress 
increment is fixed in direction, independently of the direction of 
strain increment de. It follows that there is an infinite number of 
strain increments corresponding to the same stress increment; that 
is, we have a condition of instability. If such an instability should 
be avoided stress-strain relations of the type (8) cannot be used. 

Another remark which should be made is that the locking surface 
should have no portion which is plane. Indeed, if a finite portion of 
the locking surface is plane, then it is easy to see that an infinitesimal 
increment in the stress deviator may produce finite increments in the 
strain deviator, a situation which is not permissible in stress-strain 
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ds 
, 


sds) |F 


= 


Figure 6. 








relations. From Figure 6 it is seen that to the strain deviator vector 
é there corresponds the stress deviator vector 5 normal to the plane 
portion of the locking surface. A small change ds of § will change 
the direction of the stress deviator vector from that of 5 to that of § 
+ ds, so that the endpoint of @ cannot any longer be on the same 
plane portion of the locking surface. Hence de must be of finite length 
whereas ds is infinitesimal. 


Ill. The Basic Problem of Locking 


Comparing the stress-strain relations of the ideal locking material 
with those of the rigid plastic one of the total strain theory, it is con- 
cluded that the ideal locking material does not follow from the rigid 
plastic one by a simple interchange of the stresses and strains. Of 
course, the distortional equation of state of the ideal locking material 
follows from that of the rigid plastic one by the simple interchange of 
stresses and strains. The same procedure, however, is not valid for 
the volumetric equation of state because ¢, = 0 for |on| < |km| and 
ém = indefinite for ¢,, = k,, is not the volumetric equation of state of 
the rigid plastic material (unless, of course, we assume |k,,| = ©). 
This observation and the remark that the differential equations of 
equilibrium for the stresses and the compatibility equations for the 
strains do not have the same form, show that the solution of problems 
with ideal locking materials cannot be achieved by a formal repetition 
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of the procedure used for rigid plastic materials and an interchange of 
the stresses with the strains. 

The basic problem of the ideal locking theory is to find the values 
of the boundary displacements for which the body will lock so that an 
increase in the external loads will produce no additional displace- 
ments. These values of the boundary displacements are the locking 
displacements. They correspond to the collapse load of the rigid 
plastic theory. The problem of determining the locking displace- 
ments can be solved either by means of the theorems introduced by 
Prager' or more directly by integration of the appropriate differential 
equations. This second method will be illustrated in the next sec- 
tion. 

To the locking displacements there corresponds a locking mode in 
the body. Part of the body locks volumetrically and part of it locks 
distortionally. In general, volumetric and distortional locking may 
coexist over the same region of the body. In some cases, however, 
the two types of locking cannot coexist in the same region. Indeed, 
suppose that we have a problem of spherical or axial symmetr~ in 
which v = w = 0 so that the only nonvanishing displacement com- 
ponent is u and this component depends on the coordinate r only. 
This is, for example, the case in the problems of the hollow sphere, 
and of the long hollow cylinder. In such a case the two equations 
ém = k, and g = k give two differential equations for u(r) which have 
different solutions. It follows that it is impossible for both solutions 
to be valid over the same region of the body and therefore in this case 
of locking the body will be divided into two parts, one of which locks 
volumetrically while the other locks distortionally. The calculation 
of the locking displacements and of the locking mode will be given in 
the next section for the problem of the hollow sphere. 


IV. The Hollow Sphere 
(a) Strains and Displacements 


We consider a hollow sphere of inner radius a and outer radius b. 
The displacements of the inner and outer surface are denoted by u, 
and wu, respectively. We shall assume that the ratio u,/m = s is 
fixed. We want to find the value of u, for which locking occurs. 

Because of spherical symmetry we have ¢, = « and the strain- 
displacement relations are 
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& = u/r, e- = du/dr (9) 


Volumetric locking will occur at a point when at this point 


€- + 2e, = 3km (10) 
which, with the help of eqs. (9), gives the differential equation 
(du/dr) + 2(u/r) — 3k, = 0 (11) 


The solution of this equation is 
u = (C/r?) + kar (12) 


Distortional locking will occur at a point when at this point g = k. 
As function g we shall select the octahedral shearing strain yp so that 
the distortional locking condition has the form 


2 /s 
I= n= 5 | ( — @)* + (g — «)* + «© — a) | =k (13) 


Because of the existing spherical symmetry this condition becomes 
&-~-¢ = 3k/2+/2 = ky (14) 
where k, is a new material constant. With the help of eqs. (9), eq. 
(14) becomes 
(du/dr) — (u/r) +k, =0 (15) 
the solution of which is 
u=Cyr—krinr (16) 
Obviously, expressions (12) and (16) cannot both be valid over the 
same region simultaneously. It is, however, possible that one ex- 
pression is valid over part of the hollow sphere, while the other is valid 
over the rest of the sphere. Let us assume that expression (16) is 
valid between r = a andr = ¢, and that expression (12) is valid be- 
tween r = tandr = b. Then we have 
U, = Ca — kalna 


from which C, can be calculated, so that finally we find 
u = (u,/a)r + kyr In (a/r) (17) 


fora Sr St. Similarly we find 


up = (C/b?) + kab 
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from which C can be calculated, and thus we find 
u = up(b?/r?) — k,, b?/r? + kar 


fort srs b. 
The strain components are as follows: 
é& = (u,/a) + hk, In (a/r), €- = (u,/a) + ki ln (a/r) — ky (19) 
forasrst, 
fo = Uq(b?/r*) — k,,(b?/r*) + ky 
(20) 
€, = —2u)(b?/r*) + 2k,,(b?/r*) + k, 


fort srs b. 
The following boundary conditions are now available for the de- 
termination of the quantities u,, u, and ¢. Atr = ¢ we must have: 


€, + 26, = 3k,, (i) 
when using expressions (19); 
&—¢= k (ii) 


when using expressions (20); and (iii) the same values for u when 


using eqs. (17) and (18). The first condition gives 


Ua = —k,a In (a/t) + k)(a/3) + ak» (21) 
the second condition gives 
Uy = Ky (t*/3b*) + kmb (22) 


and finally the third condition is identically satisfied. Equations 
(21), (22), and relation u,/u = 8 determine the locking displacements 
Ue, Ua, and the corresponding value of ¢ (radius of the boundary 
separating the two locking regions). 

It can be shown quite easily that the assumption a S r S ¢ for 
distortional locking and t S r S b for volumetric locking is the only 
possible one. We have to show that the region of distortional locking, 
if existing, must be bounded by the inner boundary; and that the 
region of volumetric locking, if existing, must be bounded by the 
outer boundary. Indeed, let us consider the elastic-locking material 
remembering that the ideal locking material is an idealization of the 
elastic-locking material; that is, of a material which is first elastic 
and then it locks. 
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If we consider the solution of the elastic hollow sphere under in- 
ternal and external pressure we find (see ref. 5); 


€ — ¢, = 3B/r* (23) 
Sen = € + Ze = 3A (24) 


where B and A are constants independent of r but dependent on the 
boundary conditions. From eq. (24) we see that volumetric locking 
must occur over the entire portion of the sphere which has not al- 
ready locked distortionally as soon as 3¢,, becomes equal to 3k». 
Similarly, eq. (23) shows that distortional locking must start at the 
inner boundary because there, forr = a,¢, — ¢,isamaximum. Hence, 
as the loading proceeds, either (i) ¢, — ¢«, becomes equal to k; at r = 
a first, or (ii) ¢, becomes equal to k,, over the entire hollow sphere 
first. 

In case (ii), volumetric locking occurs over the entire sphere and 
no distortional locking is possible. In case (i), distortional locking 
will occur at the inner boundary, and as the loading proceeds, the 
region of distortional locking will extend in size until either the entire 
sphere will become distortionally locked, or until for some value of t 
the portion of the sphere which has not locked distortionally will lock 
volumetrically (because ¢,, finally becomes equal to k,, at this region). 

Therefore, reverting again to the ideal-locking material, we con- 
clude that distortional locking, if it will appear at all, will appear at 
the inner boundary, whereas volumetric locking, if it wil! appear at 
all, will appear at the outer boundary. 


(b) Stresses 


In the region of distortional locking, a S r S t, the following rela- 
tions are valid: 


gn = 0 (25) 
because of no volumetric locking, and 
Sy = A(Og/Oe,;) (26) 


because of distortional locking. The octahedral shearing strain can 
be written also in the form 


2 ; _ 
7 2 V3 Ve? + 9" + Coe” ak (27) 
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7 a 
Oe 4 V3 V e,,? + és’ + ¢,,* 
from which it follows that 
sy = (V/2d/ki ey 





(29) 


From eq. (25) it follows that eq. (29) can be written in the form 


0, = (V2/ki)(e — &m), 9 = Fy = (V/2A/ks)(€) — &m) (30) 


from which we find 
, — 09 = (V2A/ki)(e, — &) 
On the other hand 
30m, = 0, + 2o, = 0 


Hence, eqs. (31) and (32) give 


a, = (24/2/3)(A/ki)(e — 9) 
which with the help of eq. (14) becomes 
o, = —(24/2/3)r 
Then eq. (32) gives 
oy = (V/2/3)r 
The equation of equilibrium is 
(do,/dr) + 2(¢, — o9)/r = 0 
which, because of eqs. (34) and (35) becomes 
(dd\/dr) + (3A/r) = 0 
Integrating we find 
A= C/r* 
from which it follows that 
a, = —(20/2/3)(C/r*), 0% = (02/3)(C/r*) 
for the region a Sr St. 


(31) 


(32) 


(33) 


(34) 


(35) 


(36) 


(37) 


(38) 
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In the region of volumetric locking, t S r S b, we have the following 
relations: 
& = 0 (40) 
because of no distortional locking, and 
om = indefinite (41) 
because of volumetric locking. Equation (40) gives 
O, = 09 = F, = Om (42) 
and equation of equilibrium (36) becomes 
de,/dr = 0 (43) 
from which we find 
o, = C; (44) 


At r = t the value of o, from eq. (44) must be the same as the 
value of ¢, from eq. (39). Hence 


—(24/2/3)(C/t®) = Ci (45) 
Recapitulating we write 
asrst:o, = —(2v/2/3)(C/r), oy = (v/2/3)(C/r*) (46) 
b<r=b: o, = —(2v/2/3)(C/t), 0 = —(20/2/3)(C/#) 


It is seen that there is always a stress discontinuity for o, at r = t. 


(c) Discussion 


If o,_ is the radial stress at r = a and c, is the radial stress at r = 
b, then from eqs. (46) we find 


Cra/Trm = t*/a* (47) 


As t has already been determined from the displacements, eq. (47) 
determines the ratio of the inner to outer pressure corresponding to 
the locking displacements. We observe that the values of ¢,, and 
7» are not determined; only their ratio is fixed. We remark also 
that ¢,, and ¢,, must always have the same sign because ¢*/a* is always 
positive. 





26 A. PHILLIPS 


From eq. (47) we also see that if the ratio ¢,,/¢,) is known, then ¢ is 
also known. The value of ¢ together with the value of uv,/u and 
eqs. (21) and (22) gives the values of the locking displacements. 

From eqs. (21), (22), and (47) it is seen that the radius ¢ of the sur- 
face separating the two locking regions is a function of the ratios 
Ugi/ Up OF Ora/e. By changing either of these two ratios we can 
move the surface separating the two locking regions, so that an ele- 
ment of the hollow sphere may change from a condition of distortional 
locking to a condition of volumetric locking and vice versa. Such a 
change will be followed, of course, by a change in the stress and strain 
distribution. 

Thus, it is seen that at locking there exist no corresponding single 
set of values (u,:, %:) but rather a continuous range of sets of values 
(tai, Vox), With each value of the ratio u,,/u; corresponding to a defi- 
nite value of the radius ¢. Once the hollow sphere has locked for a 
given value of ¢ and for the corresponding values of u,; and u,, we 
can change the set (tq:, Up, t) to a new set (u’a:, uy, t’) satisfying 
eqs. (21) and (22), and the hollow sphere will still remain locked. As 
t cannot become smaller than a or larger than b, it follows that there 
are limits imposed on the selection of the locking displacements. 
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Synopsis 

The concept of the ideal locking material, which has been introduced by Prager, 
is considered for the case of materials which have two equations of state, one for 
the change in volume and the other for the change in shape. Possible stress- 
strain relations are discussed. The concepts of volumetric locking and of distor- 
tional locking are introduced, and it is shown by means of an example that there 
are cases in which these two modes of locking are incompatible over the same 
region of the body. The fundamental differences between locking materials 


and ideal plastic materials are pointed out. The problem of locking of a hollow 
sphere is discussed in detail. 
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An Apparatus for Measuring the Dynamic 
Mechanical Properties of High Polymers 


at Low Frequencies 


R. H. SHOULBERG, F. H. ZIMMERLI, and O. C. KOHLER, 
Rohm and Haas Co., Philadelphia, Pennsylvania 


1. Introduction 


The wide interest currently shown in the dynamic mechanical 
properties of polymers, both in bulk and in solution, obviates the 
need for any lengthy introduction to measurements of this type. 
Several authors have treated the theory and methods for deter- 
mining these properties; excellent reviews and extensive bibli- 
ographies can be found in the recently published Rheology.' 

Morrisson, Zapas, and De Witt? have described a torsional pen- 
dulum apparatus for measuring the response of a sample to the 
application of a stress whose magnitude varies sinusoidally with 
time. This apparatus has virtues in simplicity of construction and 
operation, in its capability for handling both solid and liquid samples, 
and in its capabilities for scanning intensively the range of frequency 
between 10~‘ and 10 cycles/sec. 

The apparatus described in this report was patterned after that of 
Morrisson, Zapas, and De Witt, but has extensive modifications. 
In the development of the apparatus, certain difficulties were en- 
countered which should be. of interest to other investigators of 
dynamic mechanical properties. 

The two principal objectives achieved by the modifications have 
been the introduction of continuous automatic recording, on a single 
chart, of the applied stress and the sample strain; and the introduc- 
tion of a constant impedance network which, in effect, allows accurate 
measurements of stress and strain over a large range of specimen 
rigidity without the need to alter the specimen geometry. In addi- 
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tion to these features, the apparatus has been made quite compact, 
rapid changes in the temperature of the sample and its surroundings 
can be effected, and the stress can be easily changed both in amplitude 
and frequency. 

The apparatus has been used to obtain data on a variety of solid 
samples; its use with liquid samples is just beginning. Some detailed 
results on a sample of poly(isopropyl methacrylate), in bulk, will be 
presented. 


II. Description of Apparatus 


The apparatus, shown in Figure 1, combines the basic design of 
Morrisson, Zapas, and De Witt? with the unique recording ability 
of a German Lumiscript RLt4 recorder manufactured by Hartmann 
and Braun. The torsion pendulum system is basically the same as 
previously described.? The recorder was nicely adaptable for the 
measurement of the sample strain because visible reflections from a 
high pressure mercury are lamp provide an immediately available, 
yet permanent, record on sensitized chart paper. Normally the 
image is positioned by a system of mirrors mounted on galva- 
nometers within the recorder. Since it is the oscillation of one of these 
mirrors which causes the image to sweep back and forth, the oscillation 
of the strained sample can be recorded if the sample is fitted with a 
properly designed mirror. Figure 2, a photograph of the sample 
suspension, shows the result of combining these features. The mirror 
for recording the strain is mounted on the upper edge of the torque 
producing coil, which is rigidly attached to the upper edge of the 
sample. 

Calculations of dynamic properties are made from the trace of the 
torsion pendulum and from the trace of a recorder-mounted galva- 
nometer which measures the torque-producing current passing through 
the coil. These traces are recorded simultaneously on a common 
time base on a single chart since the light beams can overlap without 
interference. The two traces need not have the same axis of sym- 
metry, although by rotation of either the equilibrium position of 
the sample or of one or more of the mirrors, they can be made to 
have coincident axes. The traces should, of course, be coincident 
on the time axis, which means they should lie on the same horizontal 
line. This can easily be done by raising or lowering the entire sample 
suspension. Differences in phase and amplitude as frequency and/or 
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Drive Coil 


Sample 





Fig. 2. Enlarged view of specimen suspension. 


temperature are changed can thus be easily measured. Recently, 
the authors found that Lethersich* used a similar scheme. He 
employed a drum camera for recording, whereas in the present 
apparatus the lack of a need for film development has simplified the 
collection of data. 

The sample suspension is comprised of a rectangular self-supporting 
coil of 120 turns of 32 gage copper wire hung in the field of a per- 
manent magnetron magnet. A nonlinear magnetic field obtained 
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by making the pole pieces flat and the inner core round exactly com- 
pensates for the nonlinearity obtained in recording a rotating beam 
on a flat chart. The coil is suspended from a 37 gage copper wire 
about 3'/, in. long; electrical lead-in ribbons to the coil are pieces of 
0.0007 in. gold foil about '/1. in. wide. 

As noted, the mirror is attached to the upper end of the coil; 
the lower end of the coil supports a stainless steel tube which in 
turn holds the upper sample clamp. The lower sample clamp is 
guided in a support which allows vertical freedom of motion but 
restricts rotation. This support, which weighs about 20 g., can be 
buoyed up in a pool of mercury if necessary, so that tensile stresses 
in the specimen will have a negligible effect on measured values of the 
shear modulus. (A treatment of this effect is given in detail by 
Timoshenko.*) 

The sinusoidal driving torque is generated in a sine-wave potenti- 
ometer (Precision Potentiometer RVP3, Technology Instrument Co., 
Acton, Mass.) which is driven at 16 different speeds between 0.0004 
and 1 cycles/sec. by means of either a 3 r.p.m. or a 60 r.p.m. Barcol 
synchronous motor driving an Insco geared speed reducer. The 
output shaft of the reducer also drives the recording film in the 
recorder via a flexible shaft; therefore, all records have the same 
cycle length—about 3 in. 

Current applied to the torque-producing coil is supplied to the drive 
coil via an intermediate amplifier. Currents as high as 12 ma. can 
be applied to the coil. The amplifier is designed to load the 20,000 
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Fig. 3. Resistance network for apparatus. 
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ohm sine wave potentiometer negligibly and to provide the output 
current from a high impedance source. The amplifier feeds a resist- 
ance network which is shown in Figure 3. By suitable adjustments 
of R, and R,, according to a schedule, the drive coil is made to ‘“‘see”’ 
an almost exactly constant resistance. (As a first approximation it 
is clear that if the sum of R, + R, is made constant in all measure- 
ments, the drive coil will be in series with an effectively constant 
resistance because of the high resistance of the amplifier.) The 
current galvanometer coil, again because of the high resistance 
of the amplifier, always sees an external resistance of about 1250 Q. 
Electromagnetic damping in the circuit is thus held constant. With 
this design each signal can be adjusted to produce oscillating am- 
plitudes approximately equal to the recording film half-width, 
although the specimen rigidity may vary by a factor of a thousand. 

The sample is surrounded by an oven which has been designed to 
combine adequate control of the sample temperature, rapid change of 
temperature when desired, and modest space requirements. It is 
constructed of two sections of Transite pipe arranged concentrically, 
with the space between them filled with glass wool. The sample 
chamber is heated by an auxiliary heater (25 2) and a control heater 
(125 2), each of which is composed of four sections distributed sym- 
metrically about the axis of the oven and shielded from the sample 
and thermocouple by anodized aluminum semicylinders. The 
auxiliary heater is controlled by a thermostat (300°C. maximum) 
which is set slightly below the desired temperature. The control 
heater is operated from a Duenna controller (manufactured by 
Instruments for Research and Industry, Cheltenham, Pa.). It is 
actuated from an iron—Constantan thermocouple placed in the 
sample chamber. A circular brass chamber through which cooling 
water can be circulated can be mounted above the oven to prevent the 
driving coil and suspension wire from being heated when the oven is 
operated at high temperatures. 


Ill. Treatment of Data 


The treatment of data generally follows that of Morrisson et al., 
with some modifications; repetition of the information given therein 
will be avoided, wherever possible. 

The basic equation of motion for the torsional pendulum is 
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d*9 n’ dé G’ ; 
(5) + (+2) 24 (9 + ha) 0 = kde (1) 


where J is the inertia of the suspension plus the anchored specimen, 
6 is the angle of rotation of the suspension, ¢ is the time, 7’ is the 
dynamic viscosity of the specimen, A is a geometrical factor, B 
is an instrument damping coefficient, G’ is the dynamic shear modulus, 
ko is an instrument stiffness coefficient, k, is the torque produced by 
unit current through the drive coil, and 7,, is the instantaneous cur- 
rent through the drive coil. 

The inertia of a thin rectangular sample with one end anchored 
is given by 


I, = me*/36 (2) 


where m is the mass of the sample and c is the large dimension of the 
specimen’s cross section. In the apparatus described herein the 
inertia of the suspension is 18.5 g.-cm.*; J, is usually about 0.03 
g.-cm’ and can be neglected. 

The damping coefficient of the instrument, B, includes air damping, 
eddy current damping within the drive coil wires themselves as they 


move in the magnetic field, and resistance damping caused by cur- 
rents set up in the external circuit by the back EMF produced by 
the coil moving out of phase with the driving current. The former 
effects will probably be negligibly small; the latter effect can be 
significant and variable unless care is taken to minimize it, as can 
be done by having the coils “‘see” a large constant external resist- 
ance. (The frequencies involved are such that all impedances are 
almost purely resistive.) In the present instrument B is about 1 
dyne-cm. /rad./sec. 

The stiffness coefficient of the instrument, ko, includes not only the 
effects of the suspension wire, but also the effects of any magnetic 
impurities in the copper wiring of the drive coil, since these will, of 
course, cause the coil to seek some equilibrium position in the field. 
In this apparatus magnetic impurities of the order of 1 part in 50,000 
are sufficient to make ky twice as great as would be the case in their 
absence. (Pugh* has recently presented some data on magnetic 
effects in 99.999% pure copper.) These impurities cause no difficulty 
if they are accounted for, but since even magnetic dust particles on 
the coil become important in a sensitive movement, ky should be 
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determined frequently. In this apparatus ky is about 250 dyne-cm./ 
rad. 

For a sinusoidal input current to the drive coil it can be shown that 
the dynamic shear modulus, G’, and the dynamic viscosity, 7’, of the 
sample are given by: 


G’ = A E (‘*)(‘) cos 6; = ko — le] (3) 
Veg, 0 
n =A k (‘*)(‘*) sin 6; — B| (4) 
WwW beg, 0 


where ig, is the current amplitude in the drive coil, i,, is the current 
amplitude in the current galvanometer, 4 is the amplitude of angular 
motion of the specimen (a maximum of 4° in this apparatus), 5, is 
the lag angle between input and output, and w is the circular fre- 
quency. In calibrating the apparatus, J, ko, and B can easily be 
found by well-known methods. The group k¢(tcg,/0o)(tae/teg) 
includes the current-torque coefficient k, and the quantities %,,,, 
ic, and %. If a6 v. dry cell in series with a 1000 © resistance is used 
as a current source in place of the amplifier and the traces of @ and 
icg are equalized in a steady-state measurement in the absence of a 
sample by adjusting the various resistances of the impedance net- 
work, then igc,/tep, can be calculated and a value of k; can be found 
which, although not truly the torque produced by unit current 
through the drive coil, is invariant as long as the magnetic field 
strength and the optical system remain unchanged. In such a 
measurement the electromagnetic torque is balanced by the torque 
resulting from twisting the suspension coil, so that 


ko = K(Geq,/ Oo) / (Gaee/ teas) (5) 


No dynamic effects are involved and a precise measurement of k, 
can easily be made. 


IV. Experimental Results 


A number of bulk polymers have been investigated in the apparatus. 
For purposes of demonstrating the type of data which has been ob- 
tained, the results of tests on three materials will be reported. These 
comprise limited measurements on steel and polyisobutylene and 
more extensive findings on poly(isopropyl methacrylate). Before 
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presenting these data, however, several things should be mentioned 
concerning the accuracy of the data and the selection of sample 
geometries. 

The ratio of the amplitudes of the two traces can usually be deter- 
mined within 1%, unless the sample is so stiff that its amplitude 
cannot be made roughly equal to the film half-width. The phase 
angle can be determined within about 1 degree. It will be recognized 
from an inspection of eqs. (3) and (4) that the accuracy of calculations 
of G’ and n’ will vary widely, depending upon the specimen, upon the 
care used to evaluate ko, k,, J, and upon the accuracy of the compo- 
nents of the impedance network. It is desirable to make the cross 
section of the specimen (relatively) large to avoid excessive tensile 
stresses; however, the dimensions should also be chosen so that the 
specimen is neither so stiff at low temperature that its twist amplitude 
cannot be made large enough for accurate measurement, nor so 
soft at high temperature that most of the rigidity of the system results 
from ko. It is also essential, especially with soft materials, to make 
the ratio of sample length to sample width as large as possible in 
order to reduce the importance of any stresses caused by clamping; 
this point will be discussed in detail shortly. It may not be possible 
to satisfy all these requirements with a single specimen over large 
frequency and temperature ranges. With Plexiglas, for example, it 
has been desirable to introduce one change of sample geometry be- 
tween 26 and 155°C.; however, for many samples a single specimen 
is sufficient. The specimens used to date have varied from 2 to 90 
mils in thickness, from '/, to '/, in. in width, and from 1 to 3 in. in 
length. 

Cold-rolled shim steel (0.0021 in. thick) was tested at room tem- 
perature. The real part of the complex shear modulus was found to 
be 11,200,000 psi at both 0.5 and 0.008 cycles/sec., in good agree- 
ment with values of the (static) shear modulus given in various engi- 
neering handbooks. No phase lag could be detected, indicating 
very low internal damping. 

Measurements on a piece of poly(isobutylene) obtained through the 
courtesy of Dr. Lawrence A. Wood of the National Bureau of Stand- 
ards gave, in cur initial calculations, values of G’ at 26°C. which 
were higher by an almost constant factor than those reported by 
Philippoff* on material from the same batch. The values of the loss 
factor, on the other hand, were in good agreement. This suggested 
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either an error in the measurement of the specimen’s dimensions or an 
error caused by the sample holder constraining the movement of 
the cross sections at both ends of the specimens. The latter effect 
has apparently not been discussed in connection with the measure- 
ment of dynamic mechanical properties, and in this particular case 
was the major cause of discrepancies between the two sets of data. 
For this reason it will be discussed in greater detail. Timoshenko*® 
has shown that if a rectangular sample is twisted by movements which 
in some fashion restrain the warping of the cross sections, the stiff- 
ness of the piece is increased, essentially because the normal stresses 
which restrain the warping also contribute to the strain energy of the 
bar when it is twisted. If, for example, a narrow rectangular speci- 
men is twisted by a torque M at each end, with a torque 2M applied 
at the mid-section in the opposite sense, the mid-section, because of 
symmetry, cannot warp. The effective length of the specimen, 


l’, is given by 
5(1 + ») *) 
’ => SC 
Vv =] (1 r ; (6) 


where / is the actual length, v is Poisson’s ratio, and c is the width of 
the piece; c is assumed to be much greater than the thickness. 
Obviously, for long narrow specimens this effect is small. Since 
l/c for our poly(isobutylene) was only about 2.5, it seemed desirable 
to determine whether any analogous effects were important. We 
first examined a specimen of low-density polyethylene, on which we 
were working, and found no significant effects.. We later tested a 
piece of the NBS poly(isobutylene) and indeed found a significant 
effect. In this latter test we varied //c and calculated apparent 
values of G as a function of l/c at fixed frequency (1 cycle/sec.) 
and temperature. If we assume that 


Vv’ = Kl — ye/l) (7) 
where y is a constant, then for any two tests 


L (1 — pe/h) Gel, 
ly’ ~ ip (l— We/h) Gils 8) 


If we normalize a series of measurements with respect to test 1, then 


Gi/Ga = (1 aad ¥ c/l,)/K (9) 
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where 


K=1- ¥ c/l, (10) 


The results of our tests are shown in Figure 4, from which it appears 
that K = 0.91 and y =~ 0.45. Considering that the sample was 
unclamped and reclamped after each change in length, it is gratifying 
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Fig. 4. Dependence of apparent shear rigidity upon specimen geometry for poly- 
(isobutylene). 
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Fig. 5. Real part of complex shear modulus of poly(isopropyl methacrylate). 
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Fig. 6. Dynamic mechanical properties of poly(isopropyl methacrylate). 


that the data fall so well on a straight line. The smaller effect of end 
clamping with polyethylene is probably due to the fact that, with a 
given clamping technique, the amount of constraint of the cross 
section is smaller the stiffer the material. Therefore, this source 
of error is most important with the softest materials. After correct- 
ing for this effect and making more precise measurements of the 
sample thickness, we found substantial agreement with the results of 
Philippoff. 

Poly(isopropyl methacrylate), an amorphous bulk polymer whose 
glass temperature (81°C.) makes it quite convenient to study in 
both the rubbery and glassy states in this machine, was studied over 
fairly wide temperature and frequency ranges. The particular 
sample tested had a weight-average molecular weight of 320,000. 
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The real part of the complex shear modulus is shown at several tem- 
peratures in Figure 5. Figures 6a, b, and c show the real and imag- 
inary parts of the reduced shear modulus and the loss factor as 
functions of the reduced frequency, a7w. The real part of the 
complex modulus was given greater weight in determining the ar’s 
because of the higher accuracy of most of the data. By trial and 
error it was determined that if the characteristic temperature, 7’,, 
of Williams’ is chosen as 120°C., then good agreement is shown in 
the range of validity between values of ar as predicted by the Wil- 
liams-Landel-Ferry equation 


—8.86 (T — T,) 


] = 
6 or 1016+ T —T, 


(11) 


and the values of ar found by matching the G’ curves at different 
temperatures; a comparison of such values is given in Table I. 


TABLE I 
Comparison of ay Values for Poly(isopropyl Methacrylate) 





Log ar from Log a7 as used 
, ae eq. (11) in Fig. 6 





127 0 (0) 
124 - 0.27 0. 
115 18 

104 2.59 

95 

87 . 
68 2. 8. 
27 14. 





One expects increasing discrepancies at temperatures more than 
about 50°C. below 7’,, and this indeed is found to be the case. How- 
ever, the generally good agreement of the a7’s with those of eq. (11), 
as well as the other data mentioned, appears to demonstrate that the 
apparatus can be used to investigate materials in its design range 
without gross errors. 


The authors wish to thank Dr. H. Markovitz and Mr. L. J. Zapas for their aid 
to the authors, Dr. B. S. Garrett for his continued interest, Mr. M. 8. Bechta for 
aid in the mechanical design of the equipment, Mr. E. A. Taylor for aid in the 
amplifier design, and Mr. J. A. Shetter and Mrs. E. H. Cherry for collection of 
data presented herein. 
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Synopsis 


An apparatus is described for low-frequency measurements (0.0004 to 1 cycles/ 
sec.) of the dynamic mechanical properties of polymers. The apparatus is basi- 
cally a modified version of the Morrisson-Zapas-De Witt machine at The Mellon 
Institute. The modifications consist principally in (1) using reflected light beams 
from a mercury are source to trace both input and response curves on a photo- 
sensitive film, thus providing a common time base for recording without resorting 
to the use of any intermediate electrical system, and (2) using a constant imped- 
ance network to independently vary the recorded amplitude of the input and out- 


put traces, thus allowing accurate measurements of amplitude ratios over the 
entire frequency and temperature ranges (25-250°C.) without altering the 
geometry of the specimen. The calibration of the instrument and data treat- 
ment are discussed, and examples of the measurements made on solid samples 
are given. 
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Stress Analysis for Linear Viscoelastic Materials* 


E. H. LEE, J. R. M. RADOK,7{ and W. B. WOODWARD, t Division 
of Applied Mathematics, Brown University, Providence, Rhode Island 


1. Introduction 


In recent years methods have been developed for the determination 
of stress distributions in linear viscoelastic materials. Such solutions 
display features in marked contrast to the analogous distributions 
for elastic bodies, and these differences can have important signifi- 
cance technologically in the application of polymers. The purpose 
of this paper is to show that a variety of types of mathematical 
problems arise in determining stress distributions in linear visco- 
elastic bodies, and to present examples of two types of solution with 
comments on features which can be of significance in design. It is 
felt worthwhile to present such examples in detail, since the earlier 
literature has not discussed specific practical examples with considera- 
tion of design features. 

The general problem of stress analysis is to determine the stress 
components o;;(z;,t) throughout a body, such as that shown sche- 
matically in Figure 1, when the loading and constraint of the surface 
are prescribed, as well as the body forces f,(z,,t). Cartesian 
coordinates xz, are used, with i taking on the values 1, 2, and 3; 
the usual vector and tensor suffix notation are adopted. Thus, 
oy(%x,t) signifies the variation of the stress components on, on, ..., 
as functions of position in the body (2;, x2, 73) and the time t. Time 
arises in an essential way because of the viscoelasticity of the materials 


* Sponsored under Department of Defense Contracts, NORD 11496 and 16471 
between the Bureau of Ordnance, Navy Department, and Brown University. 

+ Present address: Department of Aeronautical Engineering and Applied 
Mechanics, Brooklyn Polytechnic Institute, Brooklyn, New York. 
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Fig. 1. Boundary conditions for the body V. 


under consideration, and the boundary data prescribed are in general 
functions of time, as are the resulting dependent variables. This 
applies even for the quasi-static situation, for which inertia forces 
are negligible compared with applied forces. Figure 1 indicates that 
the boundary traction 7’; is prescribed over the part S, of the surface 
S, and the displacement u, over the part S,, and that the regions S, 
and S; may vary with time. Other mixed combinations of traction 
components and displacement components can also be prescribed as 
in elasticity theory; for example: prescribed shear traction and 
normal displacement. 

The stress components o,; must be chosen to satisfy the constitutive 
equations for the linear viscoelastic material under consideration. 
The analysis is here limited to infinitesimal strain, so that the strain 
components are represented in terms of the displacement vector 
according to: 


jy = 1/2[(Ou,/dax;) + (Ou;/Ozx;) | (1.1) 


For linear isotropic viscoelasticity, two independent pairs of visco- 
elastic operators are needed, analogously to the two independent 
elastic constants of isotropic elasticity. Choosing these to relate 
shear and dilatational behavior respectively, and writing s,, and ey; 
for the deviator or shear components of the stress and strain tensors, 
respectively, the constitutive equations take the form: 
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P(8y) = Q(es) (1.2) 
P' (ou) == Q’ (ex) (1.3) 


where the operator pairs P,Q and P’,Q’ can take any of the various 
forms, differential, integral, or complex algebraic, by which linear 
viscoelastic relations can be represented,' and where the usual repeated 
suffix summation convention is utilized: oy = on + o22 + «3. 

The equations of motion must also be satisfied: 


(Ooy/Ox;) + fi = p(O?u,/Ot?) (1.4) 


where p is the density and f; the body force per unit volume. 

In addition to the field equations, compatibility of the stress and 
displacement variables with the prescribed boundary conditions 
requires that: 


T; = gyn; ON Si (1.5) 
u, prescribed on S, (1.6) 


where n, is the outward normal unit vector. 

The eqs. (1.1) to (1.4) form a system of linear partial differential 
equations in the independent variables x, and t. Boundary condi- 
tions are given by eqs. (1.5) and (1.6), and in general initial conditions 
over the body V at zero time will be required. Here consideration is 
limited to initially undisturbed bodies. 

For the specia) circumstances that the form of the body V does not 
vary with time, apart from the infinitesimal displacements associated 
with the strain, and that the regions S,; and S, over which load and 
displacement conditions are prescribed also are invariant with respect 
to time, application of the Laplace transform? removes the time from 
the equations and boundary conditions and reduces the viscoelastic 
problem to an equivalent elastic problem, known as the associated 
elastic problem.* This then brings to bear on the problem the 
extensively developed methods of the theory of elasticity, and Lap- 
lace inversion of any particular stress component so evaluated 
provides the stress distribution for the original viscoelastic problem. 
An example of this approach is given in Section 2 of this paper. 

If the form of the body varies with time, such as the annihilation of 
material in a burning out rocket grain, this method does not apply 
since the Laplace transform is an integral over time from zero to 
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infinity, which is not defined for the part of the material which only 
exists for a finite time. Nor does the method apply if S,; and S, 
vary with time, for then at a particular surface point the prescription 
of surface traction 7, for example, will apply only up to a certain 
time t; when the boundary between S, and S; traverses the point, and 
thereafter u, will be prescribed. -Thus the Laplace transform, being 
an integral from zero to infinity in time, cannot be determined for 
either traction or displacement, one of which would be needed to 
specify the associated elastic problem. 

When the Laplace transform method cannot be used, a tentative 
method of approach is available in terms of the sequence of elastic 
problems for the instantaneous boundary values at each time 1, 
but a check by an alternative approach is called for. This can take 
the form of application of the Laplace transform in an inverse way 
when enough additional information has been deduced to make its 
application possible. Such an example is the solution of the contact 
problem for viscoelastic bodies.® 

At the present time no general method of approach appears to be 
available when the Laplace transform is not applicable, but specialized 
problems can be treated by particular methods. An example is 
given in Section 3 of this paper by the determination of stresses in a 
thick walled cylinder of viscoelastic material subjected to internal 
pressure with annihilation of material associated with a continuously 
increasing inner radius. 


2. Stresses in an Externally Reinforced Viscoelastic Cylinder 


As an example which can be treated by application of the Laplace 
transform, consider the thick walled viscoelastic cylinder shown in 


Fig. 2. Reinforced viscoelastic cylinder subjected to internal pressure. 
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Figure 2, reinforced by an external elastic sheet and subjected to 
internal pressure. The cylinder is assumed to be long, and the 
deformation plane strain with no axial motion. Inertia forces will be 
neglected, which will be permissible if the load durations considered 
are long compared with the time for stress wave propagation across 
the cylinder. The resulting quasi-static analysis is obtained by 
substituting p = 0 in eq. (1.4). Denoting the Laplace transform of a 
function f(t) by f(p), where p is the transform parameter,’ eq. (1.1) 
becomes : ; 


&1j(Le,P) = 1/2[(0%,/dzx;) + (O%,/O2x,) | (2.1) 


and eq. (1.4), for the case needed here of zero body forces and negli- 
gible inertia forces, becomes: 


08,,/0x, = 0 (2.2) 

The boundary conditions take the form: 
T, = an; on S, (2.3) 
a, = prescribed on S, (2.4) 


The constitutive eqs. (1.2) and (1.3) are modified by application of 
the Laplace transform, since the operators become algebraic multi- 
pliers which are functions of the transform parameter p: 


P3y = Qty (2.5) 
Pty = 0’ ex (2.6) 


For differential operators, which are associated with the common 
finite spring-dashpot models for viscoelastic materials, the 0/0t 
in the viscoelastic operators is simply replaced by p in the transformed 
operators, for initially undisturbed conditions.’ 

The system of equations and boundary conditions (2.1)—(2.6) 
constitutes the associated elastic problem for the determination of 
3:j(%x,p). The usual elastic constants corresponding to eqs. (2.5) 
and (2.6) are the shear modulus: 


G = Q/2P (2.7) 


and the bulk modulus: 


K = 0'/3P’ (2.8) 
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The other corresponding elastic constants can be deduced in the 
usual way from these (ref. 6, p. 10), so that solutions in terms of, 
say, Young’s modulus and Poisson’s ratio can be utilized if this 
representation is more convenient. Bars will be retained over the 
elastic constants deduced from viscoelastic operators as a reminder 
of their significance. They represent transforms of the corresponding 
viscoelastic operators, for example, Young’s modulus E£ is the quotient 
of the transformed operator pair which relates simple longitudinal 
stress and the corresponding strain for the viscoelastic body. 

The particular problem of the thick cylinder under consideration 
can be solved directly in cylindrical polar coordinates (r,6,z) in terms 
of the well known general Lamé solution (ref. 6, p. 56): 


3, = A’ — (B’/r’) (2.9) 
Seg = A’ + (B’/r’) (2.10) 
where A’ and B’ are constants in the elastic case or arbitrary functions 
of the transform parameter for the associated elastic problem. They 
are determined by the boundary conditions at the inner and outer 


radii, r = a, r = b, respectively. The internal pressure II(i) is 
prescribed, and for this tractton eq. (2.3) takes the form 


Il(p) -_ Grr, r=¢a (2.11) 


The boundary condition at r = b must express the equilibrium be- 
tween the stress in the reinforcement and the radial pressure on the 
cylinder, and the equality of the circumferential strain of the cylinder 
and the reinforcement in contact with it. The reinforcement is 
assumed to be thin so that membrane theory can be used with the 
neglect of bending influences. If F is the transformed circumferential 


tension in the reinforcement, then radial equilibrium at r = b demands 
that: 


Fd = —#,,b dd (2.12) 


Using Hooke’s law for the circumferential strain in the case of zero 
axial displacement (ref. 6, p. 23), and the suffix R for the elastic 
constants of the elastic reinforcing material of thickness h, the 
equality of the circumferential strain requires that: 


P (1 — vp*) _ aee(1 — 9) — WL + Par 
, ee E 





ig(b,p) = (2.13) 
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where E is Young’s modulus and » Poisson’s ratio. Note that 
because the reinforcement is elastic, Zz and vz are constants, and not 
functions of the parameter p. Elimination of F between eqs. (2.12) 
and (2.13) leads to the boundary condition: 





. P a(l — 3’) b 214 
aes ~ ’ = 2. 
“1+? —E 86 





where the constant a = Eph/(1 — vg")b contains the properties of 
the reinforcement. This is a more involved type of boundary condi- 
tion than prescribed traction or displacement, but since it is the 
transform of a boundary condition in the viscoelastic case which does 
not change in form throughout the loading duration, it falls within the 
scope of the transform method of solution. 

Substituting the boundary conditions (2.11) and (2.14) into the 
general solution, eqs. (2.9) and (2.10) permit A’ and B’ to one evalu- 
ated, giving the required particular solution: 














et int (~ + 1) ~ (= - i) 
nore Ff Fle 
r= a(l — p*) (= ) (" ) sities 
aa) fees Saath 
ai(1 + ») — E \a’? a? 
a(l — ») (= i 1) _ (= + 1) 
a aW(1 + 9) — E \r? a (2.16) 





a(l — p*) b? b? 
-~ Ban + 3) a= ~ I 
aW(l + 9) — E \a’ a’ 
In order to evaluate the stress distribution for a particular example, 


the internal pressure is taken to be suddenly applied at ¢ = 0, and 
maintained constant thereafter of magnitude Ip, which determines :* 


0 = IIo/p S (2.17) 


The operators for a particular viscoelastic material for the cylinder 
must be inserted, and a material elastic in dilatation and behaving as a 
Maxwell body in shear is assumed, thus giving for the special cases 
of eqs. (1.2) and (1.3): 

[A(0/dt) + Blsy = (0/Odt)ex (2.18) 


(2.19) 





= Cew 
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where A, B, and C are material constants. The shear behavior 
(2.18) corresponds to the Maxwell viscoelastic model shown in 
Figure 3. The corresponding transformed relations, equivalent to 
eqs. (2.5) and (2.6) are: 











(Ap B)31; = per (2.20) 
Gu = Ci, (2.21) 


The equivalent G and K according to eqs (2.7) and (2.8) determine 
the corresponding E and » (ref. 6, p. 10) in the form: 


RK 3Cp ‘ (AC — 1)p + BC 
y= 


~ (2AC + 1)p + 2BC (2AC + 1)p + 2BC 


and substitution into eqs. (2.15) and (2.16) determines the transforms 
of the stresses for the particular problem considered: 


r b? b? 7 
(1 == = aA )( 240 + 1) — 3aA 4 p* 
r 4 
b? b2 
rn + (1 _ “) 2BC — aB (s "2 +4AC + 1) p — 2aB°C 
0 


b? 2 
i 3 ales )—-aatle 
a a 
b? b? 
+ [—“ 2BC — aB\3 a 4AC + 1} |p — 2aB°C 
a a a . a 


(2.22) 


(1 + pr = aA )(240 + 1) + 3aA “| p* 
r r 


b?\, b? \ 
Il +i)\{1+ 2BC — aB —3 + 4AC+1 p — 2aB°C 
0 
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Fig. 3. Maxwell model for the assumed viscoelastic properties in shear. 


Such a form of transform, consisting of the ratio of polynominals in 
p, can readily be inverted by the method of partial fractions’ to 
give the variation of the stress distribution in the viscoelastic case. 
This procedure involves the evaluation of the roots of the quadratic 
denominator, which would lead to cumbersome expressions in the 
general algebraic form of eqs. (2.22) and (2.23), so that particular 
material constants have been substituted to permit numerical 
evaluation. 

Typical values have been selecved. The bulk modulus K = C/3 
is taken to be the same as for water: 10° p.s.i. Young’s modulus for 
instantaneous elastic response, which with K determines the constant 
A, is taken to be FE = 10° p.s.i., and this corresponds to A = ‘/; X 
10-*. The only time which appears in the problem is the relaxation 
time of the Maxwell behavior, r, = A/B time units, and so for con- 
venience time measurement is taken in terms of this unit, so that 
™ = land A = B. b/ais taken to be 2, and h/b to be '/3s (this 
corresponds roughly to h = 1/,;” for b = 4”). The reinforcement is 
assumed to be steel with Ey = 3 X 10’ p.s.i., ve = 1/0/11, which 
determines a = 10°. 

With these numerical values, 


Il, [333p* + 656p + 320] + (b*/r®) [147p* + 144p] 
p 921p? + 1232p + 320 





o,, = -— 


for example, which in terms of partial fractions can be written 


0 


II 
or, — = = | (.ss16 a 
Pp 





.005282 .2233 ) 
p+ 9849 p+ .3528 


2 


b 
+ = (1.506 - 
: 





001320 _.05583 )] 
p+ 9849 p+ .3528 
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This can be inverted directly to give: 


005282 (muon) 4. 2238 
9849 3528 


001320 05583 |, _ wy] 


Or = —Ily | {2616 + 


(1 = e~-wayt 


b2 

— .1596 _ 1 =» —.9849%\, 2 
va 1 geag “~~ “B508 
with similar relations for og. These functions are plotted in Figures 
4 and 5 to give the stress distributions as functions of the radius at 
several values of time. 
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Fig. 4. Variation of the distribution of radial compressive stress. 





It is seen that although the internal pressure I) remains constant, 
the stress distribution for the viscoelastic cylinder varies very markedly 
with time. This is in sharp contrast with the elastic situation for 
which an invariant stress distribution would be obtained under 
constant pressure. In general the viscoelastic nature of a material 
causes variations of the stress distribution with constant surface 
loading. 

Figure 4 shows that the radial compressive stress at the outside, 
r = b, increases gradually from about half the internal pressure to the 
full internal pressure for long durations of loading compared with the 
relaxation time of the material in shear. This is associated with the 
relaxation of the shear strength of the cylinder material, according to 
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its Maxwell behavior (Fig. 3), while it is constrained by the re- 
inforeement. The shear strength relaxes to zero, and the visco- 
elastic material behaves as a liquid under uniform hydrostatic pressure 
so that the full internal pressure is transmitted to the reinforcement. 
The circumferential stress is then also uniform hydrostatic pressure, 
of magnitude IIp. Since the tension in the reinforcement F is directly 
proportional to the externa) radial compressive stress in the cylinder 
according to eq. (2.12), F increases gradually on application of the 
constant internal pressure. Thus initially the viscoelastic cylinder 
supports part of the internal pressure, but this is gradually trans- 
ferred to the steel reinforcement. 


10+ 


——_ 





Fig. 5. Variation of the distribution of circumferential compressive stress. 


Figure 5 shows that initially hoop tension occurs adjacent to the 
bore of the cylinder. The magnitude, and even the sign of this, 
depends on the stiffness of the reinforcement and the radius ratio 
b/a. This circumferential tension changes to compression as the 
pressure is maintained, and finally the limit of uniform hydrostatic 
compression is reached when the shear strength has relaxed to zero. 

These features can have an important influence on the design of 
such a cylinder. If, for example, a pulse of pressure is to be applied 
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for a limited time, the delay in the buildup of tension in the re- 
inforcement may be utilized to ease the design problem for the re- 
inforcement by choosing the relaxation time of the cylinder material 
to prevent full transmission of the internal pressure before it decreases. 
The linearity of the problem permits the analysis for a pulse of stress 
to be obtained from the present solution by superposition. 

The change to biaxial compression caused by the viscoelasticity of 
the cylinder material may be useful in preventing the initiation of 
fracture cracks. The initial hoop tensile stress could be reduced by 
applying the internal pressure gradually, and this again is a feature of 
viscoelastic response which provides freedom to meet design criteria 
which is not available in the case of elastic materials. It is the 
additional degree of freedom associated with the variation in time 
which permits this. 


3. Stresses in a Cylinder with Annihilating Inner Radius 


As an example of a viscoelastic stress analysis problem which is not 
amenable to treatment through application of the Laplace transform, 
the cylinder problem similar to that treated above, but with an 
annihilating inner radius: r = a(t), da/di > 0, is considered in this 
section. The problem is to be treated within the confines of in- 
finitesimal displacement theory, and the motion of the inner boundary 
r = a(t) is not to be confused with the particle motion there. At any 
time ¢, the material which initially corresponded to r < a(t) has been 
annihilated, for example, burned and removed as gas, and the particle 
displacement at the boundary is u(a(t,),t;) which is considered 
infinitesimal in magnitude and negligible compared with the radius 
a(t,). As pointed out in Section 1 of this paper, the Laplace trans- 
form cannot be used since material at radius r; only exists for a finite 
time t, given by r, = a(t), so that the transform of stress or dis- 
placement at this radius is not defined. No general methods of 
attack on this kind of problem are at present available, so that a 
specialized procedure is sought. 

In order that the direct approach detailed below can be applied, 
the problem is restricted to the situation, which is common for 
viscoelastic materials, in which shear strains dominate the deforma- 
tion, so that dilatational strains can be neglected and the theory of 
incompressible motion applied. For the plane strain situation with 
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no axial motion, incompressibility prescribes the radial displacement 
to be of the form: 


u(r,t) = k(t)/r (3.1) 


since with no volume change the outflow of material across each radius, 
2xru, must be the same for all radii at a given instant ¢. k(t) is an 
arbitrary function. Equation (3.1) corresponds to the strain com- 
ponents: 


é = u/r = k(i)/r’ (3.2) 
é- = Ou/Or = —k(t)/r’? (3.3) 


The viscoelastic relation (1.2) can be conveniently written in terms 
of stress differences: 


Play — 7) = Qles — &) = Q(2K(t)/r") (3.4) 


For the coordinate system (r, 6, z) all shear components of stress and 
strain are zero by symmetry, so that the single suffix notation for 
components of stress and strain is satisfactory. In cylindrical 
coordinates, the equilibrium equation (1.4) without body force and 
inertia terms is: 


r(00,/Or) = a, — a; (3.5) 
so that o, can be eliminated between eqs. (3.4) and (3.5) to give: 
P(r(O0,/Or)) = Q(2k(t)/r?) (3.6) 


The time operators and space derivative can be permuted since r 
and ¢ are independent variables, and integration of eq. (3.6) with 
respect to r gives: 


P(o,) = —Q(k(®))/r? + (3.7) 


where k,(¢) is another arbitrary function of time. Equations (3.5) 
and (3.7) then give: 


P(e) = Q(k())/r*? + nO (3.8) 


This general solution is analogous to the transform solution (2.9), 
(2.10), but contains the viscoelastic operators P and Q. As before 
the two arbitrary coefficients are determined by the boundary con- 
ditions at r = a(t), andr = b. 
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The equilibrium condition for the reinforcement (2.12) applies for 
nontransformed variables, and the circumferential strain compati- 
bility can be written down directly froin (3.2): 


k(t) _ F (1 — vp’) 


~ int one” (3.9) 


Combination of these leads to the boundary condition: 


k(t) b (1 — vp’) 

— = ee Cy =b ; 

7 “ae Wiha damn (3.10) 
Elimination of o, between eqs. (3.7) and (3.10) determines k;,(?) 
in terms of k(t): 


1 /P 
k(t) = — ie (7 - @) (k(t)) (3.11) 
where C = (b/h)(1 — v,*/Eg) is a constant expressing the properties 
of the reinforcement. This provides a new relation for the radial 
stress by substituting eq. (3.11) in eq. (3.7): 


k 
Qk) (3.12) 


From this k(t) must be determined to satisfy the boundary condition 
at the surface of the cavity: 


o, = —Il(#), r = a(t) (3.13) 


In general, this combination of relations to determine k(t) is not easy 
to treat, since p(c,) in eq. (3.12) is a time operator for fixed radius, 
and eq. (3.13) gives only the instantaneous value o, as the moving 
boundary passes through the radius a(t). Thus the operator P 
cannot be applied directly to II(¢). 





| 
P—— 1+ @, + 


Fig. 6. The Kelvin or Voigt model for assumed viscoelastic properties in shear. 
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This difficulty does not arise in the case of the Kelvin delayed 
elastic material corresponding to the viscoelastic model shown in 
Figure 6, for then: 


P=1, Q=4H(0/dt) +J (3.14) 


where H and J are material constants. This case is considered below, 
and substitution of —II(¢) on the left-hand side of eq. (3.12) for 
r = a(t) gives a differential equation for the determination of k(t): 


1 1 \ dk 1 1 1 
3 (aor 4 5) at (J lar - =| + aa) k = M(t) (3.15) 


k(t) is obtained from eq. (3.15) with the initial condition k = 0, 
since for the Kelvin material there is no instantaneous strain when 
the internal pressure is first applied. The relations detailed above 
then determine the entire stress and strain history for the viscoelastic 
cylinder. 

For many choices of a(t) it will be necessary to soive eq. (3.15) 
numerically, but a solution in terms of tabulated functions can be 
obtained by taking: 


a(t) = ao/(1 — xt)” (3.16) 


which approximates constant velocity annihilation for r/aj) < 2. 
The corresponding displacement vs. time plots for various ratios b/d, 
are shown in Figure 7, where the time is measured as a proportion of 
the annihilation time & for which a(t) = b. It is seen that for 
b?/a,.? = 2, a significant approximation to constant velocity annihila- 
tion is achieved by this representation. 

For this boundary motion, eq. (3.15) takes the form: 


1 1 II(t) 


—+\- + 
“(G-D-D 
ay” b? to 





a ai) (3.17) 
Tr\ 1 — — 
to 


where 7» = H/J is the delay time for the Kelvin law (3.14) and 
Tr = CH((b?/a?) — 1) has the dimensions of time and is associated 
with the viscosity of the cylinder material and the elasticity of the 
casing. It can thus be thought of as a retardation time associated 
with the casing. 
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Fig. 7. Displacement vs. time plots for annihilation of the inner boundary accord- 
ing to eq. (3.16) for various ratios b?/a,?. 


Equation (3.17) can be integrated directly using the integrating 


—te/tR 
e/™ (1 - -) 
to 


which transforms the left-hand side into a complete differential, 
giving: 


factor: 








dt (3.18) 


—t/re (1 ~*~ 
. lo ‘We 
ys f £ \t + We/re) 
(1-7) 
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on 
~I 


For constant internal pressure II(t) = Ilo, this integral can be evalu- 
ated in terms of the exponential integral tabulated in ref. 7 for integral 
values of t/o,. Intermediate results can then be obtained by 
interpolation. 
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Fig. 8. Variation of the radial stress at the reinforcement with time for several 
values of the annihilation time. 


The radial stress at r = b determines the stress in the reinforcement, 
and is given by eq. (3.12) or eq. (3.9) and eq. (2.12), to be 


o, = — k(t)/b°C (3.19) 


The variation of ¢,(b,r) as a proportion of the internal pressure II» 
is shown in Figure 8, where 7 is the time ¢ expressed as a proportion 
of the annihilation time %&. The calculation was carried out for 
equality of the delay time 7» of the cylinder material and the casing 
retardation time rz, and for several values of m = to/re. The curves 
of Figure 8 thus demonstrate the time histories of the stress in the 
reinforcement for different annihiliation times for a given cylinder 
and casing. 

For m other than ©, all curves in Figure 8 pass through the origin. 
This arises since the Kelvin material, Figure 6, does not exhibit 
instantaneous elastic response, so that on first application of the 
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pressure the displacement is zero instantaneously, and no stress is 
generated in the reinforcement. For m— o, that is, long annihila- 
tion time, viscoelastic flow takes place for finite time ¢, but zero 
t/t, since 4 — ©. Thus finite stress appears in the reinforcement 
while t/t = 0. 

For m other than zero, all curves in Figure 8 pass through B, 
which indicates the full internal pressure I, on the casing at t = t, 
This must be so since no viscoelastic material is left inside the casing 
at this time. For m = 0, the solution changes discontinuously from 
C to B in Figure 8 as detailed below. 

For very long annihilation time, m — ©, delayed elastic flow 
develops fully along OA in Figure 8, and at each time thereafter, 
because of the extremely slow motion of the boundary compared 
with the delay time of the viscoelastic material. Thus the Kelvin 
material behaves essentially as an elastic material, and the curve 
AB represents a series of quasi-static elastic solutions for different 
values of the radius of the internal hole given by eq. (3.16), with no 
other time influence. 

For m = 0 annihilation occurs so rapidly that no retarded elastic 
flow has had time to occur before the material disappears, so that no 
stress occurs in the casing until the internal pressure Iy is applied to 
it at C, causing a discontinuous jump to B. 

The curves of Figure 8 show that as the annihilation time is 
reduced, corresponding to decreasing m, or correspondingly as the 
delay time is increased, the buildup of stress in the reinforcement is 
less for a given proportion of burn out. This interplay of the visco- 
elastic characteristics of the material and the stress distribution 
because of time influences in the loading procedure offers the pos- 
sibility of design choices to benefit from such influences. For example, 
in the practical case of rocket burn out, the gas pressure falls after 
reaching a maximum, and the delayed transmission of the internal 
pressure through the annihilating cylinder would provide a reduction 
in the maximum pressure for which the casing would have to be 
designed. 


4. Discussion 


The example detailed in Section 2 indicates the generality of 
approach which is possible for viscoelastic stress analysis problems 
which fall within the scope of the Laplace transform method of treat- 
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ment. Quite general viscoelastic operators can be accommodated 
to this procedure, with only the difficulty of increasing algebraic 
complexity, and full use can be made of results from the theory of 
elasticity. When this approach is not applicable, a general method of 
solution of the system of partial differential equations which rep- 
resents the viscoelastic stress analysis problem is not at present 
available, and particular approaches lead to a variety of types of 
mathematical problems. The example given in Section 3 leads to a 
nonlinear differential equation, but more general operator equations 
will arise in other situations such as integro-differential equations, 
for which at the present time special approaches must be devised. 

The examples evaluated show the influence of the material visco- 
elastic characteristics on the variation of stresses, and illustrate the 
bearing which this may have on the design of equipment using 
materials of this type. 
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Synopsis 

The formulation of the theory for determining stress distributions in linear 
viscoelastic materials is presented. The particular situations in which the time 
variable can be removed and the problem transformed into an equivalent elastic 
problem by application of the Laplace transform are stated, and the difficulties 
to be expected when this is not possible are pointed out. An example of each type 
of problem is presented, showing the general approach available in the first case, 
and the special procedure which must be devised for the latter. The resulting 
variation of the stress-distributions with time under constant surface loading 
which is associated with the viscoelastic characteristics of the material is discussed 
from the standpoint of the flexibility it offers in the design of components contain- 
ing such materials. 
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The Concentration Dependence of the Newtonian 
Viscosity of Prolate Ellipsoids 


JOHN G. BRODNYAN, Rohm & Haas Company, Philadelphia, 
Pennsylvania 


The flow behavior of suspensions of rigid uncharged spherical 
particles has been investigated quite extensively both theoretically 
and experimentally. Einstein! theoretically determined the depend- 
ence of the relative viscosity of spheres on the volume fraction ¢ 
for dilute suspensions. Mooney’ has extended this result to higher 
concentrations by a functional method. Mooney’s equation: 


2.56 
Nret = OXP ' \ (1) 


1 — ko 


reduces to Einstein’s equation at low concentrations and does success- 
fully fit the experimental data over the whole range of concentration 
with the proper choice of k. The crowding factor k is an interaction 
constant. Mooney has shown that it is most probably equal to 1.35 
or less in the case of polydispersity. An upper limit of 1.91 is expected. 

All the data’ readily available to the author on the viscosity of 
spheres has been plotted in Figure 1. Mooney’s equation with 
various values of k is drawn on this graph. It is seen that this 
equation, and therefore the functional method, has been successful 
in extending Einstein’s equation. 

The flow behavior of ellipsoids has not been so extensively studied. 
Very little data for concentrated suspensions is available. However, 
there have been a number of theories proposed for the dependen e of 
the relative viscosity of dilute suspensions of ellipsoids on the volume 
fraction and axial ratio (p). These equations, which follow, are 
dealt with in more detail in Frisch and Simha’s review.‘ 
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p* l l | 
= 1+" = 
me + le + 5 E- 2p — k) si In2p—k+1 . 


p>5 (5) 


In general the equations are of the form: 
Nret = 1 + kid + kof (p)o (6) 


Equations (2), (3), and (4) were derived by Kuhn and Kuhn* 
while eq. (5) was derived by Simha.* These equations are analogous 
to Einstein’s equation, having been derived for infinite dilution, and 
reduce to it for p = 1. 

The most interesting problem is the viscous behavior at finite 
concentrations. Therefore, an attempt has been made to extend 
the “dilute solution’’ theories to high concentrations using the method 
Mooney found useful in treating the problem of spheres in con- 
centrated solutions. 


Conditions to be Fulfilled and Assumptions 


In attacking this problem we should expect the following condi- 
tions to be fulfilled: 

(c) The final viscosity should be independent of the sequence of 
stepwise additions of partial volume fractions of the ellipsoids to 
the suspension. 

(b) The equation should reduce to one of the four dilute solution 
equations as ¢ approaches zero. 

(c) The equation should reduce to Mooney’s equation for p = 1. 
In the following derivation the following assumptions are made: 
(a’) The only interaction is that due to crowding, i.e., only first- 
order interactions are to be considered. (b’) The ellipsoids have a 
constant axial ratio p. 


Derivation 

Assume the ellipsoids are added in two volume fractions ¢, and 
@. The addition of the first fraction will increase the viscosity by a 
factor H(¢:, p) = m/no where H must reduce to eq. (6) at low values of 
¢:. The second fraction gives an increase again; H(i) = m2/m. 
In this case the first solution acts as the solvent. Due to the crowding 
effect, i.e., the interaction, the volume fraction available to ¢, is 
reduced to 1 — kg. Therefore, yx has the form h(¢:, p)/1 — kdu. 
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However, the ¢, volume fraction mutually crowds the ¢ volume 
giving H(¢:, p) = H (yx) and Yo = A(gi, p)/1 — hbo. Next we 
must consider the fact that the ¢; and ¢,. volume fractions have each 
ellipsoids interact with the rest of the members of its own volume 
fraction. This means that yo. = h(¢:, p)/1 — k(gi + 2) or Yu 
= h(¢», p)/1 — ko. In the same way rx = A(¢i, p)/1 — ke. 
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Fig. 2. (@) M. A. Lauffer, J. Am. Chem. Soc., 66, 1188 (1944). (O) M. A. 
Lauffer and W. M. Stanley, Chem. Revs., 24, 303 (1940). (x) G. Oster, P. M. 
Doty, and B. H. Zimm, J. Am. Chem. Soc., 69, 1193 (1947). (+) J.T. Yang, J. 
Am. Chem. Soc., 80, 1783 (1958). 


Since the product H(Y.2)H(¥u) is the relative viscosity of the 
suspension of total concentration (¢; + ¢:) we have: 





H(¢: + ¢:, p) = H(vw)H (vu) = a(? ag H ny 


The above equation is satisfied by: 


H(@, p) = exp ‘ 


h(¢, p) 7) 


1 — ke 


Also it is seen that ¢ must be present only to the first power. 
If h(¢, p) = kid + kef(p)d, then eq. (7) reduces to the general 
form of the dilute solution theories as ¢ approaches zero. Also, if 
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S(p) = Owhen p = 1, and k, = 2.5, Mooney’s equation is the limiting 
case for p == 1 and all our initial conditions are fulfilled. 

If we assume for simplicity k,f(p) to be of the form a(p — 1)’, 
which is the same as eq. (3), we can plot (In ,.:/¢) — 2.5 vs. (p —.1) 
on a log-log plot and obtain b from the slope and a from the intercept. 
This was done for the experimental data’»* available to the author 
and is given in Figure 2. The final empirical result is then: 


2.56 + 0.399(p — = 
1 — ko 





Nret = EXP 1 (8) 

The values for the constants are very close to those calculated by 
Kuhn and Kuhn. In fact the values compared in Table I are seen to 
agree within 2%. 








TABLE I 
Experimental Calculated 
Constants values values 
a 0.399 0.4075 
b 1.48 1.508 





Comparison with Experimental Data 


The data on poly-y-benzyl-1-glutamate solutions given by Yang® 
are probably the best to use for comparison with eq. (8). This is 
true for three reasons: (1) his material has been proven essentially 
monodisperse; (2) his data are shown experimentally to give vis- 
cosities in the Newtonian region and this has not been shown for 
the other data; and (3) his data have also been taken in a higher 
range of volume fractions where we would notice deviations from 
dilute solution behavior. 

From his reported molecular weight and the known density and 
width of the particles an axial ratio of 93.3 was computed. In 
Figure 3 the lines for three equations with p equal to 93.3 have been 
plotted as well as Yang’s data. The equations plotted are eq. (4), 
the result if eq. (4) had been treated in the manner given above in the 
derivation, and eq. (8). 

It can be seen that eq. (8) fits the data better than either of the 
other two equations and is certainly much better than the dilute 
solution eq. (4). However, the fit is far from exact. Nothing has 
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been said about the value of k used, because the volume fractions 
are so small that it would not have mattered if we neglected con- 
sideration of the excluded volume effect. 


Evaluation of / 


As we go from spheres of axial ratio one to ellipsoids of axial ratio 
greater than one,.it becomes obvious that the cubic close packing 
assumed by Mooney to obtain k = 1.35 is no longer possible. How- 
ever, as the concentration of a suspension of ellipsoids increases, 
we would expect an ordering to take place, with the ellipsoids lining 
up with their Jong axes parallel. In fact, Vold® showed theoretically 
that this is the case if van der Waals forces are the chief forces of 
interaction. We could not expect as dense packing as would occur 
with spheres but a possible picture can be visualized in Figure 4. 
In this case p = 1/b. The volume of a unit cell would be b’/, and 
we see that one ellipsoid should be contained in this cell (four quarter 
ellipsoids). The volume of this ellipsoid would be (2/6)b*l. There- 
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fore, the volume fraction taken up by the ellipsoids when we approach 
densest packing would be: 2/6 = dmaz- 

Since k = 1/¢mez we have k = 1.91 for ellipsoids with a large 
value of p. However, we have seen that for p = 1 we have a value 
of k equal to 1.35. We, therefore, expect k to approach 1.91 rapidly 
and asymptotically as p increases. 


Discussion 
The choice of k2f(p) and the form of the equation seem to be correct 
as a first approximation. However, the experimental viscosity 
does not seem to increase quite as fast as eq. (8) would predict. It 
is not quite clear what alterations have to be made in the theory, 
and any modification will have to await further experimental data. 
Similarly the choice of k will also have to rest on experimental proof 


unavailable at present, but 1.91 probably will suffice as a first approx- 
imation. 





Figure 4. 


Since Yang* has shown that the viscosity of rigid ellipsoids as a 
function of rate of shear is given by Scheraga’s"” numerical tabulation 
of the Jeffery-Peterlin-Saito treatment, the viscosity of a con- 
centrated solution of ellipsoids may now be determined to a first 
approximation knowing only the axial ratio p and the volume frac- 
tion ¢. The rotary diffusion constant can be estimated from Per- 
rir’s' equations if the length of the major axis is known. 
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Synopsis 

Kuhn and Kuhn’s theory for the viscosity of ellipsoids in dilute solutions has 
been extended to concentrated solutions using Mooney’s functional method. The 
resulting equation is: 


m exp {250+ 0.2001 — 1)844 
Nrel = OX] } i my 


The constants 0.399 and 1.48 were obtained from the available experimental data 
and are found to be within 2% of the theoretical values. The dependence of the 
relative viscosity on volume fraction (@) is seen to agree with the experimental 
results presented by Yang. 
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An Analysis of an Absolute Torsional Pendulum 


Viscometer 


ELLIOT A. KEARSLEY, National Bureau of Standards, 
Washington, D.C. 


Introduction 


At present, the only accurate absolute method of measuring the 
viscosity of Newtonian liquids is based on Poiseuille flow through 
capillaries. In fact, accuracy better than 0.1% requires such élab- 
orate care that only the viscosity of water has been directly measured 
with this accuracy. Viscosities of other fluids are customarily 
measured with relative viscometers calibrated either directly or 
indirectly against water. 

It is of interest to develop an independent absolute method of 
measuring viscosity of Newtonian fluids comparable in accuracy to 
the technique of Poiseuille flow, i.e., a method such that all errors 
can be bounded in some practical sense not involving comparison 
with fluids of known viscosity. 

A viscometer originally considered by Helmholtz' consists of a 
spherical shell suspended from a torsion wire and filled with fluid 
to be tested. The sphere is excited in torsional oscillation and the 
subsequent motion observed. The liquid filling the shell will affect 
the oscillation in two ways: (1) damping will be enhanced due to 
consumption of mechanical energy by viscous heating of the liquid; 
(2) the period of the torsional pendulum will be lengthened by the 
inertial effects of the liquid carried along with the shell. 

In principle, a measurement of either the period or the damping of 
the torsional oscillation will determine the viscosity of the fluid. 
Practically, the values of the various parameters of the pendulum- 
torque constant, radius of shell, mass of shell, density of fluid-—will 
determine which measurement specifies the viscosity with greater 
accuracy. Although torsional pendulum viscometers have been used 
69 
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in the past,’ they have always been designed with parameters such 
that the damping rather than the period was a sensitive function of 
the viscosity of the test fluid. In view of the highly developed 
methods of measuring time intervals now available, it appears that 
the measurement of period might be the more precise method. This 
paper reports on an investigation now in progress whose goal is an 
answer to the following question: What should be the values of the 
design parameters of an oscillating sphere viscometer such that a 
measurement of period will give the most precise absolute value of 
viscosity for test fluids of a certain nominal viscosity range? 


Calculation of the Viscous Torque 


The calculation of the period as a function of the viscosity of the 
fluid is done in two steps: (1) the fluid flow problem is solved and the 
viscous torque on the shell is calculated; (2) the calculated viscous 
torque is inserted into the dynamic equation of the shell and the mo- 
tion of the shell is then found. 

Figure | is a schematic diagram of the pendulum with the spherical 
coordinates used in the flow problem indicated. The flow problem is 
solved in the approximation of small angle of rotation of the shell so 
that quadratic inertial terms can be ignored. The physical meaning 
of this approximation is that the effects of fluid flow in planes con- 
taining the axis of rotation of the pendulum are ignored. If the 
pendulum is to serve as an absolute viscometer, this error must 
subsequently be estimated. 

The Navier-Stokes equation in spherical coordinates, ignoring 
quadratic inertial terms and taking account of symmetry, is 


ee EA ee ee, a 
ot =r ft el: e+ ss aed r? sin? 6 (1) 


where w is the fluid velocity perpendicular to planes containing the 
axis of rotation (the ¢ component of velocity), and v is the kinematic 
viscosity. This is a linear homogeneous equation in w and the time 
dependence can be separated out by assuming* 


w = ve™ (2) 








* As is usual, when real physical quantities are set equal to complex expres- 
sions, the real part of the complex expression is implied. Notice that 2 is in gen- 
eral a complex quantity. 
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The boundary conditions applicable are then 
v = 120R sin @, whenr = R (3) 


v is continuous and finite within the sphere, 0 is the angle through 
which the sphere is turned from equilibrium at ¢ = 0. 
The solution of this boundary value problem is 


iN0R sin zx — zx cos 2x 


v= 





Fe sin z — zcosz (4) 
z2=iVi0/yR, x«x=r/R 
The complex parameter z is closely related to boundary layer thick- 
ness. 
Finally the viscous torque due’ to this flow pattern acting on the 
shell is 


; ; 2? sin z : 
L= (8/3) yiaoe™ ( $$ — 3) (5) 
sin z — zcosz 


LLL 








Fig. 1. The spherical coordinates used in the flow problem. 
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where x is the coefficient of viscosity; that is, density times kinematic 
viscosity. 


The Dynamic Equation of the Shell 


The dynamic equation of the shell is obtained by equating the 
torques acting on the shell to the rate of change of angular momentum 
of the shell, viz.: 


10*@e™ = Kee™ + L (6) 


where J is the moment of inertia of the shell and K is the torque 
constant of the suspension wire. Substitution of eq. (5) into eq. (6) 
and some manipulation results in the following equation, 


—10? + i2(8/15) wR%yz* Gan _ “) +K=0 (7) 
1 — zctgz 2? 

Up to this point the analysis has not been anything new. The remain- 
ing problem, to solve this equation for Q, is in a sense only arithmetic. 
Because © is contained in z, this is a complex transcendental equation 
in 2. Previous authors have developed various approximate solu- 
tions of this equation valid for certain ranges of the parameters, 
but in general not valid for values of the parameters such that a 
measurement of period will accurately determine the viscosity. 


Computations 


Before setting up the equation for solution on a computer, it is 
advisable to put it in nondimensional form so as to reduce the number 
of computational parameters. Accordingly, 2 is broken up into 
real and imaginary parts and the following quantities are defined: 


Q = (24/T) + 2 


In = (8/15) 2R'p moment of inertia of fluid in rigid body 
rotation 

a=I/I, ratio of moments of inertia of shell to 
fluid 

w = VK/I circular frequency of empty shell (8) 

t = (V/a w/2r)T dimensionless period 

n = v/R*/ a wo dimensionless viscosity 


D = 8T/x decrement per period 
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Using these definitions the computational problem reduces to 
finding a one-parameter family of curves. In particular, after de- 
fining a new function 


5 15 
&(z) = —— =i (9) 
1 — zctgz 2 
eq. (7) can be written in extremely compact form, 
@(z) = a + 1/n?2* (10) 


Mathematically, it is necessary to study the function #(z) and then 
to look for roots of eq. (10). For given values of a and 7, eq. (10) 
defines certain complex values of z. Then, 


ten (s'0(6)) 








Im(z‘*) 
1/n? = Re(z*b(z)) — aRe(z*) (11) 
1/r = —ylm(z’) 

Re(z*) 

Im(z?) 


in which Re and Im mean “the real part of” and “the imaginary part 
of,” respectively. 

For a given set of parameters, a and 7, there is one primitive com- 
plex root of the equation and an infinity of real roots. The former 
corresponds to a damped oscillation and the latter to purely damped 
components of the motion. The complex root determines the period 
of the oscillation, and the real roots are of interest only in that they 
may affect the experimental measurement of period by causing a 
“zero drift” error. 

In practice, the IBM 704 computer was used in the following way. 
First various complex values of z were inserted in eqs. (11) and the 
associated parameters were computed. After sufficient data were 
developed in this way, it was possible to make a good guess at the roots 
of eq. (10). The procedure was then reversed; values of a and 9 
and an estimated root, z, were fed to the computer. The root z was 
then corrected by Newton’s method and then r, D, ete. were calcu- 
lated. A typical line of output from the computer contained col- 
umns labelled a, 7, Re(z), Im(z), 7, D. 
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Figure 2 graphs some results of the computations, 7 as a function of 
» for various values of a. As might be expected, r ranges from +/a 
for vanishing 7 to ~/a + 1 for infinite n, corresponding respectively 
to the periods of the empty shell and the shell and fluid rotating as a 
rigid body. Figure 3 is included for completeness, and shows the 
decrement per period as a function of 7 for various values of a. 
The fact that the decrement for a equals zero approaches a finite 
limit for vanishing viscosity is a mathematical phenomenon due to the 
fact that D is not an analytic function of either a or 7 about the 
point a = 7 = 0. For a = 0 and vanishing viscosity the period 
becomes infinite and the energy loss per unit time vanishes in such a 
way that the product is finite. The slightest inertia of the shell will 
have a dominant effect for small enough viscosity, so that for any 
finite a the energy loss per period will vanish for vanishing viscosity. 
Figure 4 shows the damping against the logarithm of 7 and illustrates 
this fact more clearly. 

Figure 5 shows the dimensionless period as a function of log 7 
for various values of a. Notice that, independent of a, there is a 
flattening of the curve for 7 between 10~' and 10-*. 
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Fig. 2. Dimensionless period as a function of dimensionless viscosity for three 
values of ratios of moments of inertia, a. 
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Fig. 3. Logarithmic decrement in nepers as a function of dimensionless viscosity 
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Fig. 5. Dimensionless period as a function of logarithm of dimensionless viscosity 
for three values of a 


Estimation of Errors 
The most crucial error seems to be in the determination of the radius 
of the sphere. From eq. (8) we have, 


y= R°V a won =R°VK/Iyn (12) 


oR/r \R~ 2h * Ro)” (: 
In this notation a subscript on a partial derivative indicates the quan- 


tity held constant as the derivative is taken. In taking the partial 
derivative of the moment of inertia with respect to radius, it is nec- 


Hence 
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essary to specify how the moment of inertia actually will be com- 
puted. Since the inside radius of a sphere is difficult to measure, it 
is best to find the mass of the fluid by weighing and then calculate the 
moment of inertia from 





Ip = (2/5)MR? (14) 
Then, 
To calculate 0n/OR it is easiest to transpose eq. (8) to give 
2 —s 
T = ~< ‘= eV 1,/K T (16) 


Va w 


Taking a derivative of this equation with respect to R and holding T 
constant, one gets, after some substitutions, 





_1 ,1/0r% 2a 2) 
0-5 +2 (Fo R da (17) 
and hence, 
On = 2a(0r/0a) — + g 
OR R@r/dn) (18) 
Finally, 





R [ov 2a(d0r/0a) — ") 

= = {1 

v () ( + 07/0 log 7 (19) 
This is an equation for the ratio of relative errors in vy and R using 


a measurement of period to determine viscosity. A similar equation 
for the ratio of relative errors in » and T' is, 


7 (%) ne, Se - 
v \OT/e dr/dlogn (20) 


The quantities on the right of eqs. (19) and (20) can all be calcu- 
lated graphically or numerically from Figures 2 and 4. The ratio of 
relative errors in vy and R is sensibly 2 for y less than 10~'; for higher 
values it rises rapidly. The ratio of relative errors in vy and T has a 
minimum as a function of y for fixed a. AT min(a) determines a curve 
in the 7, plane which is shown in Figure 6, superimposed on the curves 
of Figure 5. Also shown are contours for (7'/v)(0v/OT)z equal to 
10 and 30. 
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A further error calculation has yet to be completed; the uncertainty 
in viscosity due to the effects of secondary flow. As was pointed out 
earlier, the inclusion of quadratic inertial terms in the Navier-Stokes 
equation leads to fluid flow in planes through the axis of rotation of the 
pendulum (see Fig. 7). This in turn has a second-order effect on the 
torque by distorting the boundary layer of fluid and therefore intro- 
duces a correction term to eq. (5). 

















0.0 4 4. 7’ 1 rt 
-§ a4 -3 -2 -| 0 
Log % 


Fig. 6. The curve determined by minimum ratio of relative errors in period and 
viscosity and two curves of constant ratio of relative errors in period and viscosity 
superimposed on Fig. 5. 


Carrier and DiPrima* have calculated a correction term for the 
torque on a sphere in torsional oscillation immersed in fluid. The 
same technique can be applied to the interior flow problem of interest 
here. This calculation is now in progress, and although it is clear 
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that an exact solution to the problem is virtually impossible, the 
analysis can be used to give an upper bound on the correction term. 

By analogy to the Carrier and DiPrima problem, it can be argued 
that the correction term to the torque will be proportional to the 
square of the amplitude of the oscillation and will be a monotonically 
increasing function of the radius of the sphere. These considerations 
alone are enough to show the qualitative effect of the secondary flow 
ecrrection. 


a 





Fig. 7. Secondary flow in the spherical torsional pendulum viscometer. 


It was mentioned earlier that eq. (10) defines a set of purely damped 
components of motion as well as a damped oscillation. Experi- 
mentally, the former will appear as a zero-point drift in measuring 
the period of the oscillation. In practice, it will be necessary to allow 
these components to become ignorably small before beginning to 
measure the period. Because the oscillating component of motion is 
itself damped, the initial amplitude of motion must be large enough 
to allow a practical determination of period after this waiting period. 
The initial amplitude of the oscillation must, therefore, be a compro- 
mise—large enough for an accurate period measurement, small enough 
to keep the secondary flow correction within bounds. In general, 
the larger the ratio of moments of inertia, a, the smaller the damping 
of the oscillation and the easier this compromise is to make. 
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Determination of Best Parameters 


The question of the best set of parameters for the torsional pendu- 
lum viscometer is now beginning to unravel. The following proce- 
dure can be adopted to answer the question. 

(1) From the studies of relative error, a best point (point of mini- 
mum error) can be found on the r,y plane of Figure 6. 

(2) Equations (8) are then sufficient when combined with an esti- 
mate of viscosity and density of the test fluid to select a set of physical 
parameters such that r,y, and a correspond to the best point deter- 
mined by step (1). 

(3) Various practical requirements must then be reviewed and may 
force some modification of the choice of operating point. Typical 
of the practical requirements are reasonable radius of sphere, and rea- 
sonable strength of suspension. 

Finally, although the problems of thermostating, air resistance, 
damping in suspension, etc., have not been discussed here, it is clear 
that they must be considered in a practical design of a torsional pendu- 
lum viscometer. They do not appear to be insurmountable’ prob- 
lems. 


Dr. Andreas Schopf of The American University and Mr. Stanley Prusch of 
the Computation Section of the National Bureau of Standards contributed the 
formulation and coding of the problem for the IBM 704. 
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Synopsis 


Both the period and decrement of a hollow sphere, filled with fluid and used as 
a torsional pendulum, depend on the viscosity of the fluid. It appears that a 
measurement of period might give an absolute measure of viscosity comparable 
in accuracy to measurements by Poiseuille flew. The relative errors in viscosity 
due to relative errors in measurement of period and radius of sphere are calcu- 
lated. A technique is described for selecting the best set of physical parameters 
for such a viscometer. Final results await the solution of a problem in secondary 
flow. 
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Viscoelastic Behavior of Polyethylene and 
Polypropylene 


JOSEPH A. FAUCHER, Research Department, Union Carbide 
Chemicals Company, South Charleston, West Virginia 


There have been frequent statements in the literature to the effect 
that the method of reduced variables, which has been so successful 
with amorphous materials, is inadequate to describe the mechanical 
behavior of crystalline polymers. These objections rest largely on 
two grounds: first, that the structure of crystalline polymers changes 
with temperature (so that the convenient time-temperature equiv- 
alence is invalid), and second, that these materials are nonlinear in 
their viscoelastic behavior. Although these objections are well- 
founded, the result has been that the mechanical behavior of crystal- 
line polymers has been left in a purely empirical state. In order to 
escape from this undesirable situation, we have tried to retain the 
essential features of reduced variables by determining under what 
range of conditions the theory might be valid for these materials. 
This paper is a description of these conditions and a summary of 
results obtained for the particular cases of polyethylene and poly- 
propylene. 

With respect to the first objection above it is known that many 
crystalline polymers maintain the same structure over a fairly wide 
temperature range: i.e., their per cent crystallinity vs. temperature 
curves have a long plateau and break over rather sharply near the 
melting point; consequently, if one stays in the plateau region the 
first objection will not apply. Two samples of polyethylene were 
used for the present study: DYNH, a low-branched material, and 
Marlex-50, which is essentially linear. Their crystallinity-tempera- 
ture curves are given in Figure 1 (data derived from density measure- 
ments). The plateau for Marlex-50 extends as far as 70°C., while 
for DYNH it is somewhat less, reaching to about 20 or 30°C. Similar 
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Fig. 1. Crystallinity-temperature curves for polyethylene. 
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Fig. 2. S/R vs. time. Marlex-50 at 40°C. 


measurements were not made on polypropylene, but it is clear from 
the Arrhenius plot of Figure 12 that the plateau extends to nearly 
120° in this case. 

The second objection—that of nonlinearity—can be met by a 
restriction of another kind, namely in the amount of strain. Under 








VISCOELASTIC BEHAVIOR OF POLYETHYLENE 83 


sufficiently low strain any material behaves linearly. In this con- 
nection Alfrey' has pointed out a convenient means of testing for 
linear behavior. It is based on the following formula for stress relaxa- 
tion which can be derived from the “infinite Maxwell Model”: 


S(t) = R{*2 cH(2)(1 — e-*)dinz (1) 


where S is the time-dependent stress, 2 is the rate of strain, and H(x) 
is the distribution of relaxation times. It is clear that S/F is a func- 
tion of time only. Consequently, from a series of stress-time curves 
at various strain rates one should be able to plot S/R vs. time and ob- 
tain a single master curve. No data of this kind for crystalline 
polymers have appeared in the literature, so Marlex-50 was tested at 
a number of strain rates varying from 0.1 to 25%/min. The plot of 
S/R vs. time in Figure 2 shows that a single curve is not even ap- 
proximately realized. The curve at 0.1%/min. strain rate has the 
least total strain and hence best represents the “‘ideal’’ master curve, 
if one indeed exists. It is clear that the larger the absolute strain, 
the greater is the deviation from the ideal curve. Even the values 
at 1% strain (e.g., 60 sec. on the 1% curve) are appreciably off. 
Accordingly, relaxation data for Marlex-50 were taken at a maximum 
of 0.5% strain. From the figure it is apparent that at this low strain 
deviations from the master curve will be quite small. For DYNH 
a similar series of curves shows that significant deviations do not 
occur till 1% strain is reached; this material is only about half as 
crystalline as Marlex-50. The stress relaxation data shown below 
for DYNH were taken at 1% strain. Polypropylene was not tested 
in the manner just described, but was assumed to be similar to Mar- 
lex in this respect. The relaxation data for it were taken at 0.5% 
strain. 


Experimental 


The polyethylene samples used were Marlex-50 and DYNH; both 
are commercially produced and have been frequently described in the 
literature. Crystalline polypropylene was prepared in the laboratory 
using a Ziegler-type catalyst; its density was0.912. Assuming values 
of 0.936 for crystalline polypropylene and 0.850 for amorphous poly- 
propylene, this material has a weight per cent crystallinity of about 
75%. The “amorphous” sample used for comparison was made by 
taking a low-crystalline polypropylene and extracting it successively 
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Fig. 3. Stress relaxation of Marlex-50. 


with acetone (to remove oils) and ethyl ether. The fraction dissolv- 
ing in ethyl ether was precipitated by methanol, vacuum dried, and 
made into plaques. Its density was 0.856, indicating a crystallinity 
of about 8%. 

The mechanical properties of these materials were investigated by 
the method of stress relaxation. A sample of known dimensions was 
extended at about 5%/min. to a total strain of 0.5% (1% for DYNH) 
on the Instron tester and held while the stress was measured over a 
time of 10 to 2000 sec. This was carried out at 10° intervals from 
—60° to above 100°C. The tensile modulus Z(t) can then be calcu- 
lated by: 


E(t) = S() V/A Al (2) 


where S(t) is the measured stress in dynes, A is the cross-sectional 
area in sq. cm., / is the sample length and Al is the extension. The 
data for the four materials used are shown in Figures 3-6. 

The time-temperature equivalence which has been applied to 
amorphous polymers?” relates E(t) at two different temperatures 
by: 


(To/T)Er(art) = E,,(t) (3) 


where 7; is an arbitrary standard temperature and (art) is the “‘re- 
duced” time. The factor T,/T is introduced to correct for the effect 
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Fig. 6. Stress relaxation of amorphous polypropylene. 


of rubber elasticity which increases with absolute temperature. 
Since our materials are about 75-80% crystalline, it is not at all clear 
that rubber elasticity makes any contribution in this case. Conse- 
quently, we have decided to omit this factor (which is not a large 
correction in the temperature range here). The temperature reduc- 
tion used then is 


Er(art) = Ez,(t) (4) 


(For the case of amorphous polypropylene the 7,/T factor has, of 
course, been retained, and the data of Fig. 6 are corrected by this 
amount as is indicated on the ordinate.) 

When the temperature reduction is carried out, one converts the 
series of lines in Figures 3-6 to a single master curve in each case. 
These relaxation curves, using 20°C. as the standard temperature, are 
shown together in Figure 7. It is apparent that the crystalline 
materials have very similar relaxation behavior; indeed, polypro- 
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Fig. 7. Relaxetion curves of polyethylene and polypropylene at 20°C. 
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Fig. 8. Distribution of relaxation times for polyethylene and polypropylene. 


pylene and Marlex-50 are nearly indistinguishable over a long range 
of times. On the other hand, amorphous polypropylene has a re- 
laxation curve much like polyisobutylene, which has been rather 
thoroughly investigated in recent years.‘ 

From the relaxation curve we can calculate the distribution of 
relaxation times, H(t), which is defined by the integral equation: 


E(t) = f72 H(z)e-""dln x (5) 
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The symbols used throughout this paper are those proposed by the 
Committee on Nomenclature of the Society of Rheology.’ H(t) 
can be derived from experimental data by the relation: 


H,(t) = —E(t) (d log E/d log t) (6) 


where the subscript 1 denotes a first approximation. The curves of 
log H(t) vs. log t are shown in Figure 8; as expected, the crystalline 
polymers are all quite similar. Apparently such flat curves extending 
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Fig. 11. Shift factors for Marlex-50. 


over many decades of log time are characteristic of crystalline ma- 
terials; at present there is no theory to predict this behavior. On 
the other hand, amorphous polypropylene shows the well-known 
slope of —'/: in the transition region (log t = —6 to 0) required by 
the Rouse theory,® while the tailing off at higher time seems to indi- 
cate considerable polydispersity in the sample.’ 

Since relaxation spectra of crystalline polymers are still rather 
novel, it seemed desirable to test the distribution function just 
obtained in some manner. This was done in the case of Marlex-50 
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using the spectrum to calculate a stress-strain curve according to the 
equation: 


S(e.) = Rf22 cH(x) (1 — e-*/*)dilnz (7) 


where « is the fractional strain and R is the strain rate. To perform 
the integration an average value of H was taken in each decade of log 
time. The results with H,(t) as given by eq. (6) were compared 
with an experimental stress-strain curve at 1%/min. strain rate. The 
agreement was not very impressive, so a second approximation to 


H(t) was made according to: 
d log E d log E\? E @ilogE 
nao = —2 EE , yp (Ak), oe ee 
d log t d log t 2.303 d log t? 
(ref. 8, eq. 20) where y is Euler’s constant, y = 0.5772. The second 
approximation gives a better calculated stress-strain curve, as can be 





(8) 
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Fig. 13. Shift factors for amorphous polypropylene. 


seen in Figure 9. A third-order approximation would probably give a 
still better result, but the data hardly warrant such tedious calcula- 
tions. The point is, that since we have used data from —70 to 
+70°C. to get the calculated curve, it seems that the time-tempera- 
ture equivalence concept does indeed apply in this case. 

As is well known, the magnitude of the a, shift factors can be used 
to determine the activation energy of the stress relaxation process. 
In Figures 10-12 are shown the plots of log az vs. 1/T for the crystal- 
line polymers DYNH, Marlex-50, and polypropylene. In each case 
a reasonable straight line is obtained. The slopes of these lines 
(times 2.303R) give the heats of activation—about 30 kcal. for poly- 
ethylene and 50 keal. for polypropylene. For amorphous polymers 
this activation energy is for a process of viscous flow; at present it 
is not clear just what process is involved with the crystalline materials. 
An elucidation of this point would go far to explain the rather large 
difference in activation energy of polypropylene and polyethylene, as 
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well as the rather interesting fact that DYNH and Marlex-50 
have the same activation energy, although their crystallinities differ 
by a factor of 2. 

It should be noted for the polyethylenes that deviations from the 
straight-line plot begin to occur at those temperatures where crystal- 
linity starts to decrease, i.e., about 30°C. for DYNH and 80°C. 
for Marlex-50. For polypropylene, deviations begin at about 130°, 
indicating that the crystallinity stays constant up to this point, 
although this was not actually verified by density measurements. 

The log ay vs. 1/T curve for amorphous polypropylene is shown 
in Figure 13, and the heat of activation derived therefrom is given 
in Figure 14. The maximum observed at —10 to —20° probably 
represents the glass transition for this material. 

In conclusion its seems appropriate to mention the work of Ké 
on polycrystalline aluminum between 150 and 225°C.* in which the 
method of stress relaxation was used with very small strains. The mas- 
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ter relaxation curve obtained was very similar to our Figure 7 
(Marlex-50 and crystalline polypropylene) and the plot of log a, vs. 
1/T yielded a straight line whose slope indicated an activation energy 
of about 34 keal. It seems likely that the relaxation curves of Figure 
7 (excepting amorphous polypropylene) are typical of crystalline 
materials generally and are quite insensitive to details of chemical 
structure, unless these latter may influence the position of the curve 
on the log E and log ¢ axes. 
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Synopsis 


An attempt is made to apply the theory of reduced variables to the crystalline 
polymers polyethylene and polypropylene. It appears that this can be done if 
suitable restrictions are made in temperature range and amount of strain. Data 
were obtained by the method of stress relaxation. Curves at different tempera- 
tures can be superposed to give a single master curve, and the shift factors so ob- 
tained fall on a straight line in an Arrhenius plot. The distributions of relaxation 
times are calculated for both polymers and are found to be very similar. The two 


are compared with amorphous polypropylene, whose behavior is much like poly- 
isobutylene. 
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Thermodynamics of Transport Processes 


in Liquids 


R. E. NETTLETON,* Department of Chemistry, The Rice Institute, 
Houston, Texas 


To discuss momentum and energy transport in an infinite, isotropic, 
nen-uniform fluid, it is convenient to divide the fluid into sub- 
macroscopic volume elements, each containing several million 
molecules. The non-equilibrium state of one such volume element, 
hereafter called the ‘‘system,” is described by the local density and 
temperature and by additional variables which fall into two classes. 
Type I variables describe the degrees of advancement of chemical 
reactions,' the occupation numbers of internal rotational and vibra- 
tional levels, degree of molecular association, and local concentration 
of free volume. Type II variables are the space-gradients of the 
temperature 7’ and Eulerian velocity of flow, u, and first and higher 
time derivatives of these variables. A non-equilibrium thermo- 
dynamic treatment using Type I variables has recently been given by 
Davies,? while Prigogine (see ref. 3, chapter 7) has developed the 
thermodynamics of continuous systems, using Type II. We shall 
show that these two descriptions are related through the equations 
giving the rates of change of Type I variables in such a way that in 
the limit of low frequency, the first of these two descriptions reduces 
to the second. 

The demonstration of this result employs the Onsager-Casimir 
reciprocity theorem.’ This theorem asserts that when the irreversible 
part of the rate of entropy production is given by a sum of products of 
fluxes J; and conjugate forces X; and no magnetic field is present, 
then in the linear force-flux equations J; = >) Ly X,, we have Ly = 

> | 
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L» whenever the corresponding forces X,; and X, are both odd or 
even with respect to time reversal, while Ly, = — L, when one force 
is odd and the other even. Proof of this theorem*~ from the prin- 
ciple of microscopic reversibility assumes the J; are time derivatives 
of Onsager variables which vanish in equilibrium and that the 
X;, are linear sums of such variables. We shall proceed to interpret 
thermodynamically the relaxation-rate equations which can be 
derived from physical models of the system on the basis of the 
fact that they contain terms having a mathematical form analogous 
to that usually assumed for the J; and X,. 

We consider first a model for shear relaxation. Frenkel* proposed 
that a liquid should behave like an elastic solid when sheared suf- 
ficiently rapidly so that the molecules have no time to rearrange them- 
selves through viscous flow. Thus we can envisage the system 
divided into many microscopic subregions, in some of which the 
molecular planes are sheared and in some of which the shear has 
relaxed, and introduce a strain tensor ¢‘; to represent the average 
elastic strain over all subregions. «‘; relaxes to a diagonal tensor 
with the Maxwell relaxation time 7, as the system approaches me- 
chanical equilibrium, and thus o‘, = «'; — '/3*, 5‘; relaxes to zero. 
This fact can be represented by the equation: 


a*; —=— — 1 a; + B‘, (1) 

v1 
The tensor B‘; = /{Vu + vu'}‘; — '/V-uéd‘,; is the diver- 
genceless rate-of-deformation. If all variables vary with time as 
e* a‘, and B‘, are both proportional to w, so that as w — @, the 
itiddle term in eq. (1) is negligible and ¢‘; = B';, in accord with 
Frenkel’s hypothesis, which is thus seen to demand the presence of 
the term B‘;, To examine the significance of this term, we must 
make a thermodynamic interpretation of all the terms in eq. (1). 
Since o‘, is an Onsager variable, we interpret the components of a‘; 
as thermodynamic fluxes, by analogy with the general theory dis- 
cussed above. We suppose that the term linear in o‘; is proportional 
to a thermodynamic force 2‘; = —o‘;/w7; conjugate to o*;, m 

being a phenomenological coefficient. 
The thermodynamic interpretation just made allows us to apply 
the Onsager theorem, provided we make the additional assumption 
that the rate of entropy production contains the term 0"; ¢’;, a sum of 
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products of the fluxes and forces identified in eq. (1). From the 
Onsager theorem, since the flux ¢‘; depends linearly on — B‘, which is 
thermodynamically conjugate to the pressure P‘, divided by the 
number density n, and which is odd with respect to time reversal, 
then P’, must depend linearly on 2", according to the equation: 


P', = P3; 4+ no; (2) 


where P is thermodynamic pressure given by the equilibrium equation- 
of-state. 

Since 2‘, is proportional to the divergenceless elastic strain tensor, 
eq. (2) indicates that the shear modulus of the liquid is: \, = n/y7:. 
As w > 0, a‘; > 0, and Q*, — —B*,/y,. Substitution of this result 
into eq. (2) shows that the low frequency shear viscosity is: 


n = ‘/dorn (3) 


Thus we see how P*, in general depends on the Type I variable 
0‘, but at low frequencies depends on the Type II variable B‘,, 
the transport coefficient 7 being determined from our thermodynamic 
interpretation of eq. (1). 

One can calculate 7; with the aid of a recent theory by Mooney’ 
and ), from the expression pc? = B.. + */3\, for the velocity c of 
longitudinal hypersonic acoustic waves, were p is the mass density. 
B. is the high-frequency bulk modulus, discussed below. The 
velocity c in turn can be calculated from the thermal conductivity 
with the aid of a theory of heat conduction due to Debye,* which 
assumes that it is the high-frequency sound waves which carry thermal 
energy throughout the fluid. Calculated values of 7; and 2n/): 
have been found to agree as to order-of-magnitude in carbon tetra- 
chloride, acetone, and chloroform.’ 

To calculate B.., we must discuss the relaxation of compressional 
stress, which we do in a manner mathematically analogous to the 
discussion of shear relaxation. In non-associated liquids, sound 
absorption is attributable to a lag in the energy transfer between 
translational and vibrational degrees of freedom.” If ¢ is the 
fraction of molecules in the system in the lowest excited vibrational 
level, ¢ will relax to an equilibrium function {°(p,7) of the density, p, 
and temperature with a relaxation time 7; which can be calculated 
from a theory of Schwarz and Herzfeld. If §& = ¢ — %, the 
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Onsager variable ¢ will relax to zero, which can be expressed mathemat- 
ically by the equation: 


g= —+¢4rvu (4) 

Te 
The presence in eq. (4) of the term linear in V-u = —)/p= —4, 8 
being the condensation, is explained by considering the limiting form 
of eq. (4) for frequencies w >> 7,~'. At such frequencies, ¢ has no 
time to relax and thus assumes a fixed value, {. Under these condi- 


tions if wave propagation is adiabatic, then § = f) — ¢° = —(0f*/ 
ds),8, « being the entropy per molecule. Since £ and ¢ are both 
proportional to w, we find that when w >> 7.~', § = —TI's, so that 


Ir = (0¢°/ds),. The presence of the term in V-u also suggests a 
thermodynamic interpretation of the terms in eq. (4), analogous to 
that made for the terms in eq. (1). We can suppose that & is a flux 
and that 6 = —£/yr, is the conjugate force, u(p,7') being a phe- 
nomenological coefficient. Ifthe rate of entropy production is assumed 
to contain the term #£, we can apply the Onsager theorem to eq. (4) 
and conclude that because the rate equation for £ is linear in —V-u, 
the diagonal elements in P‘,/n must be linear in 4, i.e., that: 


P= P+ nl@ = P — nTt/un, P=); PF. (5) 


If we put E = iwt in eq. (4), solve for &, and substitute in eq. (5), 
we obtain the complex bulk modulus B(w) : 


iwnT?/p 


Biw) = B, + 6) 


where B, is the equilibrium adiabatic bulk modulus. When w >> 

v2, eq. (6) reduces to B. = B, + nI?/ure, As w—> 0, we obtain 

the bulk viscosity 4. = lim Im[B(w)/w]) = nI?/y. These results 
0 


have been used to calculate both B.*° and '* for non-associated 
liquids, the calculated results being in qualitative agreement with 
experiment despite the fact that most liquids exhibit a distribution of 
relaxation times. Such a distribution gives rise to a whole set of 
normal coordinates obeying rate equations mathematically similar 
to eq. (4), while B(w) contains a linear sum of complex terms like 
the term on the right in eq. (6), there being one term in the sum for 
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each relaxation time. A thermodynamic treatment of systems 
having many relaxation times has been developed.” 


The author gratefully acknowledges support from a research grant made by 
the Robert A. Welch Foundation. 
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Synopsis 


Two types of variables are distinguished in the local description of a non- 
uniform fluid. Type I variables determine the molecular structure of the fluid 
and degrees of advancement of chemical reactions, while Type II variables are 
the gradients of temperature and velocity which appear in the hydrodynamic 
description of a continuous fluid. The rate equations for Type I variables also 
depend on Type II variables, so that at low acoustic frequencies one can solve 
for the former in terms of the latter. When specific terms in the rate equations 
are identified as forces or fluxes in the sense of irreversible thermodynamics, one 
can use Onsager’s theorem to calculate the dependence of the pressure on Type 
I variables, which in turn enables one at low frequencies to determine the de- 
pendence of pressure on the fluid velocity gradient and to calculate the shear and 
bulk viscosities. 
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Secondary Stresses in Anisotropic 


Amorphous Media 


CHARLES MACK, Research Department, Imperial Oil Limited, 
Sarnia, Ontario, Canada 


I. Introduction 


The presence of secondary stresses was first observed by Poynting,' 
who found that steel and copper wires elongated and increased in 
volume when subjected to a torsional torque. This result is, not in 
accord with the predictions of the classical theory of elasticity. 
Secondary stresses have also been demonstrated in liquids under 
laminar flow conditions by Weissenberg. Their effect is easily 
shown with solutions of high polymers. If a cylinder is rotated in 
such a solution, the liquid will creep up along the rotating cylinder, 
the more so the greater the speed of rotation. A simple liquid 
would be thrown off the rotating cylinder if the centrifugal forces 
were great enough. Such behavior indicates that, besides a shearing 
stress, there is a normal stress parallel to the cylinder axis, which is a 
tensile stress. 

In deriving a theory of secondary stresses, Weissenberg? assumed 
that such liquids have elastic properties. He associated the sec- 
ondary stresses with the shear modulus and the shearing stress 
with the viscosity. Mooney’ treated the problem in a similar manner 
with the difference that he applied his theory of nonlinear ‘‘super- 
elasticity” to the elastic behavior.‘ 

Rivlin’ discarded elastic effects altogether and postulated that the 
stress components on an element of the material depend only on the 
state of flow. He associated the shearing stress with the viscosity 
and the normal stresses with a ‘“‘normal stress coefficient’? which has 
the dimensions of viscosity. A similar treatment was given by Reiner.® 

To test the various theories, Roberts’ measured the normal stresses 
produced in a series of Newtonian and non-Newtonian liquids. He 
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found one normal stress to be a tension. The other two were found 
to be equal. The normal stress difference was proportional to the 
shearing stress, within the experimental error, at values up to 10° 
dynes/cm.’, in agreement with Weissenberg’s theory, but in disagree- 
ment with the theories of Mooney, Rivlin, and Reiner. At greater 
values of the normal stress difference, the results were incompatible 
with all theories. Philippoff* has discussed various theories and 
has partly confirmed Roberts’ results by data obtained on poly- 
isobutylene solutions.°® 

Although these experiments prove the existence of secondary 
stresses, they do not indicate whether their presence is associated 
with the elastic or viscous behavior of liquids. What is known is 
that secondary stresses are absent in deformed isotropic bodies. 
They must, therefore, be associated with anisotropic behavior as a 
result of deformation, irrespective of whether the deformation is 
elastic or viscous. 

Since most industrially important materials are anisotropic or 
become so during deformation, the secondary stresses will, in the 
following, be treated for anisotropic systems in general. The treat- 
ment in three dimensions is confronted with difficulties. For reasons 
of simplification, the case of simple shear will be considered, which 
reduces the problem to a two-dimensional one. Since the relations 
between stress components and strain or strain rate components are 
not known, the motion of an element relative to rotating axes is 
considered. This procedure eliminates strains or strain rates, and 
results in a representation of forces in terms of stress components 
and inertial forces only, which applies to any kind of deformation, 
whether elastic or time-dependent, and can be compared with experi- 
mental results. 


Il. Stress Distribution in Anisotropic Bodies under Plane Shear 


The stress components in an isotropic body under simple shear 
consist of the shearing stress, 7, acting tangentially on the shear 
plane, and an arbitrary hydrostatic pressure, p, normal to the plane 
(Fig. 1). The stress on the plane is 


= SP? f (1) 


By superimposing on this a system of secondary stresses, we obtain 
the stress distribution in an anisotropic body, as shown in Figure 2. 
A circle inscribed in the plane ABCD, with the coordinate axes as 
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Fig. 2. Stress distribution in anisotropic body under simple shear. 


diameters, is changed on deformation into an ellipse, the principal 
axes of which are the shortened and lengthened diameters of the 
criginal circle. Let AO and OB’ coincide with the directions of the 
axes of the ellipse; the secondary stresses can then be imagined to 
arise from a rotation of the coordinate axes. Considering only a 
state of stress, they refer to the axes of the stress ellipsoid. 

In addition to the hydrostatic pressure and shearing stress, we have 
the stresses o, and o, normal to the planes AO and OB’. The latter 
are associated with the principal stresses o, and o; normal to the 
planes bounded by the coordinate axes. The condition of anisotropy 
requires that o, # o, and a; # o;. The oblique planes are further 
subjected to tangential stresses which are equal and are denoted by 
the secondary shearing stress 7;. The resultant of the secondary 
stresses is a stress component ¢ acting normal to the plane AB but 
making an angle a with the direction of OB’. With this term the 
stress for the condition of plane shear becomes 


S&=aortr (2) 


which is similar to eq. (1), but ¢ is a stress deviator. 
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For the correlation of the various stress components, we consider 
the condition of equilibrium of forces in the prism AOB’, and 
¢ cosa — T8iN a = o, cosa + 1 S8iIN a@ 
. (3) 
o sin a + 7 CcOS a = of SIN a + 7, COS @ 
Multiplying the first equation by cos a, the second by sin a, and 
adding the results gives ¢ in terms of the secondary stresses: 


o = o, 8in’a + o, cos* a + 7; sin 2a (4) 


Multiplication of the first eq. (3) by sin a and the second by cos a 
gives, on addition and subtraction of the results, the following expres- 
sions for o and r: 


o = '/, (or + oy) + 7:/sin 2a — 7 cot 2a (5) 

tr = '/, (oz — oy) sin 2a + 7; cos 2a (6) 

In order to express ¢, and ¢, in terms of the principal stresses 

o, and o», which can be experimentally determined, the same opera- 

tion is applied to the equilibrium of forces in the prisms OC’A and 
OD’B’, and we have the following relations: 


o, = —p + '/2 (01 + os) + 7 cot 2a 


oa, = —p + '/2 (01 + ox) — 1 cot 2a 
m7 = —'/, (eo, — a2) sin 2a (7) 
Oz + oy = 1% + O2 — 2p; os — a, = 27; cot 2a 
01 — Oo = (0; + oe) cos 2a 


With these terms eqs. (5) and (6) become 


o = '/2 (01 + o2 — 2p) — 7 cot 4a (8a) 
o'= 0 — p — 7 cot 2a (8b) 
¢=o—p+rtan 2a (8c) 

t = 27, cos 2a = (0; — a2) sin 2a cos 2a (9) 


For the conditions chosen, there is no external pressure or tension 
applied, and the principal stresses can be either independent of the 
hydrostatic pressure, or one of the principal stresses can be equal to 
it. Hence 
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for a. = p: ¢ = 7 cot 2a (10a) 
for o; = p: o = 7 tan 2a (10b) 
o = —r(cosec 4a sec 2a + cot 4a) — p (10¢) 


where the last equation is obtained by combining eq. (8a) with eq. (9) 
and the last of eq. (7). The point O in Figure 2 on rotation describes 
a curve, the slope of which is determined by the trigonometric func- 
tions of eqs. (10). 

Introducing eqs. (10) into eq. (2) gives the stress in terms of the 
shearing stress alone for the various conditions. The equations de- 
scribe a state of plane stress and refer to the undeformed state. 
Strains, however, are associated with the deformed state. Both states 
are identical if the deformations are infinitesimal, or if the body is iso- 
tropic. This condition does not hold for larger strains in anisotropic 
bodies, and for a complete description of the state, the stresses have 
to be related to the deformed state. 


III. Stress in Relation to the Deformed State 


For the anisotropic case, the strain components can be obtained 
in a manner similar to the stress components by replacing the angle 
a by a’, the angle associated with the deformational rotation of the 
coordinate axes. Let ¢, «, €, ¢, and y be the strain components and 
shear associated with ¢, 01, a2, p, and 7, respectively; then the strains 
corresponding to the stresses in eqs. (10) are 


for € = '/3e: «= —y7 cot 2a’ (lla) 
for «, = '/3e: ¢€ = ¥ tan 2a’ (11b) 
for the general case: « = y(cosec 4a’ sec 2a’ + cot 4a) (lle) 


These equations, or their counterparts in terms of strain rates, hold 
only in special cases. From the standpoint of generalization, noth- 
ing is gained by such a procedure, since a given stress component can 
produce, in the extreme case, all of the strain components. In 
addition, the relationships between the stress and strain components 
are not known. This difficulty is overcome by superimposing upon 
the stress system the rotating axes of the strain ellipsoid, in a manner 
similar to that followed for the treatment of the secondary stresses. 
For the solution of this condition, use is made of the equations of 
motion. 
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For the case considered, the shearing stress r acts in the y direction, 
and the normal stress ¢ in the z direction. Since the normal stress 
is a function of the shearing stress, we can generalize by setting ¢ = 
+f(a) and the equation of motion becomes 


dr/dy + f(a)/Ox — pg = pk (12) 


where p = mass per unit volume and # = 0#/0¢t is the acceleration in 
the x direction. The term pg is the body force per unit volume and is 
written here as the product of density and acceleration due to gravity. 
Multiplication by dy gives this equation in terms of forces per unit 
area instead of gradients, and 


Or + Orf(a)/(Ox/Oy) = p(¥ + g)dy (13) 
The derivative 0x/dy arises from the rotation and is, therefore: 
(Ox /Oy)roe = '/2(Ox/Oy — Oy/dz) (14) 


A point (z,y) on rotation describes a curve, the slope of which is 
dx/Oy. Rotating the axes counter-clockwise through an angle 6 = 
90° — a@ (see Fig. 2) gives: 

(Oz/Oy)g = '/2(tan 8 — cot 8) = —cot 28 = cot 2a’ (15) 
During this rotation Oz increases and dy decreases. As the angle 28 
increases from 90° to 180°, cot 26 increases from zero to infinity. This 
is equivalent to a tension in the z direction in the absence of forces 
in the y direction, and corresponds therefore with the condition 
g~-> = 0. 

A clockwise rotation produces a curve which is normal to that of 
eq. (15) and 

(O2/Oy)%+8 = —1/cot 2a’ = —tan 2a’ (16) 
Here Oz decreases as a result of a tension in the y direction, and 
o, — p= 0. Where both rotations are operative, we obtain a curve 
intermediate in position between the two described; hence 
(Or/Oy) = '/2(cot 2a’ — tan 2a’) = cot 4a’ (17) 
yeneralizing by setting Ox/Oy = f(a’), eq. (13) becomes 
Or + Orf(a)/f(a’) = p(% + g)Oy = plz + g)dx/f(a’) 
Both functions f(a) and f(a’) vary with stress and rotation. With 
S(a)/f(a’) being constant, integration gives 
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t + th(a)/f(a’) = fo(% + g)dx/f(a’) = F, (18) 


This equation describes the dynamics of a particle relative to rotating 
axes. We have then the shearing stress 7, the force per unit area 
acting tangentially on the original frame of reference, and a force, 
described by the integral, which is the sum of the inertial force and 
body force (or combined inertial force, for short) per unit area with 
accelerations relative to the rotated coordinates. The difference 
between these forces is the apparent force rf(a)/f(a’) = «/f(a’), 
which also acts on the rotated frame. 
Combining eqs. (15-17) with eqs. (10) gives 


Il 


for o2 = p: tT — 7 cot 2a/cot 2a’= F;, 


for o, = p: rt — 7 tan 2a/tan 2a’= F, (19) 
for the general case: 
t — 7 (cosec 4a sec 2a + cot 4a)/cot 4a’ = F, 


On substituting the generalized function f(a)/f(a’) by 6, these 
equations reduce to 


r = br + (1 — b)r (20) 


where (1 — b)r is the combined inertial force per unit area. This 
equation holds also if b is replaced by a variable b’ varying with the 
stress. 

Equations (19) and (20) represent the equilibrium between all the 
forces acting on a system and are, therefore, applicable to any kind 
of deformation whether elastic or time-dependent. As such, they 
are incomplete, in that they do not refer to stress-strain relations. 
However, since they contain secondary effects resulting from the rota- 
tion of the cordinate axes, they can be applied to the tangential 
components of stress and deformation which are obtained from experi- 
ment. This procedure will be described for non-Newtonian flow. 


IV. Non-Newtonian Viscosity and Anisotropic Behavior 


Viscosity measurements are often carried out so that the tangential 
shearing stress is known and the resulting rate of deformation is 
measured, or by applying a given rate of deformation and measuring 
the resulting shearing stress. The direction of this rate of deforma- 
tion coincides with that of the tangential shearing stress, and the 
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Fig. 3. Relationship between shearing stress and shear rate for non-Newtonian 
flow. 


rate of deformation is proportional to the shear rate in the same frame 
of coordinates. 

Non-Newtonian viscosity varies with the shearing stress or with the 
shear rate. This variable viscosity is often expressed as an ‘‘apparent 
viscosity,” shearing stress over shear rate, in analogy with Newtonian 
viscosity. Actually, non-Newtonian viscosity is determined by the 
tangent to a point of the shearing stress-shear rate curve (Fig. 3). 
The tangent intersects the shearing stress axis at point B dividing 
the stress into two parts, an apparent stress AB and an intercept OA. 
These parts are, according to eq. (20), the apparent stress br and the 
intercept (1 — b)r. The non-Newtonian viscosity is then 


9’ = dr/dy = br/y (21) 


Substituting for b the trigonometric functions from eqs. (19a) and 
(19b) gives 


for o: = p: n’ = 7 cot 2a/¥7 cot 2a’ (22) 


Ul 


for 0, = p: n’ = 7 tan 2a/y7 tan 2a’ (23) 
These relationships show that non-Newtonian viscosity, in contrast 
to Newtonian viscosity, is determined solely by the secondary stress 
components, r cot 2a and 7 tan 2a and the secondary strain rate com- 
ponents, 7 cot 2a’ and ¥ tan 2a’. 

It will be seen that these strain rates correspond to the secondary 
strains given in eqs. (lla and 11b), and 7’ = a/é. It follows also 
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that for cot 2a < cot 2a’ and tan 2a < tan 2a’, b is smaller than 1, 
and the shearing stress-shear rate curve is concave towards the shear 
rate axis. In the reverse case, b is larger than 1, and the curve is 
convex toward the abscissa. For the general case we have from 
eqs. (19) 

n’ = dr/dy = r(cosec 4a sec 2a + cot 4a)/y cot 4a’ (24) 


It can be reasonably assumed that at low shear rates the directions 
of stress and deformation coincide. In this case a = a’, and the 
viscosity according to eqs. (22 and 23) will be constant, which has 
been verified by experiment (e.g., reference 9). This viscosity is, 
however, not Newtonian in the strict sense of the word, since it is 
stili the result of secondary stresses and strain rates, and may be 
termed ‘‘equilibrium viscosity.”” A liquid with such a viscosity still 
shows the Weissenberg effect, which would be absent in a Newtonian 
isotropic liquid. 

The three equations for non-Newtonian viscosity represent three 
different kinds of behavior under shear. In the case of eq. (22), 
the secondary stress effect results in a positive pressure normal to the 
direction of shear, and the Weissenberg effect is positive. Where 
eq. (23) holds, this pressure acts in the direction of shear and the 
Weissenberg effect is negative. This effect works in both directions 
for the general case of eq. (24) and is more difficult to demonstrate. 

Integration of eq. (21) with b = constant gives 


t/t> = (¥/%)" (25a) 


where tr, and 7 refer to the upper limits of the region with equilibrium 
viscosity. In the absence of this limit, the equation may be written 
as follows: 

r= By (25b) 
where Bisaconstant. This equation has been empirically derived by 
de Waele,” Farrow and Lowe," and W. Ostwald.'” 

Although 7 is a tangential shear rate, it is not the same as for 
Newtonian flow. However, with the aid of eq. (25a), the velocity 
gradient can be correlated with the stress gradient for any experi- 
mental shear condition. 

An interesting case with a variable value of b is obtained when the 
product of shearing stress and the variable b is constant: b’r = bro. 
The coefficient of non-Newtonian viscosity is then 
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n’ = dr/dy = b’r/y = bri/¥ 


and is a function of the shear rate alone but independent of the stress. 
Integration yields 


exp[(r — t0)/t0] = (¥/¥0)” (26) 
which has also been given by Ludwik."* 

If the volume of a non-Newtonian liquid remains constant during 
the shearing operation, the viscosity remains a shear viscosity. 
However, when accompanied by volume changes, non-Newtonian 
viscosity can be divided into a shear viscosity associated with the 
secondary shearing stress 7; and a volume viscosity in terms of the 
principal stresses in accordance with eqs. (8) and (9). 

For the evaluation of the trigonometric functions of a and a’, 
it is necessary to measure the normal stress which according to eq. 
(10) is o = rf(a@). With the knowledge of a and b, which is obtained 
from the flow diagram, the angle a’ can be calculated. Since a’ is 
associated with the deformation, its knowledge permits the calcula- 
tion of the flow birefringence (see Peterlin"*). 

The association of non-Newtonian viscosity with anisotropic be- 
havior leads to a new interpretation of the work function. The area 
under the curve of Figure 3 represents the work per unit volume per 
second, and for the case where the shearing stress is a power function 
of the shear rate according to eq. (25), is 


1 
1+6 


The shearing stress is the sum of two components as shown by eq. 
(20), and the work function can be also divided into two parts 


W' =JSrdy = 





TY (27) 








ry +S p(% + g)dxdy/f(a’) 


(28) 


The first term on the right side of the equation represents the work 
per unit volume per second due to the secondary stresses, and the 
second term, containing the inertial and body forces, is dissipated 
into heat. For the Newtonian case, this work is W’ = '/,ry, and is 
completely dissipated. 

Similarly for non-Hookean elasticity, the work per unit volume is 


ain, 
ee tes 
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2 b 
i+” 148 
For the case of Hookean elasticity this work is W = '/,ry and repre- 
sents only the potential energy. The inclusion of the kinetic energy 


in the work for non-Hookean elasticity indicates that the latter can 
be associated with time effects. 


W= 


y= ry +S p(% + g)dady/f(a’) (28) 


V. Stress Distribution in Anisotropic Bodies under Tension and 
Compression 


For the condition of tension or compression, the stress distribution 
as shown in Figure 2 is also valid after omitting the shearing stress 
rt. With o; being the principal tensile or compressive stress, the 
stress becomes in analogy with eq. (2): 


& =ate (29) 


The normal stress component oa is defined by eq. (4). Omitting the 
hydrostatic pressure, since it does not contribute to the deformation, 
we have from eqs. (7) 


o = '/(0; + a2) + '/2(01 — o:) cos 4a (30) 
o = o; sin? 2a + a; cos* 2a , 
From eqs. (7) we have in addition 


0; — o2 = (0; + a2) cos 2a 
hence 


o: = o; tan’a 
which simplifies the second of eqs. (30) as follows: 
o = o;(cos? 2a + sin? 2a tan’? a) = o,(1 — sin4atana) (31) 


The principal stress o; acts in the z direction and the normal com- 
ponent in the y direction, and eq. (18) becomes for the case of com- 
pression or tension with (1 — sin 4a tan a)d0x/Oy = constant: 


o; — o(1 — sin 4a@ tan a) cot 4a’ = f p(% + g)dy cot 4a’ (32) 


The coefficient of traction, corresponding to non-Newtonian 
viscosity, is consequently 


up’ = do;/dé = o;(1 — sin 4a tan a) cot 4a’/é (33) 
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Since for a given value of o, the trigonometric functions of a and a’ 
have different values for tension and compression, the results show 
that non-Newtonian liquids have different rheological properties 
under shear, tension, and compression, respectively, as should be the 
case for anisotropic bodies. 
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Synopsis 
Secondary stresses have been the subject of several investigations since their 
presence was first demonstrated by Weissenberg in liquids during laminar flow. 
They have been associated with elasticity by some investigators, and with viscous 
effects by others. In these theories, secondary stresses are assumed to arise from 
the deformation or rate of deformation in otherwise isotropic bodies. Since most 
industrially important materials are not isotropic, a theory of secondary stresses 
is given for anisotropic behavior in general. It is shown that non-Newtonian 
viscosity is the direct outcome of secondary stresses, and their importance for 
rheological behavior is discussed. In relating secondary stresses with non- 
Hookean elasticity, it is found that this type of elasticity is also associated with 
time effects. 
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Volume Changes and Dewetting in Glass Bead— 
Polyvinyl Chloride Elastomeric Composites 
Under Large Deformations* 


THOR L. SMITH, f Jet Propulsion Laboratory, California Institute of 
Technology, Pasadena, California 


INTRODUCTION 


Elastomeric materials loaded with inert filler have rheological 
properties which, on the one hand, are qualitatively like those of 
unfilled polymeric materials, but which, on the other hand, are unique 
to filled systems. Their stress-strain curves usually have several 
distinct regions which are similar to those for certain unfilled poly- 
meric materials. The curves have an initial region in which the stress 
increases with strain, and here the modulus depends primarily on the 
amount of filler and the properties of the binder. In the second 
region, the stress may be roughly independent of the strain or may 
increase at a rate markedly less than that in the initial region. The 
transition between these two regions may be either gradual and hardly 
perceptible or sharp, with a well-defined yield point. Thus, in these 
two regions, the stress-strain curves may be like those for typical 
unfilled rubbers or for plastics, or they may resemble the curves for 
rubbers at low temperatures, which cold-draw or neck. In the third 
region the stress increases more rapidly with strain than in the second 
region and continues to increase until the specimen breaks; again, this 
is similar to what is observed for certain plastics and rubbers. Thus, 
both filled and unfilled systems can have sigmoidal stress-strain curves 
which may or may not have distinct yield points. 


* This paper presents the results of one phase of research carried out at the Jet 
Propulsion Laboratory, California Institute of Technology, under sponsorship of 
the Department of the Army, Ordnance Corps (Contract No. DA-04-495-Ord 18). 

t Present address: Stanford Research Institute, Menlo Park, California. 


113 











114 T. L. SMITH 


However, filled systems, and especially those in which the binder 
wets, but does not react chemically with the filler, exhibit the sig- 
moidal-type stress-strain curve for different reasons than do rubbers 
and plastics. In filled systems, the binder adheres to the filler 
initially, which reinforces or strengthens the system. However, in 
the transition and second regions of the stress-strain curve, adhesive 
bonds between the filler and binder break, a process called dewetting, 
and vacuoles form around the filler particles. After the adhesive 
bonds are essentially all broken, a further increase in strain only 
stretches the binder further and enlarges the vacuoles. 

Perhaps the best method for studying dewetting is to measure the 
volume of the filled system as a function of elongation. When an un- 
filled rubber is stretched, its volume remains essentially constant, 
provided crystallization does not occur. However, when a &lied sys- 
tem is stretched, and vacuoles form, the volume increase may be 
large: e.g., 50 to 100%, or possibly more. The measurement of this 
volume increase by means of a sensitive dilatometer can give a precise 
value for the elongation at which the dewetting begins, and possibly 
that at which it is complete. 

Additional information can be derived from a study of the volume 
changes. Since the data show precisely the elongation at which ad- 
hesive failure begins, methods for modifying the adhesive strength can 
be investigated. Further, the temperature dependence of the adhe- 
sive strength can be studied by carrying out tests at different tempera- 
tures. If certain assumptions are made, an estimate can be obtained 
of the elongation of the binder between adjacent filler particles at the 
time when dewetting occurs. This result, when combined with the 
stress-strain properties of the unfilled binder, gives an indication of 
the stress required to break the adhesive bonds. 

Relatively few dilatometric studies have been made of the volume 
changes that occur during the stretching of filled materials. Holt 
and McPherson! made a careful study of the stretching of unfilled 
vulcanized rubber and the volume changes resulting from crystalliza- 
tion; also, a few studies were made of a gum stock filled with whiting 
and a carbon black. These investigations were extended by Jones 
and Yiengst,? who studied the volume changes that occurred when 
various pigmented-rubber systems were stretched in a dilatometer 
similar to one used by Holt and McPherson. More recently, some 
data were obtained in a similar manner by Bryant and Bisset* for 
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several filled rubber vulcanizates. Prior studies relating to vacuole 
formation in filled rubbers were made by Schippel,‘ Green,’ Vogt and 
Evans,® and DePew and Easley.’ 

The present study was made to extend previous work and to obtain 
data at several temperatures. Filled systems composed of small glass 
beads dispersed in a polyvinyl chloride—dioctyl sebacate rubber were 
stretched in a dilatometer at 0,25, and 50° C. Less precise data were 
also obtained by stretching specimens and measuring the change in 
their width and length. A simple equation was derived to relate the 
volume fraction of filler to the yield strain obtained from the dilato- 
metric data. 


EXPERIMENTAL 
A. Materials 


(1) Glass Beads 


The glass beads were obtained from the Minnesota Mining and 
Manufacturing Company; their No. 1160 beads were used in nearly 
all the composites studied. These beads were essentially spherical 
and were analyzed for size and polydispersity by means of a Micro- 
merigraph. The beads had a weight-average diameter of approxi- 
mately 40 u, and about 94 wt.-% of the beads had diameters between 
20 and 60 u. 

Limited studies were made with beads Nos. 118, 116, 114, and 112. 
From rough Micromerigraph data, it appeared that these beads had 
weight-average diameters of about 50, 80, 105, and 140 u, respectively. 
(The size of the No. 112 beads was not measured, but a diameter of 
140 u was assumed after comparing data supplied by the manufacturer 
with measurements obtained on the Micromerigraph for beads Nos. 
118, 116, and 114.) 


(2) Elastomer 


The elastomer was prepared from equal weights of a polyvinyl 
chloride plastisol (Geon 121) and diocty] sebacate and contained 2.5% 
lead stearate as stabilizer. To obiain a uniform dispersion, the solids 
were slowly added to the dioctyl sebacate while it was stirred mildly 
by a Brookfield counter-rotating stirrer. The mixture was then 
vigorously stirred for 10 min. by the Brookfield stirrer running at a 
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high speed. The resulting dispersion was cast into trays, degassed 
for about 30 min. (or until no more gas was evolved), and then cured 
for about 1 hr. at 180° C. The cured sheets, about '/; in. thick, were 
visibly homogeneous and completely bubble-free. 


(3) Composites 
In preparing the glass-bead composites, the polyvinyl chloride 
dispersion was mixed as described above, and then degassed. The 
glass beads were next added slowly while the dispersion was stirred at 


low speed. The resulting mixture was cast into trays, degassed, and 
then cured. 
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Fig. 1. Elongated-ring specimen used in dilatometer. 


In addition to sheets of the composites, elongated ringlike speci- 
mens were prepared for testing in the dilatometer. The shape and 
dimensions of the rings are shown in Figure 1. Specimens of this 
shape were chosen in preference to circular rings, since they could be 
stretched in the dilatometer without first being flattened, which might 
disrupt some adhesive bonds between the beads and binder. The 
rings were prepared by casting the uncured composite into an open 
mold machined in a small block of aluminum. After the dispersion 
was poured into the mold, it was degassed again; the excess was then 
seraped from the surface of the mold, and the mold was placed on a 
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hot plate to effect cure. Highly uniform rings, which ordinarily 
differed in weight by less than 5%, were prepared in this way. 

One disadvantage of rings is that the strain in a stretched specimen 
is not uniform, since the outside of the ring is elongated more than the 
inside. The rings used had an outside-to-inside-circumference ratio 
of about 1.23; therefore, the average circumference was used to cal- 
culate an average strain or extension ratio from the extension of a 
ring. 

The cured composites were essentially free of voids, as determined 
by density measurements, and the beads appeared to be uniformly 
dispersed. However, composites prepared from beads Nos. 112 and 
116 had a thin resin-rich surface, which indicated that these beads had 
settled somewhat after mixing. 

The densities of representative samples of the composites prepared 
with No. 1160 beads were obtained by hydrostatic weighing. These 
values are compared in Table I with the theoretical densities cal- 
culated from density values of 2.39 and 1.135 g./ml. for the beads and 
the elastomeric binder, respectively. The calculated and measured 


TABLE I 
Experimental and Calculated Densities of Composites 
Containing No. 1160 Glass Beads 











Glass- Glass- 
bead bead Density, g./ml. 

content, content, 

wt.-% vol.-% Experimental Calculated Difference 
25 13.7 1.273 1.308 —0.035 
50 32.2 1.355 1.541 —0.186 
60 41.5 1.645 1.657 —0.012 
7 52.5 1.840 1.794 0.046 
75 58.9 1.895 1.875 0.020 
80 65.4 1.975 1.955 0.020 





densities agreed within a few per cent, except for the composite having 
50 wt.-% beads. The reason for this disagreement is not known; it 
probably is the result of some error, since no other evidence was ob- 
tained which indicated voids in any samples for which data are re- 
ported. (Densities were measured on only one set of samples, and 
many samples were prepared during the study.) 
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(4) Polyurethane Foam Rubber 


To prepare the foam, a castable polyurethane elastomer filled with 
59% by volume of a water-soluble inorganic salt was first prepared. 
About 70% of the salt had an average particle diameter of 150 uy, 
and 30% had an average diameter of 154. The foam was obtained by 
extracting the salt from the filled elastomer with boiling water. The 
foam was dried carefully and chemically analyzed to verify that the 
salt had been extracted completely. 


B. Apparatus and Procedure 
(1) Dilatometer 


A dilatometer was constructed which was similar to the one used by 
Holt and McPherson! and by Jones and Yiengst.? It consisted of a 
metal tube which had rubber O-rings and screw caps on eachend. A 
piece of drill rod passed completely through the metal tube and screw 
caps. Inserted in each screw cap was a threaded stuffing box, which 
prevented leakage of the dilatometer fluid and yet allowed the rod to 
be moved up and down. The metal tube had a side arm, to which a 
calibrated glass capillary and an auxiliary reservoir were attached. 
To hold the ring test specimen in the dilatometer, a metal hook was 
attached on the inside, near the top of the metal tube, and a similar 
hook was attached to the drill rod. 

In making measurements, the top screw cap was removed, and a 
ring specimen was hung over the hook attached to the inside of the 
metal tube. The screw cap was then replaced, and the dilatometer. 
was filled with water from the auxiliary reservoir. A small ver t hole 
in the upper screw cap permitted the escape of any trapped air, after 
which the vent hole was plugged, and the auxiliary reservoir was iso- 
lated from the main part of the apparatus. To test for leaks, the 
water level in the capillary was observed while the drill rod was moved 
down and up. The drill rod was next moved to the appropriate posi- 
tion and then rotated until its attached hook moved inside the ring 
and was in the same vertical plane as the upper hook. The rod was 
slowly moved down until fiducial marks on the exposed portion of the 
rod and on a reference rod indicated that the ring was snugly fitted 
between the hooks and ready to be tested. Next, the rod was moved 
down stepwise to stretch the ring. After each step, the water level 
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in the capillary was read, and the position of the fiducial mark on the 
rod was measured by a cathetometer. This procedure was continued 
until the ring had been stretched about 200%. 

Some data were obtained with another dilatometer, which was 
similar except that the tube was made of Lucite tubing, rather than 
metal. In this dilatometer, the extension of the ring could be ob- 
served visually. Occasionally, the extension of a ring was measured 
directly by the cathetometer, and the results agreed closely with values 
obtained from the fiducial mark on the rod. However after two 
Lucite tubes had cracked, the tube was replaced by a metal one. 

To obtain data at 0, 25, and 50° C., the dilatometer was placed in 
a water bath controlled to +0.1° C. or less. A few runs were made 
at —19° C.; in these, a silicone fluid was substituted for water as the 
confining liquid in the dilatometer. Accurate temperature control 
of the thermostat was always maintained in order to prevent small 
temperature fluctuations during an experiment from changing the 
water level in the capillary. 


(2) Dimension Measurements 


Some studies were made by measuring the elongation and width 
of rectangular tensile specimens which were stretched stepwise. 
Bench marks were inked on a specimen, and the separation of these 
marks was measured by a cathetometer after each incremental elonga- 
tion of the specimen. The width of the specimens was measured by a 
micrometer. Although this simple procedure is not as precise as the 
dilatometric method, it provided reasonably accurate data. 


THEORETICAL 
A. Poisson’s Ratio 


The deformation of an elastic material in pure homogeneous strain 
may be characterized by three principal extension ratios. For 
example, suppose that an undeformed rectangular parallelepiped, 
having length Lo, width Wo, and thickness 7, is deformed by loads 
applied perpendicularly to its opposite faces. The resulting deforma- 
tion is given by the principal extension ratios \y = L/Lo, \» = W/Wo, 
and As; = 7'/7>, where L, W, and T are the new length, width, and 
thickness. Thus, the volume change is given by the relation 
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V/Vo _ ArA2A3 (1) 


where V and V, are the volumes of the deformed and undeformed 
material. 

When an isotropic material is subjected to tension or compression, 
Ae equals A3. If no volume change occurs, the following relation 
holds: 


1 = AA” (2) 
This leads to the equation 
dd2/dr, = —2/2M1 = dew/de (3) 


where «¢,, is the lateral strain, and « is the longitudinal strain. Ac- 
cording to classical elasticity, based on small-deformation theory, 
Poisson’s ratio vy equals —de,/de,, and eq. (3) shows that this ratio 
varies with the extension, even though it is assumed that no volume 
change occurs. . However, eq. (3) can be integrated to give 


In d2/In \ = —1/2 (4) 


where In A, and In \; are the Hencky, or logarithmic, lateral and 
longitudinal strains. Thus, Poisson’s ratio, defined in terms of the 
logarithmic measure of strain, is independent of the extension and 
equals '/», provided that the volume remains constant. 

Poisson’s ratio, defined in terms of logarithmic strain,* is a con- 
venient number to use in specifying volume changes that occur when 
an elastic material is subjected to large deformations. The volume 
change that occurs during a tensile or compressive deformation is 
given by the following equation, derived by differentiating eq. (1) 
after noting that \, equals As: 


d In V/Vo _ +2 din rs _ _ 5 
diny “dink om ” 





where Poisson’s ratio v, which may be a function of \,, is defined as 
follows: 


vy = —dind/dindy (6) 


* Certain valid reasons are recognized for not defining Poisson’s ratio in terms 
of logarithmic strain. However, Poisson’s ratio defined in terms of logarithmic 
strain is used in this paper only as a convenient index for rate of volume change, 
even though some other name should, perhaps, be given to this quantity. 
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B. Yield Strain 


Certain filled elastomers, like glass bead—polyvinyl chloride com- 
posites, can be stretched to a limited extent without showing any 
appreciable volume change. However, at a reasonably well-defined 
yield strain, dewetting occurs, and the volume of the sample increases 
continuously with further stretching. The yield strain can be related 
to the volume fraction of filler by an equation based on certain sim- 
plifying assumptions. 

Consider a unit cube of an elastomer filled with n’ spherical par- 
ticles, each having a radius r and a modulus which is large compared 
with that of the elastomeric binder. Assume that these spheres are 
uniformerly dispered, and that their centers lie in a cubic close- 
packed array: i.e., if the radius of each sphere were enlarged until all 
spheres touched, a cubic close-packed array would result. 

Now consider a line through the centers of the n spheres in one row. 
The length of the line is unity, the fraction of the line occupied by 
spheres is 2rn, and the fraction occupied by binder is d. Thus, 


2rn + dy = 1 (7) 
The volume fraction of spheres V,; can be shown, by considering the 
geometrical properties of cubic close-packed arrays, to be 


44/2 
= V2 (on)? (8) 





Vy 


By combining eqs. (7) and (8), the following expression is obtained: 
34/2 /s 
(2%) V," + do = } (9) 


If the cube containing the filler is stretched a distance Al, and if de- 
wetting does not occur, then the following equation can be written: 


1.105V;"* + db + Al=1+ Al (10) 


where 1.105 = (34/2/)'". Equation (10) can be rearranged to 
give 


1105V/;7+(lit+te)d=1lte (11) 


where e¢, is the strain in the binder between adjacent spheres and equals 
Al/do, and « is the strain imposed on the unit cube and equals Al. 
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Now, by substituting in eq. (11) the value of dy given by eq. (9), the 
following relation is obtained: 


e = «, (1 — 1.105V,/”) (12) 


This equation gives ¢/e, (the ratio of the over-all strain imposed on the 
sample to the strain in the binder along the line between the centers 
of adjacent spheres) in terms of the volume fraction of spheres. Now, 
assume that dewetting occurs when ¢«, reaches a critical value ¢,,., 
which is assumed to be independent of the amount of filler. Thus, 
the observed yield strain «, is given by the equation 


te = é¢ (1 — 1.105V,”) (13) 


When no filler is present, this equation shows that «, = «,,., and, when 
V, equals 0.74 (the volume fraction of spheres in a cubic close-packed 
array), that «. = 0. 


RESULTS AND DISCUSSION 
A. Preliminary Studies 


A tensile specimen of the polyurethane foam rubber was stretched 
in incremental steps, and the width and per cent elongation were 
measured for elongations up to about 50%. Figure 2 shows that -log 
W/W, is directly proportional to log L/L, and that Poisson’s ratio 
is 0.362. When similar measurements were made on a polyurethane 
elastomer, Poisson’s ratio was found to be 0.50 over the entire range 
of extensions covered. The behavior of the foam is like that of a 
filled elastomer after dewetting is complete, in that Poisson’s ratio is 
independent of extension. 

Similar measurements were made on specimens of the polyvinyl] 
chloride composites containing 13.7, 32.2, 41.5, and 58.9 vol.-% glass 
beads. (These volume per cents, and others given in this paper, were 
calculated from the known weight per cents of binder and glass beads.) 
The data from these preliminary tests are shown in Figures 3 and 4; 
data on the unfilled polyviny! chloride elastomer are also included in 
Figure 3. In small deformations, Poisson’s ratio for all polyvinyl] 
chloride composites is 0.50, within the experimental error, and, for the 
unfilled elastomer, it is 0.50 over the entire range of extensions (i.e., 
0 to 200%). For the composites, however, dewetting begins at some 
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Fig. 2. Variation in width and length during extension of foam-rubber specimen. 
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Fig. 3. Variation in width and length during extension of glass-bead composites. 
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Fig. 4. Variation in width and length during extension of composites containing 
58.9% glass beads. (Legend on first section of curve should read: » = 0.500.) 


critical extension; after a transition region, -log W/W» again increases 
linearly with log L/L, but at a reduced rate. The slope of the final 
linear portion of a curve equals Poisson’s ratio, which is seen to in- 
crease from 0.24 to 0.50 as the glass content is decreased. 

After a sample of the composite containing 41.5% glass had been 
stretched 200%, it was allowed to recover to within 3% of its initial 
length, and the measurements were then repeated. Poisson’s ratio 
was now found to be 0.24 over the entire extension range. Although 
this result does not definitely prove that the dewetting was complete 
during the first extension, it does show that the adhesive bonds did 
not re-form when the sample recovered. However, some qualita- 
tive observations indicated that adhesive bonds can re-form in a re- 
covered specimen, provided a sufficiently long period elapses. 

Some data, shown in Figure 5, were also obtained on composites 
containing 52.5 vol.-% of Nos. 112, 114, 116, and 118 glass beads. 
Again, Poisson’s ratio was constant at both small and large deforma- 
tions. The value was approximately 0.50 in small deformations and 
0.172 in large deformations, even though the composites contained 
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Fig. 5. Variation in width and length during extension of composities containing 
52.5% glass beads of different sizes. 


different-sized beads. The composites had different yield values, 
although the values obtained are not very accurate. In general, 
yield values are not very reproducible, but Poisson’s ratio in large 
deformations is usually reproducible. 


B. Dilatometric Studies 


More precise and extensive data on Poisson’s ratio and the yield 
strain were obtained by use of the dilatometer. First, an unfilled 
polyvinyl chloride ring was elongated to 150% in incremental steps in 
the dilatometer. The volume increase observed was exceedingly 
small, approximately 0.23%, and gave a Poisson’s ratio of about 0.499. 
This value is in close agreement with the value of 0.500 reported by 
Holt and McPherson! for vulcanized natural rubber at extensions up 
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Fig. 8. Volume increase during extension of composites containing 41.5% glass 
beads. 


to 300%. Other experiments made on the unfilled polyviny! chloride 
elastomer, stretched to 200% and more, also showed essentially no 
volume changes. It is thus concluded that, when the binder itself 
is extended, volume changes do not occur (at least not at 25° C.) 
either because of formation of voids or resulting from crystallization. 

Volume changes as a function of extension were measured dilato- 
metrically for the composites containing 13.7, 32.2, 41.5, and 52.5 
vol.-% glass. Measurements were made at 0, 25, and 50° C. for the 
composites containing 32.2% glass; the composites having 13.7 and 
52.5% glass were tested at only two temperatures; and the composite 
containing 41.5% glass was tested at — 19° C., in addition to the three 
other temperatures. Ordinarily, duplicate experiments were made, 
but only the results of single experiments are shown in Figures 6 to 9, 
plotted as log V/V» vs. log L/L». 

These figures show that the data from any experiment can be di- 
vided into three sections, similar to those in Figures 3 to 5. In the 
first section, which is for small deformations, little or no volume in- 
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Fig. 9. Volume increase during extension of composites containing 52.5% glass 
beads. 


crease occurs as a specimen is extended. The next section is the 
transition region, in which dewetting of the filler occurs, and the vol- 
ume increases at a progressively increasing rate as the sample is ex- 
tended. In the final section, the dewetting is believed to be complete 
and log V/V» increases linearly with the logarithmic strain. Figure 9 
shows that the rate of volume increase drops off at high extensions. 
This behavior is believed to result from water in the dilatometer diffus- 
ing into the vacuoles or other defects in the stretched specimen. 

By the use of eq. (5), Poisson’s ratio was calculated from the slope 
of the linear section of a curve. The linear section was also extra- 
polated to zero volume change to obtain a yield strain. Values for 
Poisson’s ratio and the yield strain are presented in Table I], together 
with the values for Poisson’s ratio which were obtained in the pre- 
liminary experiments by measuring the width and length of a speci- 
men. At —19° C., Poisson’s ratio and the yield strain for the 41.5%- 
glass composite were 0.259 and 0.08, respectively. The data in 
Table II show that the duplicate values for Poisson’s ratio are in close 
agreement, and that the values obtained from the preliminary experi- 
ments agree closely with those obtained dilatometrically, except for 
the composite containing 52.5% glass. Within the experimental 
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TABLE I 
Poisson's Ratio and Yield Strain for Glass-Bead Composites 








Glass-bead content 

















“aero gige Poisson’s ratio » Yield strain «, 
Bead wt.- vol.- squperevinnguiiasittitaniipamimtane _ jenaingumianeains 
No. % &% 0°C. 26°C. 50°C. ~25°C.* 0°C. 25°C. 50°C. 
1160 25 13.7 0.435 0.456 0.47 0.44 0.58 
0.432 0.463 0.45 0.75 
0.419 0.68 
1160 50 32.2 0.330 0.325 0.344 0.34 0.24 0.46 0.73 
0.326 0.3138 — 0.25 0.57 
. 0.310 0.50 . 
1160 60 41.5 0.234 0.268 0.265 0.24 0.24 0.23 0.49 
0.251 0.257 0.275 0.19 0.31 0.35 
1160 70 52.5 0.259 0.254 0.16 — 0.17 0.30 
0.265 0.09 
0.256 0.12 
1160 75 58.9 0.19 0.10% 
112 70 52.5 0.172 0.08» 
114 70 52.5 0.172 0.10» 
116 70 52.5 0.172 0.12» 
118 70 52.5 0.172 0.06(?) 





* Preliminary data obtained by measuring the width and length of stretched 
specimens. 
» Yield strains obtained from interpolated ‘‘break-away”’ points in Figs. 4 and 5. 


error, Poisson’s ratio is probably independent of temperature, al- 
though the data suggest that it may increase slightly with tempera- 
ture. Although the data on a single ring specimen are very precise, 
the reproducibility of results on identically prepared ring specimens is 
not accurately known. Thus, it is believed that the variation in 
Poisson’s ratio may be partly due to slight variations in specimen den- 
sity, which would indicate the existence of microscopic voids in the 
specimens, and partly due to nonuniform dispersion of the beads 
in the different specimens. The small experimental error in ob- 
taining the strain in a ring, partly due to the use of fiducial marks out- 
side the dilatometer, would not affect Poisson’s ratio, although it 
would lead to error in the yield strain. 

Table II shows that the yield strain depends rather strongly on 
temperature, although its reproducibility is poor. The poor re- 
producibility is not unexpected, since dewetting is a sensitive and 
highly rate-dependent process. Small differences in the microscopic 
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Fig. 10. Variation of yield strain with volume fraction of glass beads. 


structure or mechanical treatment of a specimen can easily cause 
marked variations of the strain range in which dewetting occurs. 
However, the temperature dependence of the yield strain is expected, 
since the modulus of the binder is temperature-dependent. 

To show how accurately eq. (13) predicts the variation of the 
yield strain with the volume fraction of filler, average values of the 
yield strain «, are plotted against 1 — 1.105V,;” in Figure 10. This 
figure also shows yield strains at 25° C. obtained from stress-strain 
curves measured on the Instron tester at a strain rate of 0.6 min.—'. 
The stress-strain curves had two near-linear portions: at strains 
somewhat less and strains somewhat greater than the yield strain. 
These linear portions were connected either by a curve of gradually 
decreasing slope or by a curve which passed through a maximum; the 
latter type of curve was like those for plastics which cold-draw. From 
both types of curves, the yield strain was obtained by extrapolating 
both linear portions until they intersected. 

Although the points in Figure 10 scatter considerably, the lines 
drawn are reasonably consistent with the data and have slopes of 








VOLUME CHANGES AND DEWETTING IN GLASS BEADS 131 


1.05, 1.69, and 2.77 at 0, 25, and 50° C., respectively. These slopes, 
according to eq. (13), should equal «¢,,,, the strain of the binder be- 
tween adjacent particles immediately preceding adhesion failure. 
From the stress-strain curves for the unfilled binder, it was found that 
these strains correspond to stresses of 550, 470, and approximately 
430 p.s.i. at 0, 25, and 50° C., respectively. Thus, the adhesion force 
appears to show little or no temperature dependence. 


C. Stress-Strain Behavior 


Only limited stress-strain data were obtained on the composites 
and unfilled binder. However, some data measured on the Instron 
tester at 25° C. are given in Table III. This table gives the yield 


TABLE III 
Tensile Data for Glass-Bead Composites at 25°C. 








Glass-bead 
content, o.; Ey, E,, E(V, - 
vol.-% p.s.i. & p.s.i. p.s.i. 0.26)~! 
13.7 230 0.50 450 285 473 
32.2 240 0.40 580 140 334 
41.5 146° 0.26" 560 114 351 
§2.5 90" 0.11* 820 94 437 
58.9 63" 0.06" 1050 76 503 





® Obtained by double extrapolation, as explained in text. 


stress o,, the yield strain ¢,, the initial modulus Zo, and the modulus 
after yield E,; the latter modulus was taken as the slope of the linear 
portion of the o—e curves after yield. (Values of o, and ¢, for the three 
composites having the most filler were obtained by interpolation, as 
explained previously.) It is seen that the yield stress and strain de- 
creases, and the initial modulus increases, with filler content; how- 
ever, the modulus does not increase very rapidly, considering the 
large variation in filler content. 

The modulus after yield decreases with filler content, as expected, 
since now the filler only dilutes the binder and has no reinforcing 
effect. If it is assumed that the filler only dilutes the binder, it seems 
reasonable to suggest that E, is proportional to V, — 0.26, where V, 
is the volume fraction of binder. The factor 0.26 is introduced, since 
if the spheres were monodisperse, the maximum volume fraction of 
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spheres would be 0.74, and thus 0.26 wouid be the minimum fraction 
of binder that could be present in a well-consolidated composite. 
The assumption that EZ, is proportional to V, — 0.26 means that EF, 
is zero when V, equals 0.26, which seems reasonable. With these 
several assumptions, the following equation results: 


E, = 1.35E,(V, — 0.26) (14) 


where E, is the modulus of the unfilled binder. Values of £,(V, — 
0.26) —', given in the last column of Table ITI, are seen to be relatively 
constant: none vary by more than 20% from the average, which is 
420. When this average value is divided by 1.35, a value of 310 p.s.i. 
is obtained for the modulus of the unfilled binder, which compares 
favorably with an experimental value of 386. 

Tensile data for the composite having 32.2% filler and the modulus 
for the unfilled binder, as measured at a strain rate of 0.6 min.~' and 
at a series of temperatures, are given in Tables IV and V. The data 


TABLE IV 
Temperature Dependence of Mechanical Properties for 
the Composite Containing 32.2% Glass Beads 

















Temperature, Ce Eo, E,, 
°C. p.s.i. é. p.s.i. p.s.i. 
—40 220 0.07 3200 660 
—23 240 0.15 1500 360 
-— 7 240 0.24 975 240 
- 1 280 0.30 905 184 
27 240 0.40 580 140 
43 168 0.52 320 115 
60 73 0.32 225 116 
TABLE V 
Temperature Dependence of the Modulus for the Unfilled Binder 
Temperature, E,, 

°C p.s.i. 

—48 2800 

—34 1580 

4 585 

27 386 

49 230 


71 132 
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given in Tables III to V depend on the strain rate, sometimes mark- 
edly. The most marked dependence on strain rate is shown by the 
section of the stress-strain curve for a composite in the yield or de- 
wetting region. Since a complete study was not made, the effect of 
the strain rate and the reproducibility of the tensile data are not ac- 
curately known. 

An analysis of 2) and E, in Table IV and of £, in Table V for the 
unfilled binder shows that these three quantities have the same tem- 
perature dependence. Since these quantities are determined by the 
properties of the binder, this behavior is expected. However, the 
yield stress in Table IV appears independent of temperature within 
experimental error, with the exception of values at 43 and 60° C. 
This temperature insensitivity implies that the filler-binder adhesive 
strength is temperature-independent. A similar conclusion was 
reached, but in a different way, by analysis of the data in Figure 10. 


D. Comparison With Literature Data 


Bryant and Bisset* made an interesting study of vuleanizates which 
contained up to 55% by volume of mineral fillers having particle size 
greater than 1 » (apparently in the range 20 to 40 4). They measured 
stress-strain curves and volume changes that occurred when speci- 
mens were extended in a dilatometer. These investigators also 
studied the effect of per cent filler, particle size (over a rather limited 
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Fig. 11. Replotted data of Bryant and Bisset for a pale crepe vulcanizate 
containing different amounts of whiting. 
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Fig. 12. Variation of yield strain with volume fraction of whiting in pale crepe 
vulcanizates. 


range), degree of cure, and wetting agents for pale crepe and GR-S 
vuleanizates loaded with whiting and other mineral fillers. Quali- 
tatively, their results are similar to those reported here for the glass- 
bead composites, except that they often found an extended plateau 
in the stress-strain curve following the initiation of dewetting. When 
this occurred, the dewetting extended over a larger range of strain 
than that found for the glass-bead composites, although the length of 
the plateau decreased with increasing particle size. 

The dilatometric data of Bryant and Bisset for a pale crepe vul- 
canizate containing between 30 and 45% whiting are shown in Figure 
11, replotted as log V/V» vs. log L/L». This figure shows that the 
Bryant and Bisset data, like those for the glass-bead composites, can 
be characterized by the constant values of Poisson’s ratio listed in 
Figure 11. Yield values «, were also obtained for the vulcanizates 
containing between 20 and 45% whiting by extrapolating their linear 
curves of volume change vs. elongation to zero volume change. The 
yield values are plotted against 1 — 1.105V/” in Figure 12, which 
shows that eq. (13) fits the Bryant and Bisset data much better than 
the data for the glass-bead composites shown in Figure 10. 

Jones and Yiengst® made a dilatometric study of vulcanized rubber 
filled with zine oxide and barytes. The effect of particle size was 
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studied by using five zine oxides having particle diameters between 
0.10 and 0.40 y, and, to a lesser extent, by using barytes having par- 
ticle sizes between 1 and 12 uw. Jones and Yiengst also studied the 
effect of cure time of the rubber. These factors affected the volume— 
elongation curves as might have been expected: i.e., the yield strain 
was decreased by increasing the cure time and the particle size. The 
Jones and Yiengst data, however, seem to show that samples prepared 
with different cure times and with fillers of different particle sizes had 
the same Poisson’s ratio at large extensions, where the dewetting was 
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Fig. 13. Variation of Poisson’s ratio for several filled systems after dewetting. 


essentially complete. If this is correct, Poisson’s ratio, after complete 
dewetting, depends primarily on the volume per cent filler, and not on 
the particle size of the filler or on the modulus or nature of the elas- 
tomer binder. 

To test this hypothesis, values of Poisson’s ratio for several com- 
posites are plotted in Figure 13 against the volume fraction of filler. 
Included are the average values obtained at 25° C. for the glass-bead 
composites and the data of Bryant and Bisset* given in Figure 11. To 
obtain Poisson’s ratio from some data of Jones and Yiengst,’? their 
dilatometric data on vulcanizates containing 0.074, 0.138, and 0.194 
volume fraction of 0.40 » zinc oxide were replotted, and the results are 
also included in Figure 13. This figure suggests that Poisson’s ratio 
is a linear function of the volume fraction of filler and is independent 
of the nature of the binder and filler. Since few data are shown for 
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composites containing filler of different particle sizes, it is premature to 
speculate on the effect of particle size. 


It is a pleasure to acknowledge the assistance of Mr. Dale Jeffries, who con- 
structed the dilatometer and carried out the experimental work. 
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Synopsis 

The mechanical properties of five composites containing up to 60% by volume 
of glass beads (40 to 80 u in diameter) embedded in a polyvinyl chloride—diocty] 
sebacate rubber were studied under large deformations. Elongated ringlike speci- 
mens were extended in a dilatometer at 0, 25, and 50°C., and the volume increase 
was measured as a function of extension. At small extensions, the volume re- 
mained constant; above a critical extension or yield point, however, the volume 
increased because of the formation of vacuoles around the beads. Above the 
yield point, Poisson’s ratio (defined in terms of Hencky strain) was calculated from 
the rate of volume change with extension and was found to be independent of ex- 
tension and temperature, but to decrease linearly with the volume fraction of 
glass. The yield point varied with the temperature and also with the volume 
fraction of beads, as predicted by a theoretical equation. Some stress-strain data 
for the composites are presented and discussed in relation to the dilatometric data 
and the vacuole formation which results from failure of the adhesive bonds be- 
tween the rubber and the glass beads. 
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A Mechanism for Non-Newtonian Flow 


in Suspensions of Rigid Spheres 


IRVIN M. KRIEGER and THOMAS J. DOUGHERTY,* 
Case Institute of Technology, Cleveland, Ohio 


Introduction 


Many colloidal suspensions and solutions of macromolecules exhibit 
a viscosity which decreases with increasing rate of shear. This non- 
Newtonian behavior has been attributed to mechanisms in which the 
shear stress, transmitted through the continuous medium, orients or 
distorts the suspended particles in opposition to the randomizing 
effects of Brownian motion. Variation of viscosity with shear rate is 
then a result of the lowered resistance to flow offered by the oriented or 
distorted arrangement. Kirkwood and Plock’s' explanation of the 
viscosity of suspensions of rigid ellipsoids and F. Bueche’s* theory for 
solutions involve this type of mechanism. 

A more general theory of non-Newtonian flow, based on Eyring’s 
theory of rate processes,’ has been developed by Ree and Eyring.‘ 
In this theory the viscosity is the sum of contributions of an indefinite 
number of unspecified “flow units.’ Since it is not based on any 
particular model, the Ree-Eyring theory should be applicable to all 
systems. Unfortunately, the number of parameters required to de- 
scribe the behavior of a given fluid is not fixed by the theory, nor is it 
possible to deduce the values of these parameters from other known 
properties of the system. Reference to a forerunner of the Ree- 
Eyring theory, in a paper by Tobolsky, Powell, and Eyring,', discloses 
a point of similarity between this theory and the orientation-distor- 
tion theories mentioned earlier. Both involve a competition between 
shearing and thermal forces, and in this competition the size of the 
suspended particle plays a determining role. 


* Firestone Tire and Rubber Fellow. This work forms part of the doctoral dis- 
sertation of T. J. Dougherty. 
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Suspensions of rigid spherical particles are known to exhibit non- 
Newtonian viscosity, but theories involving orientation or distortion 
are based on models which are not appropriate to isotropic rigid 
particles. The Ree-Eyring flow equation has been applied to such a 
system by Maron and Pierce. Using two flow units, one for the 
aqueous medium and the second for the spherical polymer particles, 
these authors were able to fit the viscometric data of Maron and Fok’ 
for a synthetic rubber latex. In their attempts to relate the fitted 
parameters to known properties of the system, however, Maron and 
Pierce achieved only limited success. 

An important clue to the cause of non-Newtonian behavior in rigid 
sphere suspensions is the observation by Maron and his co-workers’~" 
that the behavior of such systems is Newtonian up to relatively high 
concentrations; in none of their latex systems did they detect non- 
Newtonian behavior at concentrations below 20% by volume of sus- 
pended polymer. This fact suggests that crowding plays a vital role 
in the origin of non-Newtonian flow behavior in rigid sphere suspen- 
sions. The theory developed in this paper, therefore, takes explicit 
notice of the interactions between neighboring spherical particles. 
The resultant flow equation is compared with experimental visco- 
metric data on synthetic latexes and solutions of high polymers. 


Derivation of the Flow Equation 


Consider a suspension of rigid spheres of radius a in a continuous 
medium of viscosity »,. When a shear stress 7 is applied, causing the 
suspension to flow with an overall rate of shear 7, the individual 
spheres will tend to rotate with an angular velocity 7/2, dragging with 
them the medium adhering to their surfaces. However, because of 
concentration fluctuations arising from Brownian movement, there 
will be at any instant pairs of spheres whose separation is small, i.e., 
of the order of one diameter. If these nearby spheres were to rotate 
independently, the shear rate in the medium between the spheres 
would be inordinately high, as may be seen in Figure 1A. To avoid 
the consequent high energy dissipation rate, the nearby spheres will 
tend to rotate as a dumbbell about their center of mass until the 
shearing motion has caused them to separate. Figure 1B illustrates 
the hypothetical dumbbell rotation, while Figure 1C is a more realistic 
picture, intermediate between the other two. 
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Fig. 1. (A) Schematic diagram of flow pattern if nearby spheres were to rotate 
independently. (B) Schematic flow pattern if they were to rotate as rigid dumb- 
bell. (C) Diagram of actual flow pattern. 


The behavior described above can be treated quantitatively in the 
following approximate manner. If the separation between two 
spheres permits them to rotate independently of each other, the 
spheres will be called singlets and denoted by P,; if the pair of spheres 
tends to rotate as a dumbbell, the pair will be called a doublet and 
denoted by P:. In the absence of shearing forces, the distribution of 
singlets and doublets can be represented formally as a chemical 
equilibrium : 


ks 
2P, — P, (1) 


Here k, and k, represent the specific rate constants for the formation 
and dissociation of doublets, respectively. Imposition of a shearing 
force introduces a second mechanism for the decomposition of dou- 
blets, since they will be separated during the dumbbell rotation. 
ke 

P, ——— 2P, (2) 
Here k, is the specific rate constant for the shear-induced dissociation. 
Letting m, and n; be the number of singlet and doublet particles 


respectively per unit volume, an overall rate equation can be written 
for the net rate of doublet formation: 


dn2/dt = kmn;? _ (kp + k,)ne (3) 


The condition for a steady state is dn,/dt = 0, and hence 


No ky tr ( EY 
=—!_ =} 4 
n,? ky + Rs ky . ky 
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Since k, will be zero in the absence of shear forces, the zero-shear con- 
centrations of singlets and doublets, denoted by #; and ft respectively, 
will obey the relationship 


fin/n;? = ky ky (5) 


At this point it is necessary to postulate a relationship between th« 
viscosity of the suspension, 7 = r/7, and the concentrations of singlets 
and doublets. We assume that 


n/ns = fim + fore = fin + (fe + 2fi)ne (6) 


Here n is the total concentration of spheres, 


n = nm + 2ne (7) 


The factors f; and f. are assumed to be functions of n, but independent 
of y. At zero shear rate n, will equal #2, and hence the viscosity in 
this limit will be 


no/Ms = fin + (fo — 2fi)Me (8) 


At infinite shear rate, all doublets will be dissociated, and hence n,. 
will be zero. The viscosity 7. at infinite shear rate will then be given 
by 


No/Ns = fin (9) 
Equations (6), (8), and (9) may be combined to give 


9—% fn 
er 00 


and from eqs. (4) and (5), 


Ne n,;* *) tia 

— jpn ie : 

fz ft” ( * ky (11) 
We now assume that only a small fraction of the spheres will be paired 


as doublets, and hence that n; and #; are not greatly different from n. 
With this assumption, eqs. (10) and (11) combine to give 


— Ne k, ani 
14. (4 +) 12) 
10 — Ne ky 
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To evaluate their ratio, we make use of the fact that k, and k, are 
first-order rate constants for the decomposition of a doublet by the 
different mechanisms and hence each is reciprocally related to the 
corresponding mean lifetime. Letting 4% be the mean lifetime of a 
doublet undergoing thermal dissociation and ¢, the mean lifetime of a 
doublet dissociating solely by shearing forces, we obtain 


1 — Ne ty\~" 

~— =(14 ) (13) 

No — Neo ts 
The lifetime of a doublet dissociating under shear is the time required 
for the dumbbell to rotate to a position such that the two spheres are 
in different lamina so that they can separate. Taking the angular 
velocity of the dumbbell as 7/2 and estimating the angle of rotation 
as 1/2, we obtain 
ks | ayn /2 | 


72 |= mr! ly (14) 





ts 


(The absolute value is necessary since ¢, must be positive, and a factor 
a, is included because of the approximation involved in the value of 
the angle.) The mean lifetime under thermal decomposition is 
taken as the time required for the two particles to diffuse apart a 
distance X, 


t = X?/2D (15) 
The diffusion constant D is given by the Stokes-Einstein relation 
D = kT/6xan (16) 


We are assuming here that the diffusion constant of the particle is 
governed by » = 1r/¥, the viscosity of the medium at the prevailing 
rate of shear. Combining eqs. (15) and (16) we obtain 


ty = 3anaX?2/kT (17) 


From dimensional considerations it is evident that the magnitude of 
X, the distance the two particles of a doublet must diffuse apart before 
they can be considered to be singlets, is proportional to the particle 
radius a, and hence 


X? = aa? (18) 
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Thus eq. (17) becomes 
ty = 3ana.u®/kT (19) 


(Einstein’s analysis''! of the distribution under shear of the velocity of 
the medium about a spherical particle leads us to believe that a» 
should be close to unity.) 

The ratio of mean lifetimes is then 


ty/ts = 3apa%|y|/oukT = 3a*|r|/akT = |r| /7, (20) 


where a = a/a, and 


te = akT/3a* (21) 
Equation (13) then becomes 
= | 
1 te (1 . Lz!) (22) 
0 ~ Ne Te 


This relation between y and r is our flow equation. Except for the 
significance attached to r, and the use of the absolute magnitude of 
the shear stress, our flow equation is identical to that proposed on 
empirical grounds by Williamson’? and by Peek and MacLean."* 


Test of the Equation of Flow 


The proposed mechanism was tested by fitting the resultant flow 
eq. (22) to previously published experimental data obtained in this 
laboratory for two synthetic latex suspensions and a polymer solu- 
tion. The data of Maron and Fok described the behavior of a 
butadiene-styrene latex whose average particle diameter was 1390 A. 
as measured by conductimetric soap titration. Maron and Levy- 
Pascal’s data for a Neoprene latex of 1100 A. diameter were also used 
for this test, although the presence of an alginate-type creaming agent 
conferred somewhat anomalous flow behavior on this latex. The 
polymer solution data of Maron, Nakajima, and Krieger'* describe 
the viscosity of a solution of polystyrene of molecular weight 250,000 
in o-dichlorobenzene. For each of the three sets of measurements, the 
experimental temperature was 30°C. 

The following procedure was used in fitting the flow equation to the 
experimental data. By inspection of the data, a trial value of 1, is 
selected to correspond to a viscosity midway between estimated high 
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and low-shear limiting values. A graph of the viscosity 9 vs. r,/ 
(r + 7,) is prepared; if the selected value of r, is correct, the experi- 
mental points should lie on a straight line. If the points indicate a 
curvature, a new value of r, is selected and the process repeated until 
a linear graph is obtained. The intercept of this graph is 7. and 
the slope is 7 — 7... In this manner, the fitted values of r,, m0, and 
n. were obtained for each of the three systems at the various con- 
centrations. 

For each of the two latex systems studied, essentially the same value 
of r, was found at all concentrations. In other words, r, was found to 








TABLE I 
Values of Parameters of Flow Equations 

Volume 
fraction, Tes 

v dynes/cem.? Neo/ Ns no/, 

Data of Maron and Fok 
0.2990 46 3.070 3.532 
0.3970 ” 5.421 7.499 
0.4398 = 7.330 12.689 
0.4876 “ 10.907 24.193 
0.5390 4 18.39 88.80 
0.5603 4 26.92 167 .08 
0.5866 <4 44.67 448.1 
0.6017 < 48.0 1048.0 
Data of Maron and Levy-Pascal 
0.2393 20 2.711 3.151 
0.2821 = 3.561 4.599 
0.3263 - 4.871 7.168 
0.3713 m 7.040 13.690 
0.4147 ™ 11.37 31.39 
0.4429 = 15.65 66.95 
0.4617 ™ 30.60 165.72 
0.5115 . 86.08 1199.4 
Data of Maron, Nakajima, and Krieger 

0.02986 1000 6.205 8.619 
0.04516 1500 9.56 17.54 
0.06627 2000 22.03 35.38 
0.09625 2000 23.85 85.40 
0.12102 4000 57 .64 173.53 
0.14549 5000 52.30 398.3 
0.15948 10000 25.3 921.3 


0.17556 10000 151.7 1673.7 
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be independent of concentration for these systems. For the polymer 
solutions, however, the fitted 7, values increased regularly with in- 
creasing concentration. Table I lists the values obtained for the 
parameters r,,no, and 7. for each system and concentration. Due to 
experimental scatter in the polymer solution data, the accuracy with 
which the parameters could be determined was considerably lower 
than for the latexes. 

To test the accuracy of the flow equation, the viscosities were 
calculated at the experimental shear stresses, using the tabulated 
values of the parameters. Comparison of observed and calculated 
values is shown in Tables II-IV. For the latex suspensions, agree- 
ment between calculated and observed viscosities was within 1% for 
all except the three highest concentrations; even for these concen- 














TABLE II 
Comparison of Calculated and Observed Relative Viscosities. Data of Maron 
and Fok 
T, 
dynes / % % % 
cm.* obs. Neale. CITOF obs. Neale. CFTOF obs. Neale.  efror 
v = 0.2990 v = 0.4876 v = 0.5866 
50 3.293 3.291 0.06 17.18 17.28 0.58 404.10 238.16 41.1 
100 3.218 3.216 0.06 15.17 15.10 0.46 208.78 171.77 17.7 
200 3.157 3.156 0.03 13.44 13.39 0.37 122.21 120.11 1.72 
300 3.132 3.131 0.03 12.70 12.68 0.16 96.58 98.28 1.76 
500 3.110 3.109 0.03 12.03 12.03 0.00 79.28 78.68 0.76 
800 — — — 11.66 11.63 0.26 68.01 66.62 2.04 
v = 0.3970 v = 0.5390 v = 0.6017 
50 6.424 6.417 0.11 52.20 52.12 0.15 1089.3 527.1 51.6 
100 6.081 6.076 0.08 40.02 40.58 1.40 464.0 363.1 21.8 
200 5.802 5.810 0. = 31.55 31.56 0.03 241.8 235.0 2.8 
300 «45.689 5.697 0.1 27.96 27.75 0.75 176.9 180.9 2.3 
500 5.607 5.596 0. > 24.51 24.33 0.73 131.2 132.3 0.84 
800 — — _ 21.98 22.22 1.09 105.3 102.4 2.7 
v = 0.4398 v = 0.5603 
50 9.879 9.897 0.28 107.07 94.07 12.14 
100 8.993 9.019 0.29 71.18 71.08 0.14 
200 8.363 8.332 0.37 52.20 53.13 1.78 
300 8.082 8.042 0.49 45.39 45.55 0.35 
500 7.790 7.782 0.10 39.15 38.74 1.05 
800 7.597 7.621 0.32 34.80 34.54 0.75 
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trated suspensions the agreement is within experimental error except 
for the low shear stress range. The low-shear, high-concentration 
points are in the region where triplet and higher multiplet interactions 
are important, so that the discrepancies in this range are not unex- 
pected. (In fact, these data can be fitted quite well by using smaller 
values of r,, corresponding to an increase in the average dimensions of 
the interacting particles.) For the polymer solutions, the deviations 
between observed and calculated values range up to 5%. In view of 
the larger experimental error in these data, this is considered to be a 
reasonably good fit. 
Since the average radii of the particles is known for the butadiene 

styrene latex, it is possible to compute the factor a in eq. (21), namely 


a = 3a*r,/kT (23) 


- 


a ‘ 
° 
°o 
°o 
———.—¢. Ss WOE 
2.01 0 . a 
ee eee 
ned) 


___¥#0.8603 


2.56+ 


¥20.9390 


v:0.4676 











vr0.4396 





v*0.3970 


0s8i.. —o——0.__ o—— _— 


° 200 400 600 609 





eee v<0.2990 








| 

| 
4 

| 

| 
— 
| 


T, DYNES /cm* 


Fig. 2. Graphical comparison of observed (points) and calculated (solid line) 


viscosities using the experimental data of Maron and Fok. 
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TABLE III 
Comparison of Calculated and Observed Relative Viscosities. 
and Levy-Pascal 





Data of Maron 





T 











dynes / % % %o 

em.? Nobs. Neale. error Nobs. Neale. error Nobs. ‘Neale. eFror 
v = 0.2393 v = 0.3713 v = 0.4617 

30 2.888 2.887 0.03 9.712 9.700 0.12 88.85 84.65 4.7 

50 2.834 2.837 0.10 8.910 8.940 0.33 69.61 69.20 0.58 

100 2.779 2.784 0.18 8.088 8.149 0.75 52.31 53.12 1.54 

150 2.759 2.763 0.14 7.799 7.822 0.29 45.89 46.49 1.30 

200 2.751 2.751 0.00 7.640 7 644 0.05 42.40 42.88 1.13 

300 2.745 2.739 0.22 7.459 7.456 0.04 38.91 39.05 0.35 

500 2.744 2.728 0.58 7.311 7.296 0.20 35.79 35.80 0.03 

800 — — _— -= — — 33.88 33.90 0.06 
v = 0.2821 v = 0.4147 v = 0.5115 

30 3.984 3.976 0.20 19.34 19.38 0.21 _ —_ —_ 

50 3.851 3.858 0.18 17.09 17.09 0.00 481.8 404.2 15.9 

100 3.714 3.734 0.53 14.80 14.71 0.61 278.4 271.7 2.41 

150 3.667 3.683 0.43 13.81 13.72 0.65 214.1 217.0 1.35 

200 3.651 3.655 0.11 13.29 13.19 0.75 186.2 187.3 0.59 

300 3.633 3.626 0.19 12.72 12.62 0.78 156.2 155.7 0.32 

500 3.621 3.601 0.55 12.22 12.14 0.66 131.0 129.0 1.52 

800 —_— _ —_ 11.85 11.86 0.08 113.1 113.2 0.09 
v = 0.3263 v = 0.4429 

30 5.793 5.790 0.05 36.47 36.17 0.82 

50 5.539 5.527 0.21 30.29 30.31 0.07 

100 5.247 5.254 0.13 23.70 24.20 2.10 

150 5.139 5.141 0.04 21.37 21.68 1.44 

200 5.083 5.080 0.06 20.12 20.31 0.94 

300 5.023 5.015 0.16 18.97 18.86 0.58 

500 4.951 4.959 0.16 17.95 17.63 1.78 

800 — — — 17.08 16.90 1.05 





Using 695 A. for a and 46 dynes/cm.? for 7,, we obtain an a-value of 


1.10; as anticipated, this value is close to unity. Using a = 


1.10 in 


conjunction with the r, value of 20 dynes/cm.? from Table I, we cal- 
culate for the Neoprene latex a particle radius of 920 A. as compared 
with Maron and Levy-Pascal’s measured value of 550 A. The prob- 
able source of this discrepancy is the effect of the alginate thickener 
mentioned earlier. 
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TABLE IV 
Comparison of Calculated and Observed Relative Viscosities. Data of Maron 
and Nakajima 




















T, T, 

dynes / % dynes % 

em.? Nmeas. Neale. error em.? meas Neale. error 
v = 0.04516 v = 0.14549 
99 17.03 17.03 0.00 1291 308.7 322.9 4.60 
207 16.55 16.55 0.00 3345 259.0 251.5 2.90 
321 16.17 16.11 0.37 4817 223.7 218.6 2.28 
435 15.80 15.73 0.44 9841 148.2 155.6 5.00 
549 15.32 15.39 0.46 22685 102.5 98.5 3.90 
v = 0.09625 v = 0.17651 
117 85.71 81.95 4.39 6500 1125 1074 4.67 
308 78.81 77.15 2.11 16478 745 726 2.57 
2076 54.82 54.01 1.48 22200 641 624 2.68 
4120 44.59 43.92 1.50 30150 §22 531 1.65 
8315 38.93 39.18 0.64 50050 394 405 2.67 
TABLE V 
Concentration Dependence of Limiting Viscosities 
Low-shear High-shear 
limit limit 

Data of nN no/, p Nao /, Pp 
Maron and Fok 1.119 2.76 0.670 2.58 0.782 
Maron and Levy-Pascal 1.093 3.44 0.541 2.79 0.591 





The increase of r, with concentration exhibited by the polymer solu- 
tions implies that the polymer molecule, regarded as a spherical par- 
ticle, progressively shrinks as concentration increases, changing in 
volume over a 10-fold range. _Maron, Nakajima, and Krieger arrived 
at a similar conclusion from the concentration dependence of 7., 
which they obtained by fitting the experimental data to the Ree- 
Eyring flow equation. The dimension a computed from the extrap- 
olated zero-concentration value of r, is of the order of magnitude of 
the average end-to-end chain length of a freely-jointed polymer chain 
of 2300 styrene monomer units, corresponding to the given molecular 
weight. 
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TABLE VI 
Limiting Viscosities in the Newtonian Range 





Volume 





fraction, 

v (n/n, Jobs Na/ Ns 70/, 

Data of Maron and Fok 
0.0000 1.000 1.000 1.000 
0.0306 1.089 1.093 1.139 
0.0515 1.164 1.166 1.181 
0.1010 1.361 1.369 1.408 
0.1795 1.845 1.821 1.933 
0.2481 2.491 2.420 2.688 
Data of Maron and Levy-Pascal 

0.0000 1.000 1.000 1.000 
0.0091 1.009 1.012 1.035 
0.0391 1.106 1.136 1.166 
0.0769 1.283 1.288 1.370 
0.1160 1.513 1.489 1.642 
0. 1560 1.775 1.755 2.023 
0.1975 2.257 2.118 2.585 





Concentration Dependence of Parameters 


The parameters of the flow equation which were obtained for the 
three systems by the fitting procedure described above vary in a 
regular manner with volume fraction. Because the behavior of r, 
for the polymer solution indicates a variation in the size of the unit of 
flow, it can not be properly considered as a suspension of rigid spheres. 
The constancy of r, for the two latexes, however, has led us to attempt 
to apply functional relationships for the concentration dependence of 
the viscosity of rigid sphere suspensions to the variation of 4 and 7. 
with volume fraction v. 

Those functional relationships which have been proposed and suc- 
cessfully applied over a wide concentration range have two features in 
common. (1) At very low concentrations, they all reduce to the 
linear form 

n/n = 1+ [np (24) 


? 


where [7] is the “intrinsic viscosity;’’ and (2) at some characteristic 
volume fraction p the viscosity becomes infinite. The quantity p 
is called the packing fraction, since the approach to infinite viscosity 
is usually ascribed to the attainment of a close-packed structure. 
Two well-known equations of this form are those of Eilers" 
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is? al or 
/%: = | i ole (25) 

and of Mooney"* 
In(n/ns) = [n]v/(1 — v/p) (26) 











— LOG, (i- Av/g, } 


Fig. 3. Graphical test of eq. (29) for high-shear limiting viscosity data. Straight 
lines are calculated results. Filled circles represent data of Maron and Levy-Pas- 
cal; unfilled circles the data of Maron and Fok. 


These equations can be fitted to our values of 9 and .. over most of 
the concentration range, but an excellent fit over the entire range was 
obtained using an equation derived by a modification of Mooney’s 
functional analysis: 


n/m = (1 — v/p)-?™ (27) 

The derivation of this relationship will be published elsewhere. Be- 

cause of the adsorbed soap monolayer on the latex particles, the vol- 

ume fraction of spherical particles is \v rather than v, where \ has been 
shown to be 

h=1+6A/D, (28) 
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Fig. 4. Graphical test of eq. (29) for low-shear limiting viscosities. Filled cir- 
cles represent the data of Maron and Levy-Pascal; unfilled circles the data of 
Maron and Fok. 


Here A is the thickness of the soap monolayer and D, the average di- 
ameter of the polymer sphere. Equation (27) was modified by sub- 
stituting Av for v, and written in logarithmic form: 


In (n/ns) = — pln] In (1 — do/p) (29) 


Table V lists the values of \ given by Maron et al., and also the 
values of [n] and p obtained by fitting both high and low-shear limit- 
ing viscosities for the two latexes. The goodness of fit can be judged 
from Figures 3 and 4, which is a graph of eq. (29) together with data 
from TablesIand II. The values for [y] are quite close to Einstein’s 
theoretical prediction of 5/2 except for the low-shear value for the 
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Neoprene latex, where the presence of the thickening agent is pre- 
sumably affecting the intrinsic viscosity. 

One further point which requires elaboration is the fact that ob- 
served flow behavior of the latexes was Newtonian at volume frac- 
tions below ca. 0.20, whereas a small but not negligible difference be- 
tween m and 7, is indicated from the fitted equations. This dis- 
crepancy as shown in Table VI is a real one; the Newtonian viscosi- 
ties observed in this concentration range fall between the calculated 
no and »,, and can be represented well by 7. A possible explanation 
is that in this range the counter-ion concentration is not yet large 
enough to screen the coulombic repulsion between particles, which 
arises from the adsorbed soap anions, and hence doublet formation is 
prevented. 


The work reported in this paper was supported in part by a National Science 
Foundation Grant. 
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Synopsis 
A mechanism is postulated to account for non-Newtonian flow in suspensions 


of rigid spheres. On the basis of this mechanism, the flow equation (7 — 7,.)/ 
(nm — no) = (1 + |1r/7e|)~! is derived, where 7 is the viscosity at the shear stress 
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r, » and ».. are the limiting values of » at zero and infinite shear, respectively, 
and +, is a parameter determined by temperature and particle size. This equa- 
tion correlates well with published data on latex and polymer systems. For 
latex suspensions, the parameter +r, does not depend on concentration and for 
these systems the dependence of the relative viscosity at a fixed shear stress, 
nr, on volume fraction of suspended phase, v, is accurately described by the 
equation In », = —pl|n] In (1 — v/p), where [ny] is the intrinsic viscosity and p 
is the volume fraction at close-packing. In general, the parameters [y] and p 
depend on the shear stress. 








TRANSACTIONS OF THE SOCIETY OF RHEOLOGY _ III, 153-160 (1959) 


The Temperature Dependence of Flow 


Birefringence Parameters 


WLADIMIR PHILIPPOFF,* The Franklin Institute Laboratories 
for Research and Development, Philadelphia, Pennsylvania 


Introduction 


In flow birefringence, one question has not received conclusive 
treatment. That is: Is the flow birefringence primarily caused by 
the forces acting in the liquid, or by the orientation of the particles 
determined by the rate of shear? Up to now, the theories of this 
effect have been developed for infinite dilution, where the shearing 
stress and the rate of shear are proportional to each other, so that 
this decision is difficult to make. In preceding papers,'~’ we have 
tried to show that the force equilibrium, in particular the difference 
in principal stresses, Ap, determines the degree of birefringence, 
An, which means the application of the stress-optical law of solids to 
streaming liquids. However, there seems to be a way of proving that 
the stress is primarily responsible for the occurrence of flow bire- 
fringence; that is, by the use of measurements at different tempera- 
tures. With normal liquids that show flow birefringence, the viscos- 
ity can be changed to a considerable degree by varying the tempera- 
ture. Investigations over a range of rates of shear at different 
temperatures can therefore be made. Should the rate of shear, D, 
be the determining quantity, all these measurements should fall on a 
single straight line, when An is plotted versus D. If, however, 
the shearing stress, +, or difference of the principal stresses, Ap, 
is the significant variable for non-Newtonian liquids, plotting bire- 
fringence versus r or Ap should give a straight line, indicating a 
constant stress-optical coefficient, C [expressed in Brewsters (Br.) 
= 10-“cgs units]. The same reasoning applies to solutions. 


* Present address: Esso Research and Engineering Co., Linden, N. J. 
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The first measurements of the influence of temperature on flow 
birefringence of solutions were made by Signer and Gross* on a 7.8% 
solution of a low molecular weight polystyrene in tetralin at 10, 20, 
and 40°C. They found C = An/2r = 5030 Br., independent of 
temperature. Cerf* proposed and performed such measurements on 
solutions in order to decide whether the particles were rigid or soft: 
plotting An against r would give a straight line only for rigid particles. 


Experimental 

We have used substances that are somewhat different from the ones 
usually investigated: (a) a normal National Bureau of Standards 
oil “OB” (polybutene) of about 200 poises viscosity at 25°C., the 
viscosity of which could be changed by a factor of 100 in a compara- 
tively easily accessible temperature range; (b) a medicinal white oil 
(“paraffin oil’’) with a viscosity of about 35 ep. at 25°C.; and (c) 
a solution of 3% polyisobutylene of M ~ 100,000 in the white oil, 
with a viscosity of 129 ep. at 25°C. at low rates of shear. 

The flow birefringence instrument used has not yet been described 
in detail; it is the same one used in our preceding papers.'~7 The 
temperature could be changed between 12°C. (below this temperature 
dew formation was excessive) and 70°C. The viscosity of the 
liquids was measured either with a Cannon-Fenske viscometer or 
with our capillary viscometer. The OB-14 calibration curve of the 
National Bureau of Standards (viscosity vs. temperature for 20—- 
40°C.) was extended to 10 and 70°C., and the viscosity values used 
are interpolated from this curve. 


Results 


In Figure 1, the results of the flow birefringence measurements 
on the OB oil are shown. The extinction angle, x, was 45° in all 
cases, since the OB oil behaved like a normal liquid in this range of 
rates of shear and shearing stresses. The birefringence was positive. 

In Figure 2, results from the measurements on the white oil alone 
and on the solution of polyisobutylene in this solvent have been 
plotted together in a composite plot. The viscosities of the white 
oil and of the solution at the three temperatures are listed in the 
figure caption. The polyisobutylene solution at the higher rates 
of shear was non-Newtonian; its flow curve (viscosity vs. shear 
rate) was measured and the values for the rates of shear used were 
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Fig. 1. Flow birefringence of OB-Oil at 17.8, 24.5, 50.3, and 71.2°C. Viscosities: 
650, 207, 27.0, and 7.15 poises, respectively. 
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Fig. 2. Flow birefringence of white oil and a 3% polyisobutylene solution in 
white oil at 12.0, 25.0, and 60°C. Viscosities of the solvent: 0.7507, 0.3545, and 
0.0885 poises, respectively. Viscosities of the solution: 2.60, 1.29, and 0.295 
poises, respectively. 
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Fig. 4. Plot of An vs. r for the white oil and the 3% polyisobutylene solution in 
this solvent. 


interpolated. In the Newtonian range the measured constant 
viscosity was used to calculate the shearing stress. The birefringence 
was positive for both the white oil and the solution. 


Discussion 
Using the flow birefringence data from Figure 1 and the viscosity 
values, we have plotted An as a function of r in Figure 3. Wesee 


that a complete superposition of all the data is achieved: all the 
points fall on a single straight line. This proves that for the OB 
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Fig. 5. Extinction angle x as a function of the rate of shear D for PIB solution in 
white oil. 
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Fig. 6. Extinction angle x as a function of the shear stress for PIB solution in 
white oil. 


oil, a normal liquid of high viscosity, the birefringence depends 
specifically on the shearing stress. In the case of Newtonian liquids, 
1/,A p = 7, assin 2x = 1. 

The same procedure was followed for the white oil, using the vis- 
cosity and birefringence data from Figure 2. The result is seen in 
Figure 4 (bottom curve). Data on the polyisobutylene solution 
were also taken from Figure 2: '/24 p = r/sin 2x was calculated and 
plotted in Figure 4 (top curve). 

In Figure 5, the extinction angle for the polyisobutylene solution 
as a function of rate of shear is shown. In Figure 6, the same data 
are plotted as a function of shearing stress; in this case all the points 
lie on one single curve. 

We therefore feel that we have given conclusive proof that the two 
birefringence parameters—degree of birefringence, An, and extinction 
angle, x—are exclusively determined by the shearing stress or the dif- 
ference in principal stresses, and depend on the rate of shear only 
through the (possibly non-Newtonian) viscosity. 

The values of the stress-optical coefficient, C, for the four series of 
experiments on the PIB solution are listed in Table I. 
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TABLE I 
Stress-Optical Coefficients C of PIB Solution in White Oil at Different 
Temperatures 








Temp., °C. Cylinder C, Br Error, Br 
12 Il 1435 +73 
25 Il 1400 +75 
25 Ill 1356 +41 
60 Ill 1425 +56 
Mean 1404 +72 (5.5%) 





In Table II the stress-optical coefficients for the different poly- 
isobutylenes and solvents used are listed. There seems to be a 
systematic decrease in C with decreasing molecular weight, probably 
due to the fact that short molecules do not follow the usual coil 
statistics. The OB oil has a value of C roughly 30% lower than the 
rubbery polyisobutylene (B100). 








TABLE II 
Stress-Optical Coefficients and Refractive Indices 
Material Nes C, Br. Remarks 

Polyisobutylene 

B100 1.508 +1500 100% 

B100 1.508 +1500 0.1-9% decalin 

LM/MS 1.507 +1300 100% 

Used 1.473 +1404 3% in white oil 
OB oil 1.492 +975 100% 
Decalin 1.474 +190 — 
White oil 1.473 +550 — 





The case of the polyisobutylene solution just described is remark- 
able in the fact that the white oil solvent has a value of C = 550 
Br.—i.e., about 35% as great as that forthe solute. Notwithstanding 
this large value, the experimental value of C for the solution is 
extremely close to that of a 100% polyisobutylene (Vistanex LM- 
MS) determined previously.’ It is difficult to understand why the 
solvent seems to make no contribution to the birefringence of the 
solution. On the other hand, the index of refraction of decalin’ 
was practically identical to that of the white oil. This would suggest 
the possibility of the C’s being equal also. In Table I, the C’s at 
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different temperatures were listed together with the root-mean- 
square errors. Within the mean experimental error of 5'/2%, we 
could not detect any effect of temperature between 12 and 60°C. 
in a more than 100-fold range of Ap and An. The change in absolute 
temperature should have resulted in a 17% change which is 
about three times the experimertal error. Why this temperature 
independence should exist is not quite clear. Possibly the index of 
refraction. changes with the change in density and diminishes the 
infiuence of the absolute temperature. 

From these measurements, it is also seen that the value of C for 
the white oil, a low molecular weight liquid, is of the same order of 
magnitude as the ones for the polyisobutylene. For that matter, in a 
previous publication’ it was shown that another low molecular 
weight liquid, benzyl benzoate, has a value of about 2730 Br., which 
is much higher than that of the polyisobutylene, and cinnamic alde- 
hyde has a value of 4670 Br. 

The data have also been plotted according to a suggestion of Peter- 
lin; here the quantity plotted on the abscissa is the increase of 
shearing stress by the solute. Practically, we subtracted (neoivent 
X D) from r values. The plot gave a completely straight line up to a 
value of this function of 3000 dynes/cm.*. However, above this 
value the birefringence rose sharply with a deviation from the straight 
line value of up to 40%. This top region for the function is the one 
where the non-Newtonian behavior starts and the influence of x 
on Ap becomes important. We therefore think that the former rela- 
tion, shown in Figure 4, has the more general significance. 

According to the criterion of Cerf,® if the plot of An vs. shearing 
stress gives a single common curve, the particles should be rigid. 
However, this is probably not the case for the polyisobutylene investi- 
gated. It is therefore important to get more data in order to inves- 
tigate this point in greater detail. Also, a continuation of these 
measurements into the range of higher dilutions of the polyisobutylene 
would be important in order to decide whether the C value of 1400 
Br. remains constant as concentration changes. 

Plotting values of the recoverable shear, s, versus shearing stress 
for the PIB solution, we could approximate the curve obtained by a 
straight line in the region of low shearing stress. From the slope of 
this line, a shear modulus, G = 3500 dyne/cm.*? was found. As the 
concentration of the polyisobutylene was 3% by weight and the 
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density of the solution 0.8664 g./cc., we can calculate the molecular 
weight of the polymer, employing a formula used previously:' M 
= RTc/2G, where c is the concentration of the polymer in g./cc. 
Using the G value mentioned, the molecular weight value obtained 
is M = 92,000. This seems to be of the right order of magnitude, 
since the polymer Vistanex LM-S investigated by Leaderman" 
had a molecular weight of 85,400. The value of M for the LM-MS 
polymer measured at 100% concentration and calculated in the same 
way as above is 203,000. These values show that the proposed 
formula is apparently useful for determining weight-average mo- 
lecular weights from viscosity and birefringence data. 


The measurements on the OB oil were performed by Mr. John G. Brodnyan. 
Those on the white oil and polyisobutylene solution were by Gwynn K. Mc- 
Connell, to whom thanks are due. 

The white oil and polyisobutylene solution were supplied by the Esso Research 
and Engineering Company, Linden, New Jersey, who also sponsored the work, 
and whose help is gratefully acknowledged. 
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Synopsis 


Measurements of flow birefringence on two hydrocarbon oils and a solution of 
polyisobutylene in white oil were made at different temperatures. Plotting the 
values as a function of Ap (the difference of principal stresses) gives a single curve 
independent of the temperature. This indicates that the birefringence is pri- 
marily caused by the shearing stresses. 
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Theory and Analysis of Peel Adhesion: 
Mechanisms and Mechanics 


D. H. KAELBLE, Central Research Department, Minnesota Mining 
and Manufacturing Co., St. Paul, Minnesota 


I. Introduction 


Unbonding in peel might be aptly termed the “zipper” mechanism 
for bond failure. A bond involving a flexible member may display 
high strength in shear. As the peel angle increases from zero, the 
force required to initiate and propagate a failure interface decreases 
by orders of magnitude due to a mechanism familiarly termed ‘‘peel’’ 
or “stripback.” 

Bond analysis of overlap shear by Goland and Reissner! identifies 
two sets of bond forces: “tearing’’ forces, which are normal to the 
bond, and “shear” forces, coplanar with the bond. Hata and co- 
workers** recognize and treat the moment arm of the applied force as 
influenced by adherend flexibility and angle ef peel. An internal 
bond stress distribution in peel, for a peel angle of 90°, is calculated 
by Spies,‘ who also provides an explicit. solution for peel force as- 
suming cohesive failure of a Hookean interlayer. Bikerman® has 
more recently treated peel phenomena at a peel angle of 90° with 
explicit solutions for both Newtonian interlayer flow and cohesive 
failure of a Hookean interlayer. 

An objective of this discussion is to provide: (a) a definition of peel 
which applies where either one or both of the adherends are flexible, 
and (b) a definition general with respect to angle of peel, and modes 
or mechanisms of bond failure. 


Il. Theory 
A. Definition of “Simple Siripback”’ 


The descriptive assumptions for an idealized type of peel termed 
“simple stripback”’ are listed: 
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1. A steady-state unbonding process. 
. Properties of the flexible member: 
a. of slender rectangular cross section 
b. elastically (Hookean) deformed 
c. bonded and unbonded portions of semi-infinite length. 
8. Properties of the forces: 
a. applied force—acts in tension along the central plane of the 
fiexible member. 
b. shearing forces—distributed over the bonded interface of 
the flexible member. 
c. cleavage forces—a highly localized parallel array, perpen- 
dicular to the bond plane. 
d. all forces uniformly distributed across the bond width. 


%© 


The first assumption defines the unbonding as progressing at some 
constant velocity, other than zero, and is described by a balance of 
applied and reactive forces and moments of force. The second as- 
sumption describes a thin ribbon adherend, elastically deformed in 
tension and flexure, whose deformation is uncomplicated by end 
effects. The properties of the forces as indicated in the third assump- 
tion will be subsequently amplified upon. 


B. Deformation Curve of the Flexible M ember 


Preliminary to attacking the broad definition of peel, let us examine 
the characteristics of the elastic bending curve of the flexible member 
under tension. This deformation curve controls two important 
factors in peel: (a) the effective moment arm of the applied force, and 
(b) the curvature at the bond boundary and the resulting stress dis- 
tributions. The following discussion represents a special case of 
buckling of centrally compressed bars as presented by Timoshenko.® 
For simplicity it is assumed that the bonded portions of the flexible 
member satisfy the conditions of a built-in end and that weight and 
centrifugal effects may be neglected. 

Figure 1 presents a schematic profile of this deformation curve. 
The exact differential expression for elastic bending is: 


EI(d6/ds) = pm’ (1) 


By differentiating eq. (1), introducing the relation dm’ = ds sin 6, and 
multiplying by d@/ds we obtain: 
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EI(d6/ds) d(dé/ds) = P sin 6 dé 
or 


EI §,° (dé/ds) d(dé/ds) = PJ,” sin 6 dé 


fe 











Fig. 1. Deformation curve of the flexible member. 


By integrating and taking into account boundary conditions, namely 


that curvature and angular deflection are zero when m’ = 0, we 
obtain: 

EI [{dw\? 

3 (*) = P(1 — cos w) (2) 


which represents an independent differential expression for the de- 
formation curve. By writing eq. (1) for the case where m’ = m and 
6 = w, substituting into eq. (2), and rearranging terms, the following 
relation results: 

2EI (1 — cos w)]}” 
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Rewriting eq. (2) we obtain: 


; 2P (1 — cos >| , 
d ‘dis = . 
'w/ds | EI (4) 
where eq. (4) expresses (dw/ds), the boundary curvature, and eq. (3) 
expresses m, the moment arm, in terms of its applied force P, the 
flexible member stiffness E/, and the peel angle w. These equations 


differ from the usual circular bending approximations by the appear- 
ance of the factor 2 in both expressions. 


C. The Steady State in Terms of Forces and Moments 


A schematic cross-section of the bond boundary is shown in Figure 
2. Two coordinate systems are presented with the (2,yo) system 
representing the stationary substrate and the (z,y) coordinates 
defining the center of action and the translating force system. The 
line of action of P, defined by its w and m, lies outside the flexible 
member in the section of this view. 
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Fig. 2. External forces acting on the flexible member (schematic). 


As indicated in the assumptions, two systems of bond forces or 
potential mechanisms of failure in peel are identified. These are 
cleavage and shear. A rigorous definition of these stress distribution 
functions is not attempted here since this would involve further 
restrictive assumptions regarding the bond interlayer and failure 
interface. Evaluation of the requirements for steady-state equi- 
librium produce the following relations: 

for the forces 


P= (P. + P,)'” (5) 
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P. = Psinw = Sf, — Ff. (5a) 
P, = Pcosw = Xq (5b) 


for the moments of force 


M = Pm = M.+ M, (6) 
M,. = Pm = Zfiz (6a) 
M, = Pm, = hq (6b) 
and for the moment arms 
2EI (1 — cos w) }” 
m = | P | (7) 
m=m—h _ (7a) 
m, =h (7b) 


These reactive forces and moments are simply defined as total sum- 
mations of their positional increments. However, without sacrificing 
generality some description is possible. 

The cleavage forces f are highly localized at the bond boundary 
where the flexible member undergoes a sudden change from maximum 
bending to zero curvature. At high stripping angles regions AO and 
OB are, respectively, zones of tension and compression. The work 
of Spies‘ adequately justifies the implied assumption that the cleavage 
forces become negligible after the initial region of boundary tension 
and interior compression. The general form for the cleavage stress 
distribution function resembles that of a highly damped harmonic 
decay curve. 

The shear stresses may be uniformly distributed throughout the 
bond for inextensible flexible members. An assumption of finite 
extensibility results is some form of an exponential decay function 
with maximum stress at the boundary region. Except at very low 
angles of peel it may be shown, after the manner of Spies, that the 
maximum shear stress is negligible with respect to the maximum 
tensional cleavage stress. 

In “simple stripback”’ as defined here, w may vary from 0 to 180°. 
The reactive bond forces depend upon the magnitudes of both the 
sin w and cos w components of the applied force. At a peel angle of 
90°, neglecting bending strains, only cleavage forces exist. At some 
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low angle of peel the summed cleavage forces are negligible compared 
to shear forces. In addition, due to the offset distance h of P and 
shear forces 2q, at some low angle prescribed by eq. (7) and eq. (7a) 
the cleavage forces become compressional in region OB and tensional 
in regior VA. 

This discussion is confined to those conditions of applied force and 
peel angle where the flexible member half-thickness fh is negligible 
compared to the magnitude of m. Under these conditions the cleav- 
age failure mechanism predominates and the following relations 
defining this mechanism apply: 

by restatement of eq. (5a) 


P = (2f, — =f.)/sin w (8) 
by substitution of eq. (6a) and eq. (7) where m, = m: 
P = (Zfa,)*/(2EI(1 — cos w)] (9) 


Equation (8) states that a direct proportionality exists between the 
net summation (Zf, — =f.) of bond cleavage forces and P. Equation 
(9) expresses a second-power relation between P and the total sum- 
mation of cleavage moments of force (Zf,,). 

The sensitivity of the peel force to bond strength variables is 
readily understood in terms of eq. (9) due to the exponential relation 
between Zfz, and P. The distribution of cleavage bond forces is 
shown to be a function of peel angle by eq. (8) since at w = 90°, 
=f. — =f. = Pandatw = 180°, =f, — =f, = 0. 


D. Peel of Two Flexible Members by Cleavage 


The treatment of the mechanics of peel tests may be extended to 
include the system involving two flexible members joined by an 
adhesive bond or interlayer. The conditions relating to the “simple 
stripback” of a flexible member from a rigid substrate apply. In 
addition, it is assumed that the moment arm of forces P, and Ps, are 
large compared to the combined thickness of the flexible members 
and the bonding interlayer. The geometry of this two component 
system is presented schematically in Figure 3. 

From the geometry and conditions of steady-state equilibrium it 
may be observed that 
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(a) for the system: 


at+B=f (10) 
Pa = —Ps (11) 
(2fa)da = —(Zarde (12) 
(6) for the individual members A and B: 
(pm), = (Zfaida (13) 
(pm)s = (2fird)e (14) 


(fic), (fit), 
¥ 

















(fid, Aft), 


+ a 


Fig. 3. Peel of two bonded flexible members. 


oe * a > 


6 


From eqs. (10) through (14) it may be seen that: 
Mm, = Mz =m (15) 
From the deformation curve of an elastic member bent in tension: 


e (Zone) a (0 = sos)" 
= P = P 


By the introduction of eq. (11) and accounting for the fact that 
cos 8 = —cos a, eq. (16) may be written as: 








(16) 


A B 


(El) 2/(EIT), = (1 — cos a)/(1 + cos a) = (tan a/2)? (17) 
If Jz is known, then Eg may be calculated by rewriting eq. (17) as: 
Ex, = ((EI)4/Ixs] (tan a/2)? (18) 


If one of the members has known values of E, and J,, and if P, and a 
are measured, then the summation of moments in the bond may be 
calculated by the following expression : 
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(2fitde = —(2(EIP), (1 — cos a))” (19) 


Equations (18) and (19) suggest that any tape construction may be 
bonded to an ideally elastic substrate, and both the summation of 
tensional and compressional moments of force in the bond (Zf,) and 
the net effective modulus of elasticity of the tape in flexure may be 
measured simultaneously. 


E. Work of Peel 


Utilizing the theory of peel, the bond strength in peel may be 
expressed in energy or work terms after the fashion of Rivlin.’ The 








ANN he LNNNS 
Fig. 4. Work of peel. 


effective peel work may be defined as the product of P and a displace- 
ment AL in the direction of this force, as indicated by the following 
expression : 

W = PAL (20) 


A schematic sketch describing this relation is shown in Figure 4. 
For “simple stripback,” where the shape of the deformation curve is 
unchanged, the geometry indicates: 


L = (1 — cos w) AX (21) 


where AX is the equilibrium displacement of the bond forces. By 
substitution of eq. (21) into eq. (20): 


W = P(1 — cos w)AX (22) 


On comparing eqs. (22) and (9) it may be seen that the common term 
(1 — cos w) occurs. In eq. (22) this term is related to the degree of 
translation of bond forces relative to that of the applied force; in eq. 
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(9) it relates the magnitude of P tom. By combining eqs. (22) and 
(9) we obtain: 


W = (2f,a))*AX/2EI (23) 


Equation (23) expresses the ‘“‘effective work of peel’’ in terms of the 
summation of moments Zf,v,; independent of the geometry of the 
test (i.e., the stripping angle). If =(f,,) is independent of stripping 
angle, the work will be independent of stripping angle. 


Ill. Experimental 


Experiments designed to provide support for the theory were con- 
ducted on bonds containing a pressure-sensitive adhesive interlayer. 
The particular advantage of this system is due to the fact that an 
equilibrium bond is easily established and under constant conditions 
of rate and temperature the peel force is essentially invariant, thus 
fulfilling the condition of steady-state unbonding.’ In addition, since 
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Fig. 5. Photographic study of tape unwind (90° angle of peel). 
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the force of peel, for pressure-sensitive adhesive bonds, changes 
markedly with both rate and temperature, the apparent bond strength 
may be systematically varied by controlling these physical conditions 
of test.* 

Unless otherwise noted, the rigid adherend consists of a rigidly 
clamped cellophane surface. Bonds were established by multiple 
passes of a mechanically driven '/,” diameter steel roller loaded to 3.0 
lb. Tests were conducted at 25 + 1°C. and 50 + 5% relative 
humidity. The combined reproducibility of the bonding process and 
peel test was about +4%. 


A. Photographic Study of Roll Unwind 


The special apparatus developed by Stedry"' produced data applied 
to the test of eq. (9). Figure 5 provides a schematic view of the ex- 
perimental setup. Highspeed motion picture films were taken of the 
stationary stripping point for various rates of peel. The force P and 
the deformation curve of the flexible member were then measured. 
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Fig. 6. Schematic sketch of stripback testing device. 
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B. All-Angle Tape Stripback Testing 


A special apparatus was designed to obtain precise peel force 
measurements over wide ranges of controlled rate and angle of peel. 
This tester was designed as an accessory to the Instron Tensile Tester'* 
which provides a wide range of stripback rates (0.02 to 20 in./min.) 
and an extremely sensitive load sensing and recording system. Sup- 
plementary equipment included a temperature control oven per- 
mitting tape bonding and stripback testing over a temperature range 
from —50 to 200°C. A single test strip of 12” length provided data 
for rates of 0.02, 0.2, and 2.0 in./min. over the full range of angles. 
An additional 12” strip supplied the 20 in./min. results as well as 
additional rechecks of the data obtained for the slower rates. A tape 
can be tested in this manner in 1 or 2 hours, depending on the number 
of points measured. 

The principle of the stripback operation is illustrated diagram- 
matically in Figure 6. After the tape has been bonded to the surface 
of the stripping wheel and the angle of stripback w selected, the ten- 
sion line (one end of which is attached to the fixed top frame of the 
Instron) drives the stripping wheel at a perimeter velocity just equal 
to that of the moving crosshead. The bonded test tape advances to 
the stripping point B at the same rate that the unbonded tape de- 
parts, thereby establishing a condition of transient equilibrium which 
maintains the sighting angle ¢ and the stripping angle w constant. As 
long as the stripped tape and the tension line connecting the trans- 
lating and stationary portions of the apparatus are parallel and do 
not change length, this condition is maintained independent of the 
rate of stripback. To change the angle of stripping, the operator 
changes the lengths of either the tape or the tension line, and B moves 
to a new equilibrium position between points A and C. 

The reference idler is located along a line of centers between the 
load cell and the crosshead mount and at the intersection of two 
perpendicular planes tangent to the stripping wheel at points A and C. 
An inspection of the geometry of this system indicates the following 
equation : 

1— sing 
1 — cos @ 
relating the sighting angle ¢ and the angle of stripback w. Utilizing 
a specially clesigned protractor sighting device provided with parallax 


w (degrees) = 90 — @ + tan 
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correction, the sighting angle may be determined within +0.2° and 
the stripping angle from a plot of ¢ vs. w. 


C. Measurement of Young’s Modulus 


Young’s modulus values for the flexible members were obtained 
from tensile stress-strain curves measured on the Instron Tensile 
Tester. 


IV. Results and Discussion 
A. A Test of Theory 


A good test of the basic assumptions involved in the theory of tensile 
peel of a flexible from a rigid member is gained by applying eq. (3) 
to real systems. The data required for this test were provided by an 
earlier photographic study of the geometry of peel-back of cellophane 
tape conducted by Stedry.'' Table I presents information concerning 
the test conditions, nature of the bond, and measured values of P and 
flexible member m. m is assumed equal to the measured average 
radius of curvature of the delaminated backing. 








TABLE I 
Stripping Force and Backing Moment Arm for 90° Cellophane Tape Roll Unwind 
Rate, mxX PmxX Pm®x EX 
Temp., R. H. em./ P, 10%, 103, 105 10-5, 
°F. % Adhesive min. Ib. in. Ib. in. Ib. in.* Ib./in.* 
70-80 34-41 Rubber-resin 1 0.138 32 4.16 13.3 4.58 
blend 10 0.53 16 8.50 13.6 4.70 
100 1.28 12 15.4 18.5 6.40 
1000 «1.98 9 17.8 16.0 5.53 
74-76 24-26 Acrylate 1 0.53 22.8 12.1 27 .6 9.53 
copolymer 10 0.75 19.7 14.8 29.2 10.1 
100 1.03 15.4 15.9 24.5 8.45 
1000 1.45 10.1 14.7 14.9 5.16 
73-74 11-16 Acrylate 1 0.84 17.6 14.8 26.0 8.96 
copolymer 10 1.15 144.6 16.8 24.5 8.45 
1000 «61.50 11.3 17.0 19.2 6.63 





Equation (3) may be rewritten in the following forms: 


pm* = 2EI(1 — cos w) (24) 








ane oe 


4 
@ 


24) 
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E = pm*/(2](1 — eos w)] (25) 


Equation (24) predicts the quantity (pm?) to be a constant character- 
istic of the backing, for a fixed angle w. The calculated values of 
(pm*) tabulated in Table I for each set of stripping data substantially 
confirm this predicted constancy over a wide range of stripping rates 
and forces. The moment of inertia J of the flexible member about its 
neutral axis may be calculated by applying the following expression: 


I = b(2h)*/12 = 1.45 X 10°” in. 


where for the cellophane tape, b = 0.05” and 2h = 0.0015”. The 
values of Young’s modulus of elasticity EZ based on stripping meas- 
urements calculated from eq. (25) are tabulated in Table I. The 
tangent modulus values reported in Table II, resulting from room 
temperature (25°C.) measurements at low extensions (1.25% or less), 
are in reasonable qualitative agreement with those calculated from 
stripping data. 


TABLE II 
Young’s Modulus of Cello #600 Put 78 Tape Backings 








Test Temp., °C.: 75 25 —18 
R. H., %: . 50 _ 
% Elongation E xX 10°, E xX 10°, E X 10-4, 

Ib/in.? lb/in.? Ilb/in.? 

0.25 6.70 8.53 15.1 
0.75 5.20 6.71 11.9 
1.25 3.80 3.65 8.20 
1.75 2.50 1.85 4.15 
2.25 1.95 1.31 2.24 
2.75 1.43 0.99 1.41 
3.26 1.14 1.33 1.43 
3.75 _— 0.95 1.15 
4.25 0.90 0.64 





B. The Force-Angle Function 


If the summation of moments Zf,z; of tensional and compressional 
forces in the bond are independent of the stripping angle and cleavage 
is the controlling failure mechanism, then for steady-state stripping, 
from eq. (9), P should vary inversely as (1 — cos w). 








174 


D. H. KAELBLE 


















o° 
OD este = 200 1 /mn 
. 0 
* 
° 
. 9 
= « > a 
g 10 Co 
s s 
& 
~ % 
* $ 
ae 
y 
a 
iro 
{2} 
4 
a 
s x 
a re) 
oor 
ec 
f 1-08 
00s 4 5 n —_— 
° “a so no 1@0 200 240 


STRIPPING ANGLE ( DEGREES) 





Fig. 7. Peel force P vs. angle w for polyester tape. 


Several types of peel test results confirm this P—w relation. 


ping angle for three pressure-sensitive tape bonds. 
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Figures 
7, 8, and 9 present steady-state peel force data as a function of strip- 
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Fig. 8. Peel force P vs. angle w for cellophane tape. 
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Fig. 9. Peel force P vs. angle w for spring steel tape. 


varied for the three bonds is the composition and flexibility of the 
flexible member, as described in Table III. These data are typical of 
pressure-sensitive type bonds where the distribution of forces in the 
bond is strongly influenced by rate effects. Plotting stripping force 
P on a logarithmic ordinate scale facilitates comparison of force-angle 
plots at different rates. Changing the stripping rate and the re- 
sultant Zf,2, values should merely translate the force-angle curve 
along the log P ordinate without changing the shape of the function, 
if the inverse relation between stripping force P and (1 — cos w) 
applies. A change in the shape of this force-angle plot, at constant 
rate, indicates the deviation of experimental data from theory. 

To illustrate the effect of stripping rate upon the shape of the force- 
angle plots, eq. (9) is represented graphically in Figures 7 through 9, 
where 


Lfia,/2EI = K 


At low stripping rates, and accompanying lower stripping forces, 
experimental data confirm extremely well the theoretical inverse P vs. 
(1 — cos w) relationship over a wide range of angles. At higher 
stripping rates the deviation becomes more marked. It may also be 
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TABLE III 
Comparative Properties of Tape Bonds at a Stripping Angle of 90° 
(Adherend: cellophane. Adhesive interlayer: 
rubber-resin; thickness, 0.001 in.) 





Flexible member type 














Polyester Cellophane Steel 

b, in. 0.5 0.5 0.5 

2h, in. 1.0 x 107° 1.5 xX 10-3 1.0 x< 10-3 

I, in.* 4.16 X 107"! 1.45 k 10-” 4.16 x 10°" 

E, \b./in.* 5.0 x 105 6.8 xX 105 1.0 x 10’ 

ET, lb. in.* 2.08 X 10-5 9.85 < 10-5 4.16 x 10~* 

Peel Properties 

Rate, 

in./ P, Zhi» P, Ziti, P, Ziti, 
min. Ib. Ib. in Ib Ib. in lb. Ib. in. 
0.02 0.044 1.35 10-* 0.13 5.1 107? 0.32 16.3 x 10-3 
0.20 0.079 1.80x10-* 0.31 7.8x< 107? 0.50 20.8 x 107% 
2.00 0.39 4.00 x 10-* 0.81 12.6 10-* 1.00 28.8 x 10-' 
20.00 1.22 7.10 x 107 — 1.909 39.8 x 10-3 





noted that the deviation is uniformly more marked as the stripping 
angle is increased for all rates. A qualitative explanation of these 
deviations is available when the dual effect of tensile and flexural 
strains upon the Young’s modulus of the backing in the region of the 
stripping is considered. On the convex side of the flexed backing the 
tensile and flexural strains are additive, and at higher stripping forces 
non-Hookean, with an accompanying lower Young’s modulus. As the 
stripping angle increases, the backing length under high flexural strain 
increases. Both these effects tend to cause non-elastic buckling of 
the backing, thus decreasing the backing m, and requiring a corre- 
sponding increase in the stripping force P to maintain the conditions 
of static equilibrium required by eqs. (6) and (6a). The non-Hookean 
type flexural deformation of the convex side of the flexible member is 
not fully recoverable, and produces set or curl which increases with 
rate and angle of delamination. It may be noted that the force-angle 
functions for the more elastic spring-steel tape (see Fig. 9) more closely 
fit the theoretical curve at higher levels of P, in agreement with the 
above argument. 

Hata has given particular attention to the relation of stripping 
angle to stripping load for thin films of benzoyl cellulose, celluloid, 
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Fig. 10. Peel force P vs. angle w for benzyl cellulose from Dural (fixed angle 
method). 




















] 
' 
3 
3 iS} 
| 
& 
a 
+ OF 
rw] 
3 
oS>- 
c=) 
+ 4 
<= ° ' 2 5 Ss 5 - 
1/1- cos Ww 
Papen Type Samece No Symsou 
Uncoateo Boox 28 —x— 
Coates Boox 5376 —o— 
Orrser 100 —3— 


Fig. 11. Peel force P vs. angle w for paper bond strength (Bekk Surface Tensile 
Tester Method). 


and polyvinyl butyral bonded directly to metal or glass by solvent 
activation or combined heat and pressure.** Results of this study 
indicate good agreement between experimental results and the pre- 
dicted linear inverse relation between stripping force P and (1 — 
cos w) for bonds involving no adhesive interlayer. Plots of these data 
for the tensile peel of benzyl cellulose from Dural are presented in 
Figure 10. 

The Bekk Surface Tensile Tester measures the force necessary to 
tear the surface layer of a strip of paper which has been bonded to a 
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Fig. 12. Effect of the physical properties of the flexible member. 


flat metal bar. The effect of angle of force application upon tearing 
force has been studied by the Institute of Paper Chemistry.'* Overa 
range of angles from 36 to 90° the linear inverse relation between 
tearing force and (1 — cos w) show good agreement with reported 
experimental data for this type of peel, which involves intemal failure 
of the flexible member. Plots of peel force (P) versus [1/(1 — cos 
w)]|, in Figure 11, give linear curves for internal peel delamination of 
three types of paper over the above-mentioned range of stripping 
angles. 

Thus, this definition of peel would appear to provide a satisfactory 
explanation of the dependence of stripping force on angle for (a) p.s. 
bonds, (b) bonds having no highly elastic interlayer, and (c) delamina- 
tion involving cohesive failure. It is interesting to note that while 
the peel force, in general, increases with decreasing angle of peel, the 
work of peel as defined by eqs. (22) or (23) remains constant or de- 
creases with peel angle. 


C. Influence of ihe Flexible Member 


Additional information concerning the distribution of force in 
the bond is supplied by peel tests of a pressure-sensitive bond in which 
only the thickness and Young’s modulus of the flexible member are 
varied to produce corresponding variations in EJ. Descriptions of 
the bonds and values of P at rates of 0.02, 0.20, 2.00, and 20.0 in./min. 
(w = 90°) are summarized in Table III. The apparent summation 
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of moments of the adhesive bond Zf,, is calculated by eq. (9). The 
dependence of Zfz; upon EI is illustrated graphically in Figure 12. 
At the lowest rate of stripping, 0.02 in./min., where the flexible mem- 
ber undergoes predominantly Hookean deformation, Zf,2; would 
appear to be directly proportional to EJ. At high rates of stripping, 
due to the buckling effect of the backing, eq. (9) no longer accurately 
applies and calculated values of Zf; are smaller than the true value, 
producing the upward curvature in the curves represented in Figure 
12. These results might indicate that for simple stripback, at con- 
stant rate, (dw/ds) is a constant and Zfq; is influenced only by the 
magnitude of EJ. The Bikerman® equation for 90° peel involving a 
Newtonian adhesive interlayer may be written: 


K nb 


/2 
_ ~ '/s 
ste (12 tm ") (#1) (26) 


and for cohesive failure of a Hookean interlayer: 


0.716 o. 
P = —T7, (EN (27) 








where b is the bond width, a the interlayer thickness, K/t, a peel-rate 
factor, 7 the interlayer Newtonian viscosity, oo the adhesive stress at 
failure, and Y the interlayer Young’s modulus. For 90° peel, eq. (9) 
becomes 


P = (Zfia:)*/2EI (28) 


Combining eq. (28) with eq. (26), we find, other factors being held 
constant, that f,z, is proportional to (EI)"* for a Newtonian inter- 
layer. A similar combination of eq. (28) and eq. (27) predicts that 
Sf; is proportional to (EI) for cohesive failure of a Hookean 
interlayer. 


Conclusions 


In conclusion, it is suggested that this definition of peel may have 
the following applications: 

1. Provide a basis for a more general theoretical analysis of peel 
adhesion, involving explicitly-defined interlayer rheology and mech- 
anism of failure. 

2. Assist in the better understanding of present diverse methods 
of bond-strength testing of adhesive joints involving flexible members. 
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Synopsis 

An idealized type of tensile peel is defined which substantially represents real 
systems in the range of low peel forces. The theory applies to the steady-state 
unbonding of (a) a flexible and a rigid member and (6) two flexible members. 
Failure interface propagation is discussed in terms of two stress mechanisms. The 
first is boundary cleavage, which involves compression-tension bond stress, and 
the second boundary shear, involving bond stress in simple shear. The theory is 
general with respect to the stripping angle, and a definition of work of peel is de- 
veloped which is independent of this variable. Experimental confirmation of the 
theory is provided by photographic studies of the peel process. A new experimen- 
tal device is described which permits precise measurement of peel force over a wide 
range of controlled angles of peel, rate and temperature. Bonds involving both 
soft viscoelastic interlayers and no interlayer are evaluated and discussed in terms 
of the theory. 
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The Behavior of Jets of Viscoelastic Fluids 


FREDERICK H. GASKINS* and WLADIMIR 
PHILIPPOFF?, Laboratories for Research and Development, 
The Franklin Institute, Philadelphia, Pennsylvania 


Introduction 


Many investigations have been conducted into the behavior of 
liquid jets.'~* These have dealt principally with free-flying jets of 
simple viscous liquids, for which the primary parameters have been 
described; they are the viscosity, density, and surface tension. 
When polymer solutions or other viscoelastic liquids are being dealt 
with, the existing theoretical treatments have proven to be inadequate 
to account for the observed behavior. The investigations‘ that have 
been conducted with polymer solutions have generally emphasized the 
importance of viscosity, but have made little reference to the effect of 
the elasticity necessarily introduced. 

As part of a general experimental program into the factors con- 
tributing to the stability of liquid jets, we have had to deal with visco- 
elastic behavior. 

This paper presents experiments that were made with viscous and 
viscoelastic liquids in jets in which the marked qualitative difference 
in velocity behavior is readily observable. The experimental con- 
ditions under which these tests were made are described, particularly 
because the constant-volume driven gun developed for this application 
is uniquely advantageous to define accurately the initial velocity im- 
parted to extrude the test fluid. 

The Hagenbach treatment® of the behavior of jets of simple viscous 
liquids expelled after Newtonian flow through a capillary (nozzle) has 
been modified and extended to be applicable when non-Newtonian 
flow is present. Further, a treatment is derived for the behavior of 


* Present address: Aeroprojects Inc., West Chester, Pa. 
t Present address: Esso Research & Engineering Co., Linden, N. J. 
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viscoelastic liquids which is analogous to the Hagenbach treatment for 
viscous liquids. This novel treatment leads to an evaluation of the 
elastic energy present in the jet by analysis of the kinetic energy 
balance occurring in terms of the velocity relationships. 


2. Experimental Conditions 
2.1. Instrumentation 


The fluid jets are extruded at a known velocity by expelling a 
volume of the test fluid through a cylindrical nozzle at a predeter- 
mined rate. The instrument devised for this purpose is based on a 
design by Gill and Gavis,’ which has been modified to give a constant 
efflux volume; it is shown in Figure 1. It consists of a cylinder and a 





Fig. 1. Jet extrusion instrument. 


piston which are rigidly mounted on a massive framework. The 
piston, which transmits the expulsion force to the fluid, is actuated by 
a so-called bridge-jack which is driven by a reversible motor through a 
series of reduction gears of adjustable ratio. A ‘delayed reversal” 
switch has been incorporated to allow the automatic reversing of the 
direction of piston travel. This prevents the piston from overrunning 
and breaking the front endplate. 

The pressure exerted inside the cylinder can be determined by a 
gage which is mounted on the front endplate of the instrument. 

Cylindrical nozzles of various diameters are used with this instru- 
ment in order that a wide range in velocities may be available. Table 
I shows the dimensions of the nozzles and the theoretical velocity 
range that can be impressed. The rates of shear listed are those oc- 
curring within the nozzles. It can be seen that the theoretical veloc- 
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TABLE I 
Dimensions and Operating Ranges for Jet Extrusion Nozzles* 

N h Di Theoretical impressed Rate of shear 

Nos- Lengt sneneter velocity, v, cm./see. D, sec.—! 

zle L, d, 

no. cm. cm. Minimum Maximum Minimum Maximum 
1 5.97 0.1819 1.29 x 10' 3.26 x 10? 5.68 x 10? 1.43 x 105 
2 5.74 0.0391 2.80 x 10? 7.07 x 10 5.74 x 10* 1.45 x 10° 
3 5.98 0.152 1.84 < 10' 4.65 x 10* 9.68 x 10% 2.45 x 108 
4 2.49 0.1827 1.28 x 10' 3.24 x 10 5.61 x 10? 1.42 x 105 
5 2.54 0.0473 1.92 X 10? 4.85 x 10‘ 3.24 x 10‘ 8.19 x 108 
6 2.48 0.1210 2.92 x 10! 7.37 x 10? 1.93 x 10% 4.87 x 105 
7 2.53 0.0193 1.15 x 10* 2.91 X 10° 4.77 x 10° 1.21 x 10 
8 2.53 0.0612 1.15 x 10% 2.89 x 10 1.51 XK 10* 3.79 x 10° 
* Minimums based on gear ratio, n = 0.0629. Maximums based on gear 


ratio, n = 15.89. 


ity range is 12.8 to 291,000 cm./sec. The practical upper limit is dic- 
tated by the mechanical setup, and would certainly fall far below the 
theoretical value stated. Nevertheless, an extensive range in operat- 
ing conditions is available. 

The jets are extruded horizontally, the trajectory being measured 
directly and/or photographed over the entire path of flight or any 
segment thereof. The photography is carried out with a Graflex 
camera and either a Strobolume flash unit of 9 usec. duration or a 
modified flash unit of 2 usec. duration. 


2.2. Experimental Parameters 


By having the cylindrical capillary nozzles in series with the main 
cylinder, the conditions under which the flow occurs can be exactly 
pre-determined. The primary flow parameters are the impressed or 
“continuity” velocity, 


vo = Q/Az = (A;/A:2)(dL/dt) in em./sec. (2.1) 


where, v» is the initial mean relative velocity of jet to air, in cm./sec.; 
Q is the volumetric rate of extrusion, in cm.*/sec.; A; and A, are the 
cross-section areas of the main cylinder and the nozzle, respectively, 
in em.?; (dL /dt) is the rate of piston travel, in cm./sec.; and the 
velocity in flight evaluated using the ballistics equation which assumes 
a parabolic trajectory, 
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vp = 2V g/2y, inem./sec. | (2.2) 


Cir 


where, vy is the horizontal component of the flight velocity of the jet 
relative to the air, integrated over the range, in cm./sec.; zx is the 
horizontal distance traversed (range), in cm.; y is the distance of fall, 
in em.; g is the gravitational acceleration constant (980.7 cm./sec.’) 

In addition, two rheological parameters can be readily determined. 
One is the rate of shear occurring in the nozzle, 


D = 4Q/axr* = (4/r)vo, in sec.~ (2.3) 


where r is the radius of the nozzle, in cm., and vp is the “continuity 
velocity.”” The other is the shearing stress, 


r = rP/2L, in dyne/cm.? (2.4) 


where P is the pressure drop across the nozzle (gage reading), in 
dyne/cm.’, and L is the length of the capillary nozzle, in cm. 

Note: The value of r calculated by eq. (2.4) has to be corrected for 
kinetic energy, etc., to obtain the “true” shear stress developed within 
the nozzle; by the same token, the rate of shear according to eq. (2.3) 
must be corrected by (n + 3)/4 when non-Newtonian flow occurs, if 
n = d log D/d log r is significantly larger than unity. These evalua- 
tions are discussed in detail in reference 8. 


2.3 Test Fluids 


Dibutyl phthalate (DPB), water, and an unmodified lubricating 
oil were the pure (viscous) liquids used during this investigation. The 
polymers, all dissolved in DBP at various concentrations, were: 


Polyisobutyl methacrylate (PIBMA) supplied by the Chemical 
Warfare Labs., U. 8. Army Chemical Center. 

Cellulose acetate—butyrate copolymer (CAB), commercial grade 
(EAB-38 1-20). 

Cellulose nitrate (NC, no. 1) commercial grade (RS'/, sec.). 

Cellulose nitrate (NC, no. 2) dynamite grade. 


The physical properties of the solutions were measured. For 
0.5%-by-weight solutions and for DBP the densities were between 
1.04 and 1.05 g./em.’, the surface tensions were approximately 33 
dyn./cm., and the relative viscosities were 1.2 to 2.4 based on a vis- 
cosity of 0.1598 poise for DBP at 25°C. 
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The rheological properties of some of the solutions were determined 
in a high pressure capillary viscometer using the “capillary length 
effect.’* The results are shown in Table II. 


3. Experimental Results 


Preliminary tests with viscoelastic solutions showed that the flight 
velocities were often smaller by a considerable amount than those im- 
pressed at the nozzle. This deviation from the theoretical expecta- 
tion, which, in a first approximation, would be expected to hold, 
resulted in a careful analysis of the experimental setup. By having 
the values of vp calculated from the piston travel rate [the second part 
of eq. (2.1)], direct measurements were made of the efflux volume 
against time. This provided a check on the first part of eq. (2.1), 
and gave results which were precise to better than 1%. Therefore, 
the instrumentation was performing properly. 

The test fluids were then examined systematically in order to define 
better the relationship between the flight velocity, vy, and the “con- 
tinuity velocity,” vo. The materials were investigated at a series of 
velocities in the range of 100 to 1500 cm./sec. using one or more cylin- 
drical nozzles with results as follows. 


3.1. Viscous Liquids 


The water was expelled at 766 cm./sec. from a nozzle with a diam- 
eter of 0.1827 cm. It was found that the water sprays apart immedi- 
ately upon emergence from the nozzle so that no evaluation of v, 
could be made within the velocity range cited. 

The lubricating oil (7 = 19.2 eps. at 37.5°C.) at the same m gave 
a vy of 1050 cm./sec.; the ratio vo/vp was 0.730. The values of v, 
in this case are unreliable because the jets formed sprays which cannot 
be evaluated accurately in this way. 

The dibutyl phthalate was measured in five capillary nozzles 
(nos. 1, 4, 5, 6, and 8) with a four-to-one difference between the largest 
and smallest diameters. The results as shown in Figure 2 are that all 
of the data fit a straight line with a mean slope of 0.801 within experi- 
mental error. These data were evaluated by the method of least 
squares to get an accurate result. The divergence above 10 m./sec. 
is due to the jet’s disintegrating into a spray; below this level, uni- 
form jets with a parabolic trajectory were obtained. It should be 
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Fig. 2. Velocity relationships for dibutyl phthalate. 





noted that one of the runs was made with a nozzle of 5 cm. length in- 
serted into the cylinder. This test of the “Borda effect,” to be dis- 
cussed, is indistinguishable from the other runs with the same cylin- 
drical nozzles in the normal or extended position. 


3.2. Viscoelastic Liquids 


(a) Polyisobutyl Methacrylate in Dibutyl Phthalate. The PIBMA 
in DPB was investigated at 0.5 and 0.25% concentrations (by weight 
of solution), the results for the two concentrations being qualitatively 
the same. For the 0.5% solution four capillaries were used (nos. 4, 
5, 6, and 8). The results in Figure 3 show a marked difference from 
the results obtained for viscous liquids. None of the curves, for the 
nozzles used, fit either the theoretical ratio of (vo/vy) = 0.801 or the 
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Fig. 3. Velocity relationships for 0.5% polyisobutyl methacrylate in dibutyl 
phthalate. 


1:1 relationship. Rather there is a significant decrease in vy for a 
given vp. Also, the magnitude of the difference increases with de- 
creasing size of radius. (Note: The smaller the radius, the higher 
the rate of shear range involved; therefore, while the curves are cor- 
relatable on the basis of equivalent jet velocities, they represent dif- 
ferent initial flow conditions within the nozzles.) 

An interesting point to be observed is that the curves are generally 
linear and are displaced parallel to the theoretical curve, which repre- 
sents the DBP data. 

(b) Cellulose Derivatives in Dibutyl Phthalate. The cellulose- 
derivative solutions were investigated in nozzle no. 4 (d = 0.1827 
em). The results for 0.5% cellulose acetate—butyrate, 0.5, 1.0, and 
2.0% cellulose nitrate (commercial grade) and 0.5% cellulose nitrate 
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Fig. 44Velocity relationships for various polymers in dibutyl phthalate. 


(dynamite grade) in DBP are shown in Figure 4, together with the 
data for 0.5% PIBMA in DBP using the above-mentioned nozzle. 

As seen in this figure, the solutions of the two commercial grade 
polymers, CAB and NC #1, gave results which coincide entirely with 
the results for DBP up to the velocity at which disintegration of the 
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jet occurs. On the other hand, the highly nitrated (high molecular 
weight) cellulose in solution shows a drop in vy as a function of v that 
is even greater than that obtained for the PIBMA. Here, again, it is 
apparent that the displacement of the lines is parallel to that of the 
theoretical, despite the magnitude of the velocity difference. 


4. Discussion of Results 


The investigation of the velocities of free-flying jets gave results 
which were contrary to expectation on two counts: the Newtonian 
liquids showed a 20% increase in flight velocity compared to the im- 
pressed, and the non-Newtonian materials showed a marked decrease 
in many instances. A systematic examination of this behavior gave 
the quantitative results presented, which show that this is not a 
random occurrence. Therefore, the explanation had to lie in the 
theoretical treatment applied to the particular conditions of operation. 
The primary parameter of interest is the discharge (or contraction) 
coefficient, C, = vo/vy, so far as the theoretical flight velocity to be 
expected at a given continuity velocity level is concerned. Addi- 
tionally, the exact conditions, which result in the large velocity dif- 
ferences shown by the viscoelastic solutions, have to be analyzed and 
evaluated quantitatively. 


4.1 Discharge Coefficient for Newtonian Liquids 


The observed 20% increase of the flight velocity over the impressed 
or “continuity” velocity did not agree with the theories of viscometry, 
which in the Boussinesq’ treatment lead to a value of m = 1.12 
(equivalent to C,). This would mean a flight velocity lower than m. 
Analysis of the conditions show that the Boussinesq treatment applies 
in the case of a regular viscometer in which the active hydrostatic 
pressure head is diminished by kinetic energy to an available pressure 
head. Specifically, the parameters are calculated in this case from 
the total pressure utilized, knowing a priori that there are losses other 
than that to viscous flow. 

The calculations of Hagenbach* in his original publication, in which 
he treated the boundary conditions of our problem directly, lead to a 
discharge coefficient of ~0.8. Hagenbach’s reasoning assumed an 
exit velocity that would not be changed by the kinetic energy in the 
liquid, even though the kinetic energy is very high. This is the case in 
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our jet extrusion instrument, where the speed of the motor is not in- 
fluenced by the pressure in the liquid; also, the parameters are evalu- 
ated from the change in kinetic energy at the nozzle exit and in flight. 
Therefore, Hagenbach’s theory is the pertinent one for these condi- 
tions. 

Examination was also made of the classical statement that the dis- 
charge coefficient is caused by the liquid entering a capillary tube, 
with the flow outside of the capillary leading to the so-called “vena 
contracta.”” We checked the results of this theory, as reported by 
Birkhoff and Zarantonello,"® by reversing the cylindrical nozzle so 
that it was inserted into the cylinder with a 360° influx. This is 
called the “Borda tube effect.” The data reported in the reference 
gave a change in C, from 0.611 for the extended capillary to 0.500 for 
the inserted, presumably using water as the test medium. Our results 
show no detectable difference in the velocities whatsoever using 
capillaries of d = 0.182 cm., with one L = 5.97 em. extended and 
inserted, and another L = 2.40 cm. extended. These results show 
that the Borda tube effect may be damped out by the viscous dibutyl 
phthalate, whereas the effect is reputedly quite noticeable in the exper- 
iments described. Also, the results with three geometrical configura- 
tions establish that laminar flow is present, and no turbulence occurs 
within the tube. These aspects have been substantiated by evalua- 
tion of the Reynolds numbers involved which shows that critical 
values are not attained within this velocity range (to 1500 cm./sec.). 

The treatment of Hagenbach is repeated in the Appendix; however, 
it is given in a different way from the original in order that the result- 
ing equations can be extended to the case of non-Newtonian flow. 
The result of the evaluation is that 


Cy = vo/vr = 1/¥/2 = 1/1.2599 = 0.7937 (4.1) 


This theoretical value is equal, within experimental error, to the 
value of 0.801 determined experimentally. This means that for a vis- 
cous, Newtonian liquid, the diameter of the jet decreases by VC, = 
0.890 or by 12% as compared to the diameter of the capillary tube. 
Also, the flight velocity to impressed velocity relationship is 


Vp = vo/C, = W/2u) = 1.260 (4.2) 











192 F. H. GASKINS AND W. PHILIPPOFF 


4.2 Discharge Coefficient for Non-Newtonian Liquids 


The Hagenbach derivation in Appendix I is applicable to the behav- 
ior of Newtonian (viscous) liquids. For viscoelastic polymer solu- 
tions, the flow behavior is generally non-Newtonian. Therefore, the 
discharge coefficient must include the ““Power Law,’ which is likewise 
included in the evaluation of velocity gradient for the capillary 
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Fig. 5. Evaluation of discharge coefficient as a function of the “viscous” slope. 


viscometer flow of these materials. The Hagenbach treatment has 
been extended to include the “Power Law,” but not possible elastic 
effects, i.e., 


ov/Or = Kr" (4.3) 
This derivation, which is presented in Appendix II, leads to 
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where 
n = Olog D/O log r (4.5) 


is the ‘‘viscous slope,’’ i.e., the slope of the D or (0v/dr) vs. r curve in 
double-logarithmic plot (‘flow curve’). 

The C,, so evaluated, is the discharge coefficient which correlates 
the continuity velocity impressed on the viscoelastic fluid with the 
theoretical velocity of the jet in flight assuming non-Newtonian vis- 
cous flow of the fluid through the nozzle. The numerical evaluation 
of eq. (4.4) has been carried out, the results being presented in graph- 
ical form in Figure 5 for n between 1 and 10. It should be noted 
that this correction for the behavior of non-Newtonian fluids in free- 
flying jets is equivalent to the standard correction for the parabolic 
distribution of velocity across the radius in capillary viscometry 
(cf. ref. 8). For n — ©, a constant velocity over the entire cross 
section results and C, = 1. Non-Newtonian flow, alone, therefore, 
can change C, only from 0.7937 to 1.0000. 


4.3 Elastic Energy in Jets of Viscoelastic Fluids 


The significant decrease in flight velocity as compared to the im- 
pressed velocity, for viscoelastic fluids, could not be accounted for 
by the existing theories, which deal only with viscous behavior. 
Nor does taking into account the non-Newtonian flow of these solu- 
tions explain the large difference. As seen in Figure 5, C, varies 
between 0.79 and 0.94 under ordinary conditions; this means that 
the theoretical flight velocity for any material after undergoing 
Newtonian or non-Newtonian viscous flow should be from 6% to 26% 
higher than the impressed velocity. Therefore, the elastic energy 
developed in viscoelastic materials under flow conditions has to be 
considered. In reference 8 it was stated that “the elastic energy, 
imparted to the liquid at the entrance of the capillary, is carried out 
of the capillary by the flowing liquid and dissipated outside the capil- 
lary much in the same way as the kinetic energy.”” Equations were 
derived therein for evaluating the elastic energy as recoverable shear, 
8, or normal stress, P,, from the pressure vs. L/R relationships ob- 
tained during the flow through the tube. Under the present ex- 
perimental conditions, we have a direct manifestation of the elastic 
energy from the effects obtained outside of the tube; these effects are 
that the jet diameter increases and the flight velocity decreases. 
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Nevertheless, the continuity equations are satisfied in that volume- 
time constancy is obtained. 

The total energy transported out of the capillary tube is recovered 
as both kinetic energy and potential (or elastic) energy. Knowing 
the kinetic energy imparted to the fluid (from the theoretical flight 
velocity, vy = v/C,) and the kinetic energy utilized to propel the 
fluid through space (from the actual flight velocity, vr), we can eval- 
uate the elastic energy as the difference in these kinetic energies. 
In order to do this, one has to evaluate an “elastic discharge coeffi- 
cient” to take into account any nonlinearity in the shear stress, r, 
vs. recoverable shear, s, bchavior, since the shear modulus, G = 
7/8, is generally not a constant. 

The derivation of the elastic parameters is given in Appendix III, 
but can be summarized as follows: 

The elastic discharge coefficient is calculated from the slope of the 
shear stress-recoverable shear curve, assuming the same velocity 
distribution over the cross section as is usual in the treatment of non- 
Newtonian flow, as: 


C, = (n’ + 1)(n’ + 3)/2 (4.6) 
where 
n’ = Olog 7/0 log s (4.7) 


is the elastic slope, i.e., the slope of the 7 vs. s curve in double-logarith- 
mic plot. 

The numerical evaluation of C, shows that as n’ varies from 0 to 1, 
C, takes on values from 1.5 to 4.0; it is thus possible for four-fold 
differences in velocities to occur that account for the observations 
noted. 

Having the elastic discharge coefficient, the elastic energy can be 
evaluated as normal stress, 


Py = C.C,(vr? — ve*)p (4.8) 


where, Py, is the normal stress at the wall, in dyne/cm.*; C, is the 
fluid density, in g./em.*; vr is the theoretical flight velocity = 
vo/C,, in em./sec.; and vy is the actual flight velocity, in cm./sec. 
This eq. (4.8) probably does not represent the final formulation. 
For example, one factor that should be included is the stated effect 
of the elasticity on the flow pattern in the nozzle; in other words, the 
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change in shear stress and rate of shear distribution across the radius 
under the influence of normal stress must be calculated. 

The elastic energy in viscoelastic fluid jets has not yet been quan- 
titatively evaluated by the method described; additional data from 
the capillary viscometer and from flow birefringence in the same 
range of D are required before experimental proof of the relationships 
can be offered. 
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Fig. 6. Recoverable shear vs. shear stress for three different polymers in dibutyl 
phthalate (at 25°C.). 


Some data were obtained in the capillary viscometer using the 
“eapillary length effect’’ technique (see ref. 8). These results, which 
are giver. in Figure 6 as a plot of log s vs. log r, show the comparative 
elastic nature of three of the solutions. The commercial grade nitro- 
cellulose solution shows an s of unity (100% shear) this being attained 
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at a rate of shear (>10° sec.—') which is beyond that used in the jet 
extrusion experiments. In the range of the jet operation (up to 
50,000 sec.—'), no elastic effects were detectable in the capillary, and, 
as was shown in Figure 4, no velocity decrease was observed. The 
PIBMA solution shows a large s even in the range of jet operation, 
and this result correlates qualitatively with the large difference in the 
vy vs. v9 curve. The CAB solution has a large s at high rates of shear, 
but the values are too small to be measurable in the range of D below 
50,000 sec.~'; this again correlates directly with the jet velocity com- 
parison, which shows no difference for this solution. 

Examination shows that elastic energy is dissipated within milli- 
seconds after emergence of the jet from the nozzle. In the cases 
where the jet was analyzed, the relaxation times were always below 30 
millisec., usually much smaller. After this relaxation period, the 
jet velocity remains constant throughout the remainder of flight 
(unless of course, the jet disintegrates). The constancy of the 
velocity has been carefully checked by microscopic measurement of 
the jet diameter along the entire path of flight and measurement of 
the true parabolic shape of the trajectory in several cases, and the 
stated results held true. 


5. Summary 


The behavior of fluids in free-flying jets has been quantitatively 
described in terms of the flow parameters involved during jet ex- 
trusion. The treatment of this problem has been extended to take 
into account not only the Newtonian flow, but also the non-Newtonian 
response. Particularly significant is the development of a treatment 
analogous to the above by means of which the elastic energy developed 
within the fluid during its flow through the extrusion nozzle can be 
evaluated solely in terms of the kinetic energy relationships outside 
of the nozzle. It must be restated that this equation for evaluating 
elastic energy as normal stress probably does not represent the final 
formulation; additional experimental confirmation has to be obtained. 
Nevertheless, this is a new method for measuring elastic energy, and 
one which appears to have advantages for the investigation of visco- 
elastic materials due to the simplicity of the experimental technique 
involved. 
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APPENDIX I 


Modification of Hagenbach Derivation of the Discharge 
Coefficient for Newtonian Liquids 


The kinetic energy of the mass of emerging liquid with a certain 
velocity distribution over the radius is equated to the kinetic energy 
in flight, assuming a uniform velocity, v, over the area. The mass 
of the liquid is given as the area multiplied by the velocity and den- 
sity. The kinetic energy is determined by the mass multiplied by 
1/w*. For the capillary tube, the velocity must be integrated over 
the radius to obtain the basic equation: 


i? '/,-2arprdr-v-v? = '/srR?- pvp (A-1) 
This can be simplified to: 
2S," v*rdr = Rv,;* (A-2) 


The velocity distribution over the capillary radius is obtained by 
integrating Newton’s Law. The result is: 


: V = Umax[1 — (r/R)?] (A-3) 
Introducing eq. (A-3) into eq. (A-2), we obtain: 


2mae! fo [1 — (r/R) = op! (A-4) 
Let y = 1 — (r/R)?*, then dy = —2(r/R)d(r/R) 
and 
—Vmex' Sy y'dy = —Vmex?-y*/4|} = 1/amex? = v7" 
Umax = WA-vy (A-5) 


We now apply the continuity equation: 


Q = rR? - v9 
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where Q is the efflux volume, in cm.*/sec. and 
Q= Qn Sf," vrdr = Setmahs” [1 — (r/R)? |rdr 
= 2etmex!/y?|) = (#/2)R max 
Equating the two expressions for Q: 
Q = rR? -v9 = (4/2) R Vmax 
Vo = '/2 Umax 


Introducing eq. (A-5) we obtain: 


Vo = '/eW/4-vp 
or, for the velocity of the jet in flight: 
vp = (2/4) = 2-0, (A-6) 
The discharge coefficient C, for Newtonian flow is 
C, = vo/vp = 1/*/2 = 1/1.2599 = 0.7937 (A-7) 
APPENDIX Il 


Extension of the Hagenbach Treatment to 
Non-Newtonian Liquids 


The first expression to be considered is the relationship between 
the velocity integrated over the radius, r, and the assumed uniform 
velocity, v, to the velocity gradient: 


v—% = = (=) dr (A-8) 


The velocity gradient can be expressed in terms of the shear stress by 
introducing the ‘Power Law” (see ref. 8): 


(5) aE? of () = Krp" (7) (A-9) 


Incorporating eq. (A-9) into eq. (A-8), we obtain: 


r r\" dr 
y—» -R{'kK (Z) 
0 , TR R R 
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Let x = r/R. Then: 


. 1 : g?t 
v— UW = Krat-R [ z"dz = Kr," R (. + ‘) 


For +=0,v=%; for z=1,v=0 and —» = Kra'h 
n+ 1 
Therefore: 
v— Wo = —max"t! 
or: 
v = v9 (1 — 2**!) (A-10) 


Introducing eq. (A-10) into the Newton’s law integration (A-4) 
results in: 


2 Vmax’ Se (1 a ade xdx : (A-11) 
The expansion is: 


2 Omen’ Se (1 — 3a"t! + 32°"+? — x3"+8) xdzx 


a a 
“13 #3 ° +4 S845) |’ 


Be 
Vmax = V 6 6 2 A-12 
PF Ra 2 4 ( ) 


n+ n+4 3n+5 





For 
n—~1l, VY +~wW 


n—> ©, VY >! 


Applying the continuity equation, including the expression for v (in 
eq. A-10) gives: 


Q = 2a f,* urdr = 2eR? J.” vo (1 — xt) xdzx 


and 


1 1 
= 2 1 — "+2 = 2 Pe ae = 
Q = 2eR? vax Sc (2 — 2"t*) dx = 2Rmax [5 P | (A-13) 
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Forn—~1,[ |—'/s Therefore: 


l 1 
R19 = 2eR mex _-— 
ate _ E n+ | 


l l 2 
» = 2 max _ = Umax A- 
” : E | | n 2 | \A-t6) 


Combining eq. (A-12) with eq. (A-14) gives 


: ; 9 
vp = — A-15) 
. ZV er sa on +5 \ 
n +3 
For 
3 | as 
n-— 1, (the Newtonian case), vy = 2 .” = ~/2-% 
n> ©, 0p = % 


The discharge coefficient is 











C.=- = ———— ———————— (A-16) 
ee ae: apenas HT: 
n+3 2n + 4 3n + 5 
where 
a E log | 
~ Ld log r 
APPENDIX III 
Evaluation of Elastic Energy in Viscoelastic Fluid Jets 


The energy balance existing for viscoelastic flow through a capillary 
tube was stated in ref. (8), eq. (11), as: 


PrQ = PQ + mpv*Q + PQ (A-17) 


where P; is the total pressure applied, in dynes/cm.*; P, is the pressure 
taken up by viscous flow, in dynes/cm.’; P, is the pressure transported 
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out of the capillary as elastic energy (P?,-Q), in dynes/em.*; Q is the 
efflux volume, in cm.*/sec.; and mpv*Q is the kinetic energy, in dyne- 
em. (per sec.), (m is substituted for by C, as in the Hagenbach treat- 
ment). 

At the exit of the capillary (nozzle), the viscous energy (PQ) is 
zero. Then the equation, reduced and restated in terms of velocity, 
is: 


C,pv°Q = Cypvr*Q + PQ (A-18) 


Equation (A-18) states that the elastic energy is the difference be- 
tween the kinetic energy imparted to the fluid and that utilized to 
propel the fluid jet through space. 

The mean value for elastic energy over the capillary cross section is 


P,- eR? = f," 2ardr J; rds (A-19) 


rf @f~ am 


We approximate the general nonlinearity between +r and s by a 
Power law: 


or 








r = As" 
where 
as E log 4 
0 log s 
Then, from eq. (A-20): 
ft rds = £ As"ds = — £ sg” +1 (A-21) 
0 0 n+l 


The following two relationships established in our previous publica- 
tions:*!! 


Py = 7-s yields Py = As” + 
and 


Pu = (Pu)n (Z) 
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are incorporated into eq. (A-21) to give: 


Pu -< (Pude (zy * 
n’ +1 n’+1\R 
Therefore: 


Pe eis ee ee aif nai 
hha on +i (7) (7) % Pade| + 1) (n’ =| 








(A-22) 
Introducing eq. (A-22) into eq. (A-18) results in: 
2 
= aD cus 2 
(Pu)e E +1) n+ 5 | C.p(v vp’) 
or 
(Pure = CL (v0? — ve’) p (A-23) 
where 
‘+ 1) ‘+3 
c.<.% + 1) (n’ + 3) (A-24) 


2 


is the “elastic’’ discharge coefficient. 

In order to complete the evaluation of normal stress (i.e., elastic 
energy) we have to refer to the impressed velocity. This means that 
the contraction of the jet which would occur for viscoelastic fluids 
even in the absence of elastic phenomena, has to be taken into account 
as: 


vr = v/C, 


where vr is the theoretical flight velocity. 
The fina] equation for the normal stress is: 


(Pu)e = CL Avr* — vp*)p (A-25) 
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Synopsis 


The stability of liquid jets is well known to be described by the theory of 
Rayleigh which uses the density and surface tension of the liquid as important 
parameters. Further investigations have taken into account the viscosity, the 
combination of these three parameters being known as the Weber number. 
However, polymer solutions are now known to have an elasticity in flow; the 
mechanism of breakup of solutions of polymers has not heretofore been investi- 
gated. We have used a constant volume-driven extrusion device by means of 
which the velocity of the emerging jet could be calculated exactly, which gave us 
an additional parameter as compared with the standard pressure-driven device. 
A number of pure liquids and polymer solutions were investigated. The general 
result was that the “vena contracta’’ causes an increase of the velocity of the 
free flying jet by about 20% for pure liquids, as compared to the predetermined 
“continuity velocity.”” This can be explained by the old theory of Hagenbach. 
Viscoelastic polymer solutions, however, showed a decrease of velocity as com- 
pared either to the normal liquids or to the “continuity velocity.’ Very prob- 
ably that difference is caused by the elastic energy accumulated in the jet which 
absorbs some of the kinetic energy of the emerging liquid. The quantitative 
evaluation showed that we can determine the elastic stress of the liquid and 
therefore also the normal stress by measuring the difference between the jet 
velocity and the velocity of the solvent. 
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Communication Problems in Industrial Rheology 


Committee Recommendations for Describing Flow Phenomena 


Increased study of the flow characteristics of polymers and dis- 
persions has revealed a lack of authoritative definitions and stand- 
ardized terminology for stress-dependent flow phenomena. Con- 
fusion resulting from these deficiencies has become so acute that a 
Committee on Communication Problems in Rheology has been or- 
ganized to establish areas of agreement in describing the flow charac- 
teristics of materials. 

The Committee is sponsored by several American Chemical Society 
divisions (Agricultural and Food, Colloid, Industrial and Engineer- 
ing, Paint, Petroleum, and Polymer), the Society of Rheology, and 
the American Society for Testing Materials. 

Before delving into more controversial details, including the defi- 
nition of viscosity, the Committee considered that it must first es- 
tablish and define those broad aspects of flow concerning which no 
serious disagreement exists among practicing rheologists. It was 
decided that the following realities must be recognized as the frame- 
work upon which definitions and terms can be built: 

1. The basic measured quantities are shear stress (S), shear rate 
(D), and time. 

2. Shear-dependent properties are necessarily time-dependent, 
and confusion results from attempts to measure both phenomena 
simultaneously. 

3. The terms used to describe flow behavior are the main source of 
confusion; the concepts are fairly well established. 

The recommendation is made that flow phenomena be described 
using a broad classification of behavior based on the dependence of 
applied stress on shear rate. Three broad areas of flow behavior are 
recognized. They are based on the ratio, S/D. A constant S/D 
ratio is characteristic of a Newtonian system; a decreasing ratio with 
increased D is shear thinning, and the material it characterizes is a 
shear-thinning system; an increasing ratio with increased D is shear 
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thickening, and the material it characterizes is a shear-thickening sys- 
tem. 

The Committee recommends further that flow phenomena be ex- 
pressed in terms of the directly measured variables S and D, and that 
reference to terms such as viscosity and fluidity either be avoided or 
qualified whenever the possibility exists for their misinterpretation in 
a quantitative sense. It is essential that the dependence of S on D 
be determined under steady-state conditions. 

The usage of new words was considered by the Committee which 
sampled the opinions of academic and industrial personnel engaged in 
making and describing measurements of flow behavior. Sentiment 
was strongly in favor of reducing the number of terms used to describe 
the three main types of flow behavior, rather than adding two ad- 
ditional terms. 

The scientific community is strongly urged to adopt the terms shear 
thinning and shear thickening in future communications regarding 
shear-sensitive flow behavior. 


Raymond R. Myers, Chairman 
Donald R. Brookfield 
Frederick R. Eirich 

John D. Ferry 

Leonard Stoloff 

Ralph N. Traxler 
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Abstracts of Additional Papers Presented at the 
Society of Rheology Meeting, 
November 5-7, 1958 


Ultrasonics and Its Use in Rheology 


T. A. LITOVITZ 
Catholic University of America, Washington, D.C., and 
Nawal Ordnance Laboratory, White Oak, Maryland 


Ultrasonic spectroscopy has demonstrated the existence of volume 
viscosity in most liquids and has shown that there are many different 
molecular mechanisms responsible. In all cases in liquids, the ab- 
sorption of sonic waves has been related to relaxation processes and 
not to resonance effects. In this paper, resonance and relaxation are 
discussed in relation to the molecular basis of the ultrasonic loss (or 
volume viscosity) of liquids. 


[To be published in J. Acoust. Soc. Am.] 


Viscosity of n-Hexadecane 


A. B. BESTUL and R. C. HARDY 
National Bureau of Standards, 
Washington, D.C. 


Attempts to prepare high-purity n-hexadecane from commercial 
cetane by simple laboratory procedures were unsuccessful. Frac- 
tional distillation at reduced pressure of material previously treated 
with silica gel produced a few small fractions of about 99.5 mole-% 
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purity. The kinematic viscosities of the original material and of the 
fraction of highest purity, 99.63 mole-%, were found to be the same, 
4.4635 cs. at 20°. All other fractions had lower viscosities. The 
viscosity of NBS standard sample of n-hexadecane, 99.94 mole-%, 
was found to be 4.4642 es. or 3.4540 ep. at 20°C. n-Hexadecane is 
not likely to be suitable for use as a second calibration standard for 
viscometry until a simple, easily defined routine of purification is de- 


veloped. 


Mechanical Resonance Dispersion in Crystalline Solids 


EDWIN R. FITZGERALD 
The Pennsylvania State University, University Park, 
Pennsylvania 


Measurements of the complex shear compliance at closely spaced 
intervals in the range of 100 to 5000 cps. have resulted in the discovery 
of sharp resonance dispersions in polycrystalline metals, crystalline 
polymers, and single crystals of various materiais such as sodium 
chloride, copper sulfate, and quartz. The results indicate that the 
phenomena of low frequency resonance dispersion in mechanical com- 
pliance are of general occurrence. The data can be analyzed and 
fitted closely by means of a generalized stress-strain relation involving 
a linear combination of the strain and its first and second time deriva- 
tives. Inclusion of the third derivative leads to an anti-dissipative 
term which can be used to account qualitatively for a negative ab- 
sorption or induced emission observed over a narrow frequency region 
in some materials. The negative absorption seems to be connected 
with static clamping stresses applied to the sample, decreases slowly 
with time, and can be reactivated by restressing. 

An explanation of the resonances is suggested by calculations of 
Fermi, Pasta, and Ulam for nonlinear systems in which no tendency 
toward equipartition of energy among modes wasfound. Accordingly 
it is proposed that: (1) in crystalline solids, because of the nonlinear 
nature of the binding forces, equipartition of energy among the various 
modes of vibration does not take place; instead energy actually is 
passed back and forth among relatively few of the available modes, 
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(2) the frequency of the energy exchange, in contrast to the lattice 
frequencies themselves, may be very low, e.g., in the audio-frequency 
range, and (3) the observed resonance dispersions occur at frequencies 
corresponding to these various acoustic exchange frequencies. This 
lack of energy equipartition will not necessarily be noticed in specific 
heat measurements of solids, but may have some consequences in 
other areas. 


Viscous Motion of Dislocations in Crystals 


J. J. GILMAN and W. G. JOHNSTON 
General Electric Research Laboratory, 
Schenectady, New York 


Direct observations of dislocations in LiF crystals show that there 
is a quasi-viscous resistance to dislocation motion. This resistance 
does not involve dislocation interactions, but occurs for dislocations 
moving in dislocation-free regions of the crystals. The viscosity in- 
creases with impurity content, with neutron bombardment, and with 
decreasing temperature. Measurements have been made over a wide 
range of dislocation velocities extending from a few atom distances per 
second up to 10" atom distances per second. 

Our microscopic measurements of dislocation behavior have been 
applied to the calculation of macroscopic stress-strain curves with 
good results. 


[To be published in J. Appl. Phys.] 


Radiofrequency Spectroscopy Applications 


WILLIAM L. ROLLWITZ 
Southwest Research Institute, San Antonio, Texas 


Radiofrequency spectroscopy covers the heretofore unused portion 
of the electromagnetic spectrum used for absorption spectroscopy. 
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The measurements involved are based upon the absorption of radio- 
frequency energy by nuclei and/or free electrons when they are simul- 
taneously subjected to both a steady magnetic field and a radio- 
frequency field. The amount of energy absorbed will give a quantita- 
tive analysis. The width and shape of the absorption line will give a 
measure of the chemical and physical environment of the nuclei or 
electrons. The saturation energy will give added information on in- 
ternal environment. These measurements can be made under many 
dynamic conditions so that the chemical and/or physical changes can 
be measured as a result of varying external environment. 

It is therefore possible to apply these radiofrequency spectroscopy 
techniques to quality control and measurement, as well as to research, 
in such fields as high polymers, plastics, metals, petroleum products, 
rubber, paint, and starch. 


Dynamic Mechanical Properties of Concentrated Solutions of Poly-n- 
Butyl Methacrylate 


P. R. SAUNDERS, DONALD M. STERN, SHELDON F. 
KURATH, CHAROEN SAKOONKIM, and JOHN D. FERRY 
Department of Chemistry, University of Wisconsin, Madison, Wisconsin 


Measurements of the real and imaginary parts of the complex 
compliance (J*) have been made on solutions of poly-n-butyl meth- 
acrylate (molecular weight 3.05 x 10°) in diethyl phthalate at con- 
centrations of 30, 40, 50, and 60% by weight, using the Fitzgerald 
transducer apparatus over wide ranges of frequency and temperature. 
Except for certain anomalies mentioned below, the method of reduced 
variables provided single composite curves with temperature shift 
factors which followed the WLF equation; the parameter a, varied 
from 2.2 to 2.8 K 10~‘, and for the 50 and 60% solutions, whose glass 
transition temperatures were determined refractometrically as —67° 
and —46°C., the parameter f, was 0.024. At the lowest temperatures 
and highest frequencies, the reduced data for the 60% solution showed 
an anomaly which was tentatively identified as a 6-mechanism but 
could not be separately resolved. At temperatures and frequencies 
corresponding to the neighborhood of the maximum in J’’, the reduced 
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data for all four solutions showed an anomaly like that in the undiluted 
polymer, where the maximum shifts upward and to lower reduced fre- 
quencies with increasing temperature. This behavior can be inter- 
preted in terms of a temperature dependence of the number of entan- 
glement coupling points, with an apparent heat of dissociation of 2.1 
keal. in all solutions. Retardation and relaxation spectra have been 
calculated for each solution, and from these the concentration depend- 
ences of the monomeric friction coefficient and other parameters have 
been obtained. 


(This work was supported in part by the Ordnance Corps, Department of the 
Army.) 


[To be published in J. Colloid Sci.] 


Some Generalizations of Linear Viscoelastic Stress-Deformation 
Relationships 


HERSHEL MARKOVITZ 
Mellon Institute, Pittsburgh, Pennsylvania 


Oldroyd has developed a method for generalizing rheological equa- 
tioas of state of limited applicability so that they have correct in- 
variance properties and are of a universally valid form. In this paper 
a general Maxwell body is taken as a basis for various generalizations 
of the Oldroyd type. These results are then examined for their im- 
plications with respect to the rate of shear dependence of the shearing 
and normal stresses in steady state flow. 


Dynamic Mechanical and Creep Properties of a Concentrated 
Cellulose Nitrate Solution 


DONALD J. PLAZEK* 
Department of Chemistry, University of Wisconsin, Madison, Wisconsin 


Low frequency torsion pendulum and creep measurements have 
been made on a solution of cellulose nitrate in diethyl phthalate 


* Dr. Plazek’s present address is Mellon Institute, Pittsburgh, Pa. 
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(23.0% by weight of polymer) to supplement previously reported 
work with a Fitzgerald-Ferry double transducer. The torsion pen- 
dulum measurements were made at temperatures from — 25 to +35°C. 
within a frequency range extending from 0.07 to 9.0 cycles/sec. 
Creep runs, also made in torsional shear, were made at temperatures 
from —33 to +65°C. At the temperatures where crystallinity was 
present (5°C. and below), the deformation, 7, at long times can be 
described by the expression y = yo + 6;/*, where ¢ is the time and 
yo and £ are constants. 

In spite of the presence of crystallinity, all of the results were 
independent of thermal history and the dynamic results at — 10°C. 
and above could be and were reduced to a reference temperature 
(T> = 25°C.). The necessary shift factors, ay, are given by log 
ar = —8.48 (T — Ty) /(166 + T — To). At lower temperatures the 
dispersion is appreciably broadened by the crystallinity. 


[To be published in J. Colloid Sci.] 


On the Thermostatics of Continuous Media 


BERNARD D. COLEMAN and WALTER NOLL 
Mellon Institute, Pittsburgh, Pennsylvania 


This paper discusses materials in states of rest; that is, in states 
which are achieved in the limit of very slow measurements. When 
in such a state, every material possesses a strain energy function, and 
behaves as either a perfect fluid or a perfectly elastic solid. A strain 
energy function cannot be completely arbitrary. Its possible forms 
are limited by symmetry considerations (i.e., objectivity and isotropy) 
and by the second law of thermodynamics. The restrictions caused 
by symmetry have been understood for some time; in this paper we 
investigate the thermodynamic requirements. In the course of a 
derivation of the laws of hydrostatics and elastostatics from thermo- 
dynamics, certain inequalities are obtained. These inequalities, 
which express the thermodynamic restrictions on a strain energy 
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function, have important implications when applied to the theory 
of finite elastic strain. 


[To be published in Arch. Rational Mech.] 


Four Separate “Regimes of Flow’’ or Systems of Blood Rheology 


WILLIAM H. DUNHAM, FANN HARDING DUNHAM, and 
MELVIN H. KNISELY 
Department of Anatomy, Medical College of South Carolina, Charleston, 
South Carolina 


Blood flow patterns resulting from different properties of cells, 
combinations of cells, and plasma, during health and disease have not 
been thoroughly analyzed. At least four regimes of flow are known: 

1. Newtonian flow, which consists of a parabolic velocity distribu- 
tion of a homogenous fluid around the axis of flow. This occurs when 
plasma alone flows through a cylindrical vessel. 

2. Poiseuillian flow, the commonly described flow pattern, which 
consists of alternate concentric cylinders of cells suspended in plasma, 
and plasma, radially symmetrical around the axis of the vessel. 

The above two regimes of flow occur during health when blood 
zells are not agglutinated; during disease, blood cells often aggluti- 
nate into masses having different sizes and different degrees of internal 
rigidity. The following two types of fiow sometimes occur during 
disease: 

3. Plug flow, a plastic flow during which the “‘yield value” of the 
masses has been exceeded at their boundaries. 

4. Gravitationally layered flow, as seen with the observer looking 
horizontally, which consists of several layers of settled, stationary, 
and/or sliding cellular constituents, often agglutinated, of different 
sizes, shapes, densities, and elastic properties. And above these 
there is often a layer of moving plasma containing suspended small 
particles such as small masses, free erythrocytes, and leukocytes. 
This flow is not radially symmetrical around the axis of the vessel. 

Between branchings, segments of arteries and of veins are truncated 
cones, not cylinders (Jeffords and Knisely, 1956). Saunders and 
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Knisely (1954) found that the maximal internal diameters to which 
the terminal mesenteric arterioles can distend as they are impacted 
with masses of sludge are: for hamsters and mice, 9 micra; for rats 
and rabbits, 8 micra; for cats, 11 micra; and for dogs, 10 micra. 
These are the narrowest vessels of the vascular system. 

The above four separate regimes of flow are considered in terms 
of the forces of the plasma, influenced by the vessel walls, acting on 
the surfaces of the suspended particles. The paper was illustrated by 
lantern slides and motion pictures taken through a horizontally 
aimed microscope. As far as is known, rheologists have not ordinarily 
used horizontally aimed microscopes to study the flow of suspensions. 
The walls of small living blood vessels are nearly transparent, which 
permits direct observation of the behavior of particles in “flowing” 
suspensions. 


(This research was supported by U.8.P.H.S. Grant H-1683.) 


[Published in Angiology, 9, 317 (1958).] 


Effects of Drawing on the Mechanical Behavior of Nylon 6 Filaments 


J. H. DUSENBURY, C. WU, C. J. DANSIZER, and A. HAUSER 
Textile Research Institute, Princeton, New Jersey 


The extensional properties, bending moduli, and torsional moduli 
have been determined for a series of nylon 6 filaments ranging in draw 
ratio from 1.00 to 5.48. The extensional properties were measured 
with an Instron tester; the bending and torsional moduli were 
measured with an electrostatic vibroscope. Both the extensional 
and bending moduli increase smoothly with increasing draw ratio, 
whereas the torsional modulus is essentially independent of the draw 
ratio. These findings are in good agreement with those reported 
previously for two other melt-spun filaments: nylon 66 and Dacron. 
The actual values observed for the nylon 6 moduli, however, are about 
2/; the corresponding values observed in the case of nylon 66. During 
load-extension experiments the filaments of lower draw ratio were 
observed to go through relatively long periods of extension without 
appreciable changes in load; they were, in effect, being cold-drawn 
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under the test conditions. The values observed for the ratio of the 
bending modulus to three times the torsional modulus suggest that 
there is no preferential orientation in the undrawn filaments and 
that these filaments exhibit increasing molecular orientation with 
respect to the filament axis as the extent of drawing increases. These 
views are in general agreement with the results of x-ray diffraction 
studies carried out on the same materials. 


The ‘Effective’ Solid Volume—Viscosity Relationship of Rigid 
Spheres Versus the Dimensions of the Measuring Instrument 


STEWART R. OGILBY 
Research Division, Gocdyear Tire & Rubber Co., Akron, Ohio 


A viscosity—solids relationship that is based upon the effect of the 
solid phase in altering the “effective’’ dimensions of the measuring 
instrument is proposed. 

The “effective’’ solid volume fraction has been found to vary with 
the ratio between the shearing area of the instrument and the diameter 
of the solid particles relative to the volume of sample under shear. 
The absorption or interaction forces at the solid-liquid interface affect 
the “effective” solid-volume fraction V, in direct relation to the 
specific surface of the solid phase and the amount of shear. 

Two simple equations are used to relate the solid-volume fraction 
(V) to the “effective” solid—volume fraction over a range of shear that 
shows both pseudo-plastic and dilatant flow. 


A Method of Generalizing the Effect of Temperature and Rate on the 
Stress-Strain Curves of Nylon 


JOHN B. MILES 
Chemstrand Corporation, Decatur, Alabama 


It has been found possible to considerably simplify the description 
of the stress-strain curve of nylon as a function of temperature and 
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of the rate of elongation. The stress-strain curves obtained at a con- 
stant rate of elongation in the range of temperature from 21 to 150°C. 
can be accurately represented by one master curve with the units along 
the stress axis multiplied by a constant which depends on the te . per- 
ature. Thus the stress P is given by 


P = f(e) X 9(T) 


where e and 7’ are the strain and temperature, respectively, and f(e) is 
a function of e only, and g(7) of T only. The temperature function, 
g(T), has been determined over this temperature range. 

A similar type of relationship has been found for different rates 
of elongation at constant temperature in the range from 0.2 to 200% 
elongation per minute. 

The characteristics of the master curve, f(e), are discussed, with 
possible interpretation in terms of the mechanism of deformation. 


[To be published in Textile Research J.] 














